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INTRODUCTION 


A computer program based on state-of-the-art compressor and 
structural technologies applied to beaded shrouded discs has been 
developed and made operational In NASTRAN Level 16. 

The problems encompassed Include aeroelastlc analyses, modes 
and flutter. 

The program Is documented In the form of five NASA Contractor's 
Reports — — one Technical Report and four Updates to NASTRAN Level 
16 Theoretical, User's, Programmer's and Demonstration manuals. 

This report describes the User's manual updates. 


USER'S MANUAL UPDATES 


APPENDIX 1 -a USER'S MANUAL SUPPLEMENT 




TABLE OF CONTENTS (Continued) 

Section Page No. 

1.11 AEROELASTIC MODELING 1.11-1 

1.11.1 Introduction 1.11-1 

1.11.2 Aerodynamic Modeling 1.11-1 

1.11.3 The Interconnection Between Structure and 

Aerodynamic Models 1.11-3 

1.11.4 Modal Flutter Analysis 1.11-5 

1.11.5 Sample Problem.... 1.11-6 

1.12 CYCLIC SYMMETRY 1.12-1 

1.13 FULLY STRESSED DESIGN 1.13-1 

1.14 AUTOMATED MULTI-STAGE SUBSTRUCTURING 1.14-1 

1.14.1 Substructuring Terminology 1.14-2 

1.14.1.1 Storage of Substructure Data 1.14-2 

1.14.1.2 Identification of Substructure Data 1,14-7 

1.14.1.3 Input Data Checking 1.14-8 

1.14.2 The Substructure Operating File ($0F) 1.14-10 

1.14.3 The Case Control Deck for Substructure Analyses 1.14-11 

1 .14.3.1 Phase 1 1.14-12 

1.14.3.2 Phase 2 1.14-12 

1.14.3.3 Phase 3 1.14-12 

1.14.4 Example of Substructure Analysis 1.14-13 

1.15 STATIC AEROELASTIC AND FLUTTER MODELING 

OF AXIAL FLOW TURBOMACHINES 1.15-1 

1.15.1 Introduction 1.15-1 

1.15.2 Aerodynamic Modeling 1.15- 

1.15.3 Steady Aero thermoel as t i c ’Des 1 gn/Ana 1 ys i s “ 1.15- 

1.15.3.1 Aerodynamic DTI Data 1.15- 

1.15.3.2 Aerodynamic Output Data 1.15- 

1.15.4 Sample Problem 1.15- 

1.15.5 Modal, Flutter and Subcritical Roots Analyses 1.15- 

1.15.6 Sample Problem . 1.15- 


\ 


v (9/30/78) 


TABLE OF CONTENTS (Continued) 


Section Page No . 

2. NASTRAN DATA DECK 

2.1 GENERAL DESCRIPTION OF DATA OECK 2.1-1 

2.2 EXECUTIVE CONTROL OECK 2.2-1 

2.2.1 Executive Control Card Descriptions 2.2-1 

2.2.2 Executive Control Deck Examples 2.2-5 

2.3 CASE CONTROL DECK 2.3-1 

2.3.1 Data Selection 2.3-1 

2.3.2 Output Selection 2.3-2 

2.3.3 Subcase Definition 2.3-3 

2.3.4 Case Control Card Descriptions 2.3-6 


va (9/30/78) 


TABU ?• CONTENTS (Continued) 


Igir'^rrr '• T<-- „• • 


r. 


Section Page No , 

3.18.2 Description of DMAP Operations for Nonlinear 

Static Heat Transfer Analysis 3,18-5 

3.18.3 Case Control Deck and Parameters for Nonlinear 

Static Meat Transfer Analysis 3,18-8 

3.19 TRANSIENT HEAT TRANSFER ANLAYSIS... 3.19,1 

3.19.1 DMAP Sequence for Transient Heat Transfer 

Ana lys Is 3.19-1 

3.19.2 Description of DMAP Operations for Transient 

Heat Transfer Analysis 3.19-8 

3.19.3 Case Control Deck and Parameters for Transient 

Heat Transfer Analysis 3,19-13 

3.20 MODAL FLUTTER ANALYSIS 3.20-1 

3.20.1 DMAP Sequence for Modal Flutter Analysis 3.20-1 

3 20.2 Description of DMAP Operations for Modal 

Flutter Ana lysis 3 . 20-9 

3.20.3 Output for Modal Flutter Analysis 3.20-16 

3.20.4 Case Control Deck and Parameters for Modal 

Flutter Analysis 3.20-16 

3.20.5 Modal Flutter Analysis Subsets 3.20-18 

3.20.6 DMAP Sequence for Modal Flutter 

Analysis, Subset 4 3.20-19 

3.20.7 DMAP Sequence for Modal Flutter 

Analysis, Subset 5 3.20-21 

3.21 COMPRESSOR BLADE MESH GENERATOR 3.21-1 

3.21.1 DMAP Sequence for Compressor Blade Mesh 

Genera tor 3 . 2 1 - 1 

3.21.2 Description of DMAP Operations for Compressor 

Blade Mesh Generator 3.21-2 

3.21.3 Output for Compressor Blade Mesh Generator 3.21-3 

3.21.4 Case Control Deck, DTI Table and Parameters 

for Compressor Blade Mesh Generator 3.21-4 

3.22 STATIC AE ROTH ERMOELAST I C ANALYSIS WITH DIFFERENTIAL 

STIFFNESS 3.22-1 

3.22.1 DMAP Sequence for Static Aero th >rmoe 1 a s t i c 

Analysis 3.22-1 

3.22.2 Description of DMAP Operations for Static 

Aerothermoelastic Analysis 3.22-13 


x (9/30/78) 


TABLE OF CONTENTS (Continued) 

Section Peoe No. 

3. 2?. 3 Output for Static Aerothermoelastlc 

Analysis 3.22-22 

3.22.4 Case Control Deck, DTI Table and Parameters 

for Static Aerothermoelastlc Analysis 3.22- 

3.23 COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 3.23-1 

3.23.1 DMAP Sequence for Blade Modal F .ter 

Analysis 3.23-1 

3.23.2 Description of DMAP Operations for Blade 

Modal Flutter Analysis 3.23-H 

3.23.3 Output for Blade Modal Flutter Analysis 3.23-19 

3.23.4 Case Control Deck and Parameters for 

Blade Modal Flutter Analysis 3.23-19 

4. PLOTTING 

4.1 PLOTTING 4.1-1 

4.2 STRUCTURE PLOTTING 4.2-1 

4.2.1 General Rules 4.2-2 

4. 2. 1.1 Rules for Free-Field Card 

Specifications 4.2-2 

4. 2. 1.2 Plot Request Packet Card Format 4.2-2 

4. 2. 1.3 Plot Titles 4.2-2 

4.2.2 Plot Request Packet Card Descriptions 4.2-2 a 

4.2.2. 1 SET Definition Cards 4.2-3 

4. 2. 2. 2 Cards Defining Parameters 4.2-4 

4.2.3 Summary of Structure Plot Request 

Packet Cards 4.2-16 

4.3 X-Y OUTPUT 4.3-1 

4.3.1 X-Y Plotter Terminology 4.3-1 


xa (9/30/78) 


TAB it OF CONVENTS (Continued) 


Secti on Page No, 

6. DIAGNOSTIC MESSAGES 

6.1. RIGID FORMAT DIAGNOSTIC MESSAGES 6,1-1 

6.1.1 Displacement Approach Rigid Formats. 6.1-1 

6. 1.1.1 Rigid Format Error Messages for Static 

Analysis 6,1-1 

6. 1.1. 2 Rigid Format Error Messages for Static 

Analysis with Inertia Relief 6,1-1 

6. 1.1. 3 Rigid Format Error Messages for Normal 

Mode Ana lys 1 s 6.1-? 

6. 1.1. 4 Rigid Format Error Messages for Static 
Analysis v/lth Differential Stiffness.... 6,1-2 

6. 1.1. 5 Rigid Format Error Messages for Buckling 

Analysis 6.1-3 

6. 1.1. 6 Rigid Format Error Messages fo*- Piecewise 

L 1 near Ana lysis... 6.1-3 

6. 1.1. 7 Rigid Format Error ftessages for Direct 

Complex Eigenvalue Analysis 6.1-4 

6. 1.1. 8 Rigid Format Error Messages for Direct 

Frequency and Random Response 6.1-4 

6. 1.1. 9 Rigid Format Error Messages for Direct 

Transient Response 6.1-4 

6.1.1.10 Rigid Format Error Messages for Modal 

Complex Eigenvalue Analysis 6,1-6 

6.1.1.11 Rigid Format Error Messages for Modal 

Frequency and Random Response 6.1-5 

6.1.1.12 Rigid Format Error Messages for Modal 

Transient Response 6.1-5 

6.1.1.13 Rigid Format Error Messages for Normal 

Modes v/lth Differential Stiffness 6.1-6 

6.1.1.14 Rigid Format Error Messages for Static 

Analysis Using Cyclic Symmetry 6,1-6 

6.1.1.15 Rigid Format Error Messages for Normal 

Modes Analysis Using Cyclic Symmetry..., 6.1-7 

6.1.1.16 Rigid Format Error Messages for Static 

Aerothermoelast ic Analysis with 
Differential Stiffness 6,l-7a 

6 . 1 * T Heat Approach Rigid Formats 6.1-8 

preceding page blank not filmed 


x 1 1 1 (9/30/78) 


TABIC OF CONTENTS (Continued) 


Section Pipe No, 

Rigid Format Error Messages for Static 

Heat Transfer Analysis 6.1-8 

6.1. ?.? Rigid Format Error Messages for Nonlinear 

Static Heat Transfer Analysis........... 6.1-8 

6.1, ?. 3 Rigid Format Error Messages for Transient 

Meat Transfer Analysis 6.1-9 

6.1.3 Aero Approach Rigid Format 6.1-9 

6. 1.3.1 Rigid Format Error Messages for Modal 

F 1 utter Ana lysis 6.1-9 

6. 1.3. 2 Rigid Format Error Messages for 

Compressor Clade Cyclic Modal Flutter 
Analysis 6.1-9 


xiiia (9/30/78) 


USER'S MANUAL UPDATES 


STRUCTURA L MODELING 

1.15 STATIC AEROELASTIC AND FLUTTER MODELING OF AXIAL FLOW 

TURBOMACHINES 
1.15,1 Introduction 

The NASTRAN aeroelastic and flutter capability has 
been extended tc solve a class of problems associated with axial 
flow turbomachines . The capabilities included are; 

1. Steady state aerothermoelastic analysis of 
compressors to determines 

(a) The change in geometry between the design 
point operating shape and the "as manufactured" shaoe of the 
flexible blade to ensure the required performance (pressure 
ratio, flow rate r rpm) at the design point. (This is termed 
the "design" problem.) 

(b) The performance at off-design operating 
conditions for a given "as manufactured" blade shape. (Thi3 is 
termed the "analysis" problem.) 

(c) Displacements, stresses, reactions, plots, 
etc., at selected operating points over the compressor map. 

(d) A differential stiffness matrix due to 
centrifugal and aerodynamic pressure and thermal loads for use 
in subsequent modal analysis. 

2. Modal, unstalled flutter and subcritical roots 
analysis of compressors and turbines. 


The rotor/stator of a single-stage, or each stage 
of a multi-stage compressor or turbine is analyzed as an 
isolated structure. Two new Rigid Formats (Displacement 
RF 16 and Aero RF 9) have been developed, one each for the 
aeroelastic steady state and the oscillatory state problems 
(see Sections 3.21, 3.22, 3.23). The rotational cyclic 
symmetry (see Section 1.12) inherent in these structures about 
♦-he axis of rotation has been taken into account in designing 
the capability, so that only a representative one-blade sector 
need be idealized. 

The steady aerothermoelast ic analysis is based on 
the theory described in Volume I of Reference 1. The computer 
code of the same reference (Volume II), with minor changes, 
has been adapted for NASTRAN in the functional module ALG. 

The current NASTRk Static Analysis with Differential Stiffness 
Rigid Format has been accordingly modified to include the effect 
of centrifugal, aerodynamic pressure and temperature loads. 

The existing features of NASTRAN for Normal Modes 
Analysis using Cyclic Symmetry (Section 3.16) and Modal Flutter 
Analysis (Section 3.20) have been suitably combined for the 
modal flutter and subcritical roots analysis of the axial flow 
turbomachinery rotor/stator. 

These developments are compatible with the general 
structural capability in NASTRAN. The structural part of the 
problem is modeled as described in Section 1 of the User's 
Manual. This section deals with the aerodynamic data pertaining 
to the bladcd disc sector. The associated aerodynamic modeling 
is discussed in Section 1.15.2. 


Se.r ion 1.15,3 describes the steady aerothermo- 
elastic "desikjn/ai.alysis" formulations. 

Section 1.15.8 presents the modal, flutter and 
•ubcritical roots analyses. 

Sample problems and their solutions are presented 
in Sections 1.15.4 and 1.15.6. 
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1.15.2 Aerodynamic Modeling 

The aerodynamic model is based on a grid generated 
by the intersection of a series of streamlines and "computing 
stations" (similar to potential lines) as shewn in Figure 1. 

This arrangement also facilitates the subsequent use of two- 
dimensional , unsteady, subsonic and supersonic infinite cascade 
theories (see Section of the Theoretical Manual) in the 
flutter problem. They are used in a strip-theory manner on 
the various streamlines spanning the blade. 

The aerodynamic loads are assumed significant only 
on the blad-d portion of a bladed disc and no other part of 
the structure need be modeled aerodynamically . The data required 
to generate the aerodynamic model for the steady state aeroelastic 
analyses are specified on DTI bulk data cards, and are described 
in Section 1.15.3.1 of the User's Manual. Blade streamline 
data for flutter and subcritical roots analyses are specified 
on STREAMLi bulk data cards. 

The streamlines are defined by the intersection of 
the blade mean surface and a set of coaxial cylindrical (or 
conical) surfaces. The axis of the cylinders (cones) coincides 
with the axis of rotation of the turbomachine. The "computing 
stations" lie on the blade mean surface and divide it from 
the leading edge to the trailing edge. The choice of the number 
and location of the streamlines and the "computing stations" is 
dictated by the expected variation of the relative flov/ 
properties across the blade span, and the complexity of the 
mode shapes exhibited by this part of the structure. However, 
a minimum of three streamlines (including the blade root and 




the tip) and three "computing stations" (including the blade 
leading edge and the trailing edge) must be specified. 

The distribution of the aerodynamic parameters over 
the blade is, in general, different from that of the structural 
parameters such as stress, strain, etc. Accordingly, the 
aerodynamic model and the structural model of the blade, in 
general, may differ. The difference currently permitted in 
the two models is as seen in Figure 1 wherein the aerodynamic 
grid is shown to be a part of the structural grid. 

The x-axis of the BASIC coordinate system (Figure l)is chosen 
to coincide with the axis of rotation and is oriented in the 
direction of the flow. The location of th origin is arbitrary. 

The z-plane (BASIC) lies normal to the "mean" meridional plane 
passing through the blade, with the z-axis (BASIC) directed 
towards the blade. The aerodynamic grid can be specified in 
any coordinate system (CP) . The aerodynamic model data mainly 
related to the bladad disc problems are specified on the DTI, 

STREAML1 and STREAML2 bulk data cards. 
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1.15.3 


Steady Aero thermoe la Stic "Design/Analysis " 

An operating point on a compressor map defines a 
distribution of centrifugal force and aerodynamic pressure and 
temperature loads on the biaded-disc of the axial flow turbo- 
machine. The equilibrium, deformed shape of the elastic 
structure is reached at the end of a series of quasi-equilibrium 
states during which the loads on the bladed-disc and its geometric 
stiffness change as a function of the deformation. The operating 
point pressure ratio (given the flow rate and the rpm) , in effect, 
also changes during this process. 

Two different problems can thus be stated: 

1. Given the desired design operating point and 

the "rigid” geometry, to determine the "as manufactured" geometry 
("design" problem) that would produce the design conditions and 

2. Given the "as manufactured" geometry, to deter- 
mine the performance of the flexible blade at off-design operating 
points ("analysis" problem) . 

Rigid format Displacement 16 has been developed to 
solve these "design/analysis" problems. The value of the 
PARAMeter SIGN (* +1) selects the analysis or the design mode 
of the rigid format. Deformation of the structure as a result 
of the applied centrifugal and aerodynamic loads is used to 
revise the blade geometry each time through the differential 
stiffness loop of the rigid format. Because of the non-linear 
relationship between the blade geometry and the resulting 
operating point pressure ratio, provision is made to control 
the fraction of the displacements used to redefine the blade 
geometry. Thi.s is especially helpful in the solution of the 
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"design" problem. The fractions of the displacements used to 
redefine the blade geometry are specified via the FXCOOR, FYCOOR 
and FZCOOR parameters. The application of the aerodynamic 
pressure and thermal loads is controlled respectively by the 
parameters APRESS and ATEMP. These parameters also enable the 
inclusion of the centrifugal loads alone. 

The functional module ALG is U3ed in the rigid format 
before, within and after the differential stiffness loops (see 
Section ) to generate the aerodynamic loads. Printed output 

from this module during these three stages can respectively be 
controlled through the use of the parameters IPRTCI, IPRTCL and 
IPRTCF. This enables observation of the variation in the 
aerodynamic loads as a function of the blade geometry. 

GRID, CTRIA2 and PTRIA2 bulk data cards for the 
final blade shape can be punched out using the parameter PGE0M . 
At the end of a "design" run, these define the "as manufactured" 
blade shape which can subsequently be "analyzed" at selected 
operating points over the compressor map. In an 'analysis* 
run at any operating point, the total stiffness (elastic and 
geometric) of the bladed-disc structure can be saved via the 
parameter KT0UT for use in subsequent modal, modal flutter 
and subcritical roots analyses. 

The subsections 1.15.3.1 and 1.15.3.2 describe 
the aerodynamic Direct Table Input and the output data for 
the steady state analyses. 


1.15.3.1 Aerodynamic DTI Data 

The input data consist of an initial indication 
of the number of entries that are to be made to each of the 
two program sections (analytic meanline blade section and 
aerodynamic section), and then a data-set for each entry to 
each section. The data that are required for the interfacing 
of the output from the analytic meanline blade section to 
the aerodynamic section are included in the data-set for the 
analytic meanline section. Because partial input to the aero- 
dynamic section is generated by execution of the analytic mean- 
line section, the input for the aerodynamic section to be 
supplied directly by the user varies. This is indicated in 
the charts below by giving the variable name L0G5 for the 
file from which any data are taken that are not always 
supplied directly. 

L0G5 is the file from which input is taken that is 
generated by the analytic meanline section. When the analytic 
meanline section has been directed to produce data for the 
aerodynamic section for a particular computing station, L0G5 
becomes an internally generated scratchfile. Otherwise, L0G5 
is attached to the standard input unit and the user supplies 
the data. 

The following input data items must be input using 
NASTRAN Direct Table Input (DTI) bulk data cards. A description 
of the DTI card is in the NASTRAN User's Manual on page 2.4-105. 
The table data block name must be ALGDB. The trailer value 
for T1 is the number of logical records in the DTI table, not 
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counting the header record. This is the same as the maximum 
value of IREC used in the table. The trailer values for T2 
through T6 are all zero. Each of the following input cards 
corresponds to one logical record of the DTI table. 

Trailing zeroes need not be input. Data types, i.e., alpha- 
numeric (BCD), real and integer, must correspond to those 
specified for each data item. Data item names that begin with 
the letters I,J,K,L,M, and N are to be input as integers while 
all others are input as real numbers. Titles are input as 
alphanumeric (BCD) with the restriction that only alphabetic 
letters occupy the first character in each field of the DTI 
card. Titles may use up to nine DTI fields. 


In the following chart, one line, which may be 
continued , corresponds to one logical record of a DTI 
card. 


TITLEl 
NANAL NAERO 

The following data-set is input to the analytic 
mean line section, and will occur NANAL tiroes. 

The last record in this set is indicated with an 
asterisk. 


TITLE2 

NLINES NSTNS NZ NSPEC NPOINT NBLADE ISTAK 
(cont.) IPUNCH ISECN IFCORD IFPLOT IPRINT ISPLIT INAST 
(cont.) I RLE IRTE NSIGN 


Z INNER ZOUTER SCALE STACKX PLTSZE 


KPTS I FANGS 


XSTA RSTA 
R BLAFOR 


- Occurs KPTS times 

- Occurs NLINES times 


Occurs 

NSTNS 

Tiroes 


ZR B1 B2 PP QQ RLE 
TC TE Z CORD DELX DELY 


Occurs 

NSPEC 

Times 


S BS - Only if ISECN = 1 or 3 
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NRAD NDFTS NDATR NSWITC NLE NTE 


XKSHPE SPEED 


N0UT1 NOUT2 NOUT3 - Refers to leading 

edge station 


(cont. ) 
(cont. ) 


NR NTERP NMACH NLOSS NL1 
NL2 NEVAL NCURVE NLITER NDEL 
NOUTl NOUT2 NOUT3 NBLAD 
R XLOSS - Occurs NR times 


RTE 


DM DVFRAC ] -Occurs NDPTS 


1 

tunes 


Occurs 
for each 
station 
within 
blade or 
at trailing 
edge 

Occurs 

NRAD 

times 


RDTE DELTAD AC-Occurs NDATR times 


This group 
is used to 
generate 
LOG5 data 
for the 
aerodynamic 
section 


The following data-set is input to the aerodynamic 
section and the last record in this set is indicated 
with a double asterisk. 


TITLE3 


CP GASR G EJ 


NSTNS NSTRMS NMAX NFORCE NBL NCASE 
(cont.) NSPLIT NSETl NSET2 NREAD NPUNCH NPLOT 
(cont.) NPAGE NTRANS NMIX NMANY NSTPLT NEQH NLE NTE NSIGN 
NWHICH - Occurs NMANY times on the same card 

G EJ SCLFAC TOLNCE VISK SHAPE 

XSCALE PSCALE RLOW PLOW XMMAX RCONST 
CONTR CONMX 


FLOW SPDFAC 
NSPEC 

XSTN RSTN - Occurs NSPEC times 


Occurs 

NSTNS 

Times 
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(LOG 5) 

(cont. ) 

(cont . ) 

(LOGS) 

(LOGS) 

(cont. ) 

(LOGS) 


NDATA NTERP NDIMEN NMACH 

DATAC DATA1 DATA 2 DATA 3 - Occur* 

NDATA times 

NDATA NTERP NDIMEN NMACH NWORK 

NLOSE NL1 NL2 NEVAL NCURVE NLITER 

NDEL NOUT1 NOUT2 NOUT3 NBLADE 
SPEED- If NDATA >0 
DATAC DATAl DATA 2 DATA 3 DATA4 

DATA 5 

DATA 6 DATA7 DATA 8 DATA 9 _ 

DELC DELTA - Occurs NDEL times 


Occurs 

NDATA 

times 


Inlet 

condition 

specification 


For 

sta- 

tions 

2 

thru 

NSTNS 


WBLOCK BBLOCK BDIST -Occurs NSTNS times 


NDIF 

DIFF FDHUB FDMID FDTIP -Occurs NDIFF 

times 



Occurs 

NSET1 

times 


NM NRAD 


occurs 


TERAD I Occurs 

NRAD 

DM WFRAC -Occurs NM times I times 


NSET2 

Times 


DELF(l) DELF (2) DELF (NSTRMS) - if NSPLIT = 1 

or NREAD = 1 


** R X XL II JJ - Occurs NSTRMS times for NSTNS stations 

if NREAD = 1 
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Data Item Defintionr'. 

The aerodynamic section may be usod with any self- 
consistent unit system and, additionally, a "linear dimension 
scaling factor" (SCLFAC) is incorporated into the input so that 
some commonly used but inconsistent unit sytems may be used. 

This is principally intended to allow the use of inches for 
physical dimensions and yet retain feet'-for velocities. The 
basic dimensions used in the data are length (L) , time (T) , and 
force (F) . Angles are expressed in degrees (A), and temperatures on 
an absolute temperature scale (D) . Heat capacities (H) are also 
required. Some possible unit systems are given below, togehter with 
the corresponding value of SCLFAC , 


L 

T 

F 

D 

H 

SCLFAC 

Feet 

Seconds 

Pounds 

Deg. Rankine 

BTU 

1.0 

Inches 

Seconds 

Pounds 

Deg. Rankine 

BTU 

12.0 

Meters 

Seconds 

Kilograms 

Deg. Kelvin 

CHU 

1.0 


Note that 

some data names 

are used in more 

than 

one 


section; care should betaken to consult the correct sub-division 
below for defintions. 

a. Initial Directives 

TITLEl Thi^: is a title card for the run. 

NANAL Set NANAL = 1 

NAERO Set NAERO » 1 

b. Analytic Meanline Blade Section 

For a more detailed discussion of the input to 
this section through item XB, see Reference and . For 
this section, the dimensioned input is either in degree (A) 
or in length (L) . 
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TITLE2 


NLINES 


NSTNS 


NZ 


NSPEC 


NPOINT 


NBLADE 

ISTAK 


A title card for the analytic meanline 

section of the program. 

The number of stream surfaces which are 
defined, and on which blade sections 
will be designed. Must satisfy 

2 < NLINES s 21. 

The number of computing stations at which 
the stream surface radii are specified. 

Must satisfy NSTNS s 10. 

The number of constant-z planes on which 
manufacturing (Cartesian) coordinates 
for the blade are required. Must satisfy 

3 5 NZ < 15. 

The number of radially disposed points at 
which the parameters of the blade sections 
are specified. Must satisfy is NSPEC *• 21. 
The number of points that will be generated 
to specify the pressure and suction surfaces 
of each blade section. Must satisfy 
2< NPOINT 5 80. Generally, no lest* than 30 
should be used. 

The number of blades in the blade row. 

If ISTAK ■ 0, the blade will be stacked at 
the leading edge. 

If ISTAK = 1, the blade will be stacked at 
the trailing edge. 
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1 PUNCH 
ISECN 


IFCORD 


IFPLOT 

IPRINT 


If ISTAK » 2, the blade will be stacked at, 
or offset from, the section centroid. 

Set i punch - o 

If ISECN - 0, the blade will be constructed 
using the polynomial camber line and the 
standard (i.e., double-cubic) thickness 
distribution. 

If ISECN ■ 1, the exponential camber 

line and the standard thickness distribution 

will be used. 

If ISECN * 2, the circular arc c a r line 

and the double-circular-arc thickness 
distribution will be used. 

If ISECN ■ 3, the multiple-circular-arc 
meanline and tne standard thickness 
distribution will be used. 

If IFCORD = 0, the meridional projection 
of the stream surface blade section chords 
are specified. 

If IFCORD * 1, the stream surface blade 
section chords are specified. 

Set IFPLOT = 0 

The input data is always listed by the 
program. Details of the stream surface 
and manufacturing sections are printed 
as prescribed by IPRINT. 
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ISPLIT 

INAST 

I RLE 
IRTF 
NSIGN 


If IPRINT ■ 0, details* of the stream 
surface and manufacturing sections are 
printed. 

If IPRINT ■ 1, details of stream surface 
sections are printed. 

If IPRIFT ■ 2, details of manufacturing 
sections are printed. 

If IPRINT - 3, details of neither stream 
surface nor manufacturing sections are 
printed. (The interface data for use with 
the aerodynamic section of the program is 
is still displayed.) 

Set ISPLIT * 0 

Set INAST * 0. See the Output Data 
description (Section ) for further 
details. 

The computing station number at the blade 
leading edge. 

The computing station number at the blade 
trailing edge. 

Indicator used to sign blade pressure 
forces according to program sign conven- 
tions. For compressor rotors , if the 
machine rotates clockwise when viewed 
from the front, set NSIGN to 1; other- 
wise, set NSIGN to -1. For compressor 
stators , the two values given for NSIGN 
are reversed. 
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Z INNER, 
ZOUTER 

SCALE 

STACKX 


PLTSZE 

KPTS 


IFANGS 


XSTA 


The NZ manufacturing sections are equi- 
spaced between z equals ZINNER and ZOUTER. 
Set scale * 0.0. 

This is the axial coordinate of the stacking 
axis for the blade, relative to the same 
origin as used for the station locations, 
XSTA. 

Set PLTSZE » 0.0, 

The number of points provided to specify 
the shape of a computing station. 

If KPTS = 1, the computing station is 
upright and linear. 

If KPTS = 2, the computing station is 
linear and either upright or inclined. 

If KPTS >2, a spline curve is fit through 
the points provided to specify the shape 
of the station. 

If IFANGS = 0, the calculations of the 
quantities required for aerodynamic 
analysis will be omitted at a particular 
computing station. 

If IFANGS = 1, these calculations will 
be performed at that station. 

An array of KPTS axial coordinates (relati/e 
to an arbitrary origin) which, toghter with 
RSTA, specify the shape of a particular 
computing station. 
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RSTA 


R 

BLAFOR 

ZR 


B1 

B2 

PP 


An array of KPTS radii which , together 
with X5TA, specify the shape of a 
particular computing station. 

The stream surface radii at NLINES locations 
at each of the NSTNS stations. 

Set BLAFOR = 0.0. 

The variation of properties of the stream 
surface blade section is specified as a 
function of stream surface number. The 
various quantitites are then interpolated 
(or extrapolated) at each stream surface. 

The stream surfaces are numbered con- 
secutively from the inner-most outward, 
starting with 1.0. ZR must increase 
monotonically , there being NSPEC values 
in all. 

The blade inlet angle. 

The blade outlet angle. 

If ISECN = 0, PP is the ratio of the second 
derivative of the camber line at the leading 
edge to its maximum value. Must satisfy 
-2 . 0 < PP< 1.0. 

If ISECN si, PP is the r-'tio of the second 
derivative of the camber line at the 
leading edge to its maximum value forward 
of the inflection point. Must satisfy 
0 . 0 < PP *: 1.0. 
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If ISECN ■ 2 or 3, PP is superfluous. 

QQ If ISECN ■ 0, OQ is the ratio of the second 

derivative of the camber line at the trailing 
edge to its maximum value. Must satisfy 
0 . 0 s QQ s 1 . 0 . 

If ISECN ■ 1, QQ is the ratio of the second 
derivative of the camber line at the trailing 
edge to its maximum value rearward of the 
inflection point. Must satisfy 0.0<QQ.<1.0. 
If ISECN = 2 or 3, QQ is superfluous. 

RLE The ratio of blade leading edge radius to 

chord . 

TC The ratio of blade maximum thickness to 

chord. 

TE The ratio of blade trailing edge half- 

thickness to chord. 

If ISECN s 2, TE is superfluous. 

Z The location of the blade maximum thickness , 

as a fraction of camber line length 
from the leading edge. 

If ISECN = 2, Z is superfluous. 

CORD If IFCORD - 0, CORD is the meridional 

projection of the blade chord. 

If IFCORD = 1, CORD is the blade chord. 

DELX, The stacking axis passes through the stream 

DELY 

surface blade sections, offset from the 
centroids, leading, or trailing edge by DELX 
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S, BS 


NRAD 

NDPTS 

NDATR 

NSWITC 


NLE 

NTE 

XKSHPE 


and DELY in the x and y directions 
respectively. 

If 1SECN * 1 or 3, S and BS are used to 
specify the locations of the inflection 
point (as a fraction of the meridionally- 
projected chord length) and the change in 
camber angle from the leading edge to the 
inflection point. If the absolute value 
of the angle at the inflection point is 
larger than the absolute value of Bl, 

BS should have the same sign as Bl , 
otherwise, Bl and BS should be of opposite 
signs . 

The number of radii at which a distribution of the 
fraction of trailing edge deviation is input. Must 
satisfy 1 £ NRAD £ 5. 

The number of points used to define each deviation curve. 
Must satisfy 1 £ NDPTS £ 11. 

The number of radii at which an additional deviation angle 
increment and the point of maximum camber are 
specified. Must satisfy 1 £ NDATR £ 21. 

If NSWITC = 1, the deviation correlation parameter "m" 
for the NACA (Ajq) meanline is used. 

If NSWITC = 2, the deviation correlation parameter "m" 
for double-circular-arc blades is used. 

Station number at leading edge. 

Station number at trailing edge. 

The blade shape correction factor in the deviation rule. 
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SPEED 


See definition for Aerodynamic Section. 

The number of radii where a "lose" is input. 


NR 

NTERP 

NMACH 

NLOSS 

NL1 

NL2 

NEVAL 

NCURVE 

NL1TER 

NDEL 

NOUTl 

NOUT2 

NOUT3 

NBLAD 

R 


See definition for Aerodynamic Section. 


Radius at which loss is specified. 


XLOSS 

RTE 


DM 


DVFRAC 


Loss description. The form is prescribed by NLOSS; 
see aerodynamic section. 

Radius at blade trailing edge where the following deviation 
fraction/ chord curve applies. 

If NRAD = 1, it has no significance. Must increase 
monotonically. 

The location on the meridional chord where the deviation 
fraction is given. Expressed as a fraction of the 
meridional chord from the leading edge. Must increase 
monotonically. 

Fraction of trailing -edge deviation that occurs at location 
DM. 


RDTE 


Radius at trailing edge where additional deviation and 
point of maximum camber are specified. 
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DELTAD 

AC 

c. 

TITLE3 

CP 

GASR 

G 

EJ 

NSTNS 

NSTRMS 

NMAX 

NFORCE 

NBL 


Additional deviation angle added to that determined by 
deviation rule. Input positive for conventionally positive 
deviation for both rotors and stators. 

Fraction of blade chord from leading edge where maximum 
camber occurs. 

Aerodynamic Section 

A title card for the aerodynamic section of the program. 

Specific heat at constant pressure. An input value of 
zero will be reset to 0. 24. Units; H/F/D. 

Gas constant. An input value of zero will be reset to 
53.32. Units: L/SCLFAC/D. 

Acceleration due to gravity. An input value of zero will 
be reset to 32. 174. Units; L/SCLFAC/ T/ T. 

Joules equivalent. An input value of zero will be reset 
to 778. 16. Units: LF/SCLFAC/H. 

Number of computing stations. Must satisfy 3 £ NSTNS 
s 30. 

Number of streamlines. Must satisfy 3 s NSTRMS s;21. 

An input value of zero will be reset to 11. 

Maximum number of passes through the iterative stream- 
line determination procedure. An input value of zero will 
be reset to 40. 

The first NFORCE passes are performed with arbitrary 
numbers inserted should any calculation produce 
impossible values. Thereafter, execution will cease, 
the calculation having "failed". An input value of zero 
will be reset to 10. 

If NBL = 0, the annulus wall boundary layer blockage 
allowance will be held at the values prescribed by 
WBLOCK. 

If NBL = 1, blockage due to annulus wall boundary layers 
will be recalculated except at station I. VISK and 
SHAPE arc used in the calculation. 
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NCASE 


Set NCASE « 1. 


NSPLIT 


NSET1 


NSET2 


NREAD 


NPUNCH 

NPLOT 

NPAGE 


NTRANS 


If NSPLIT = 0, the flow distribution between the stream- 
lines will be determined by the program so that roughly 
uniform increments of computing station will occur 
between the streamlines at station 1. 

If NSPLIT = 1, the flow distribution between the stream- 
lines is read in (see DELF). 

The blade loss coefficient re-evaluation option (specified 
by NEVAL) requires loss parameter/diffusion factor 
data. NSET1 sets of data are input, the set numbers being 
allocated according to the order in which they are input. 

Up to 4 sets may be input (see NDIFF). 

When NLOSS - 4, the loss coefficients at th'i station are 
determined as a fraction of the value at the trailing edge. 
Then, NSET2 sets of curves are input to define this 
fraction at a function of radius and meridional chord. Up 
to 2 sets may be input (see NM). 

If NREAD = 0, the initial streamline pattern estimate 
is generated by the program. 

If NREAD = 1, the initial streamline pattern estimate and 
also the DELF values are read in. (See DELF, R, X, 

XL.) 

Set NPUNCH = 0 
Set NPLOT « 0 

The maximum number of lines printed per page. 

An input value of zero will be reset to 60. 

IF NTRANS = 0, no action is taken. 

If NTRANS = 1, relative total pressure loss 
coefficients will be modified to account 
for radial transfer of wakes. See Section 
V.ll, Ref. 
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NMIX 


If NMIX ■ 0, no action is taken 


NMANY 


NSTPLT 


NEQN 


If NMIX * 1, entropy, angular momentum, 
and total enthalpy distributions will 
be modified to account for turbulent 
mixing. See Section V.12, Ref. 

The number of computing stations for 
which blade descriptive data is being 
generated by the analytic meanline 
section. 

If NSTPLT * 0, no action is taken. 

If NSTPLT * 1, a line-printer plot of the 
changes made to the midstreamline 
coordinate is made for each computing 
station. If more than 59 passes through 
the iterative procedure have been made, then 
the plots will show the changes for the 
last 59 passes. The graph should decay approx 
mately exponentially towards zero, indicating 
that the streamline locations are stabilizing. 
Decaying oscillations are equally acceptable, 
but, growing oscillations show the need for 
heavier damping in the streamline relocation 
calculations, that is, a decrease in RCONST. 
This item controls the selection of the 
form of momentum equation that will be used 
to compute the meridional velocity distri- 
butions at each computing station. There are 
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two basic forms, and for each case, one 
may select not to compute the terms 
relating to blade forces. (See also 
Section V. 1, Ref. 

If NEQN = 0, the momentum equation involves 
the differential form of the continuity 
equations and hence (1-M^ ) terms in the 
denominator. Streamwise gradients of 
entropy and angular momentum (blade forces) 
are computed within blades and at the blade 
edges (provided data that describe the 
blades are given) . Elsewhere, streamwise 
entropy gradients only are included in a 
simpler form of the momentum equation, 
except that at the first and last computing 
station, all streamwise gradients are taken 
to be zero. This is generally the preferred 
option when computing stations are located 
within the blade rows. 

If NEQN = 1, the momentum equation form is 
similar to that used when NEQN * 0, but 
angular momentum gradients (blade force 
terms) are nowhere computed. This generally 
is the preferred option when computing 
stations are Located at the blade edges only. 
If NEQN = 2, the momentum equation includes 
an explicit dVm/dm term instead of the (1-M X ) 
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denominator terms. All streamwise 
gradients (including blade force terms) 
are computed as for the case NEQN * 0. 

When computing stations are located within 
the blade rows, the results will generally 
be similar to those obtained with NEQN = 0, 
and solutions may be found that cannot be 
computed with NEQN - 0 due to high meridional 
Mach numbers. 

If NEQN * 3, the momentum equation is similar 
to that used when NEQN = 1 , but (as for the 
case NEQN = 1) no angular momentum gradients 
are computed. This may be used when computinq 
stations are located only at the blade edges 
and high meridional Mach numbers preclude the 
use of NEQN = 1 . 

NLE I 

NTE I See the Analytic Section. 

NSIGN I 

HWKICH The numbers of each of the computing stations 

for which blade descriptive data is being 
generated by the analytic meanline section. 

SCLFAC Linear dimension scale factor , see page 

An input value of zero will be reset to 

12 . 0 . 

TOLNCE Basic tolerance in iterative calculation 

scheme. An input value of zero will be 
reset to 0.001. (See discussion of 
tolerance scheme in Section VI, Ref. . ) 
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VISK 


SHAPE 


X SCALE 

PSCALE 

RLOW 

PLOW 

XMMAX 


RCONST 


CONTR 

CONMX 


Kinematic viscosity of gas (for annulus 
wall boundary layer calculations). An input 
value of zero will be reset to 0.00018. 
Units: LL/SCLFAC/SCLFAC/T. 

Shape factor for annulus wall boundary 
layer calculations. An input value of 
zero will be reset to 0.7. 


Set each equal to 0.0. 

The square of the Mach number that appears 
in the equation for the streamline relocation 
relaxation factor is limited to be not greater 
than XMMAX. Thus, at computing stations where 
the appropriate Mach number is high enough 
for the limit to be imposed, a decrease in 
XMMAX corresponds to an increase in damping. 

If a value of zero is input, it is reset 
to 0.6. 

The constant in the equation for the streamline relocation 
relaxation factor. The value of 8.0 that the analysis yields 
is often too high for stability. If zero is input, it is reset 
to 6. 0. 

The constant in the blade wake radial transfer calculations. 

The eddy viscosity for the turbulent mixing calculations. 
Units: L 2 /SCLFAC 2 /T. 
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FLOW 


Compressor flow rate. Units: F/T. 


SPDFAC 


NSPEC 


XSTN, RSTN 


NDATA 


NTERP 


N DIMEN 


NMA CH 


DATAC 


The speed of rotation of each computing station is SPDFAC 
times SPEED (I). The units for the product are revolutions/ 
(60xT). 

The number of points used to define a computing station. 
Must satisfy 2 s NSPEC s 21, and also the a uxn of NSPEC 
for all stations £ 150. If 2 points are used, the station ts 
a straight line. Otherwise, a spline-curve is fitted 
through the given points. 

The axial and radial coordinates, respectively, of a point 
defining a computing station. The first point must be on 
the hub and the last point must be on the casing. Units: L. 

Number of points defining conditions or blauc geometry at 
a computing station. Must satisfy 0 s NDATA s. 2 1 , and 
also the sum of NDATA for all stations s 100. 

If NTERP - 0, and NDATA i 3, interpolation of the data 
at the station is by spline-fit. 

If NTERP - 1 (or NDATA <, 2), interpolation is linear 
point-to-point. 

If NDIMEN - 0, the data are input as a function of radius. 

If NDIMEN = 1, the data arc input as a function of radius 
normalized with respect to tip radius. 

If NDIMEN = 2, the data are input as a function of distance 
along the computing station frepm the hub. 

If NDIMEN = 3, the data are input as a function of 
distance along the computing station normalized with 
respect to the total computing station length. 

If NMA CH = 0, the subsonic solution to the continuity 
equation is sought. 

If NMA CH = 1, the supersonic solution to the continuity 
equation is sought. This should only be used at stations 
vhcre the relative flow angle is specified, that is, 

NWORK - 5, 6, or 7. 

The coordinate on the computing station, defined according 
to NDIMEN, where the following data items apply. Must 
increase n.onotonically. For dim ensional cases, units are 
L. 
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DATAl 


DATA 2 


At Station 1 and if NWORK - 1, DATAl is 
total pressure. Units: F/L/L. 

If NWORK * 0 and the station is at a blade 
leading edge, by setting NDATA f 0, the blaae 
leading edge may be described. Then DATAl ia 
the blade angle measured in the cylindrical 
plane. Generally negative for a rotor, 
positive for a stator. (Define the blade 
lean angle (DATA3) also) . Units: A. 

If NWORK = 2, DATAl is total enthalpy. 

Units: H/F. 


If NWORK = 3, DATAl is angular momentum (radius times 
absolute whirl velocity). Units: LL/SCLFAC/T. 

If NWORK = 4, DATAl is absolute whirl velocity. Units: 
L/SCLFAC/T. 

If NWORK = 5, DATAl is blade angle measured in the 
stream surface plane. Generally negative for a rotor, 
positive for a stator. If zero deviation is input, it becomes 
the relative flow angle. Units: A. 

If NWORK = 6, DATAl is the blade angle measured in the 
cylindrical plane. Generally negative for a rotor, positive 
for a stator. If zero deviation is input, it becomes, after 
correction for stream surface orientation and station lean 
angle, the relative flow angle. Units: A. 

If NWORK = 7, DATAl is the reference relative outlet 
flow angle measured in the stream surface plane. Generally 
negative for a rotor, positive for a stator. Units: A. 

At Station 1, DATA2 is total temperature. Units: D. 

If N LOSS = 1, DATA 2 is the relative total pressure loss 
coefficient. The relative total pressure loss is measured 
from the station that is NL1 stations removed from the 
current station, NL1 being negative to indicate an 
upstream station. The relative dynamic head is determined 
NL2 stations removed from the current station, positive 
for a downstream station, negative for an upstream station. 
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DATA 3 


DATA4 
DATA 5 


DATA6 


DATA? 
DATA 8 


If NLOSS * 2, DATA2 i* the iscntropic efficiency of 
compression relative to conditions NLl stations removed, 
NL1 being negative to indicate an upstream station. 

If NLOSS = 3, DATA2 is the entropy rise relative to the 
value NL1 stations removed, NL1 being negative to 
indicate an upstream station. Units: H/F/D. 

If NLOSS = 4, DATA2 is not used, but a relative total 
pressure loss coefficient is determined front the trailing 
edge value and curve set number NCURVE of the NSET2 
families of curves. NLl and NL2 apply as for NLOSS - 1. 

If NWORK = 7, DATA 2 is the reference (minimum) 
relative total pressure los9 coefficient. NLl and NL2 
apply as for NLOSS = 1. 

The blade lean aV.gle measured troni the projection of a 
radial line in the plane of the computing station, positive 
when the innermost portion of the mde precedes the 
outermost in the direction ol rotor rotation. Units: A. 

The fraction of the periphery that is blocked by the presence 
of the blades. 

Cascade solidity. When a number of stations are used to 
describe the flow through a blade, values are only required 
at the trailing edge. (They are used in the loss coefficient 
re-estimatio:. procedure, and to evaluate diffusion tactoi s 
for the output. ) 

If NWORK = 5 or 6, DATA6 is the deviation _angle 
measured in the streamsurface plane. Generally negative 
(or a rotor, positive for a stater. Units: A. 

If NWORK = 7, DATA6 is reference relative inlet angle, 
to which the minimum loss coefficient (DATA 2) and the 
reference relative outlet angle (DATA7) correspond. 
Measured in the streamsurface plane and generally 
negative for a rotor, positive for a stator. Unit:,: A. 

If NWORK - 7, DATA 7 is the rate of change of relative 
outlet angle with relative inlet angle. 

If NWORK = 7, DA TAB is the relative inlet angle larger 
than the reference value at which the loss coefficient attains 
twice its reference value. Measured in the stream surface 
plane. Units: A. 
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DATA 9 


NWORK 


NLOSS 


If NWORK = 7, DATA9 is the relative inlet angle smaller 
than the reference value at which the loss coefficient attains 
twice its reference value. Measured in the strcamsurface 
plane. Units; A. 

If NWORK = 0, constant entropy, angular momentum, and 
total enthalpy exist along streamlines from the previous 
station. (If NMIX = 1, the distributions will be modified.) 

If NWORK = 1, the total pressure distribution at the com- 
muting station is specified. Use for rotors only. 

If NWORK = 2, the total enthalpy distribution at the com- 
puting station is specified. Use for rotors only. 

If NWORK = 3, the absolute angular momentum distribution 
at the computing station is specified. 

If NWORK = 4, the absolute whirl velocity distribution at 
the computing station is specified. 

If NWORK = 5, the relative flow angle distribution at the 
station is speciiied by giving blade angles and deviation 
angles, both measured in the stream surface plane. 

If NWORK = 6, the relative flow angle distribution at the 
station is specified by giving the blade angles measured 
in the cylindrical plane, and the deviation angles measured 
in the streamsurface plane. 

If NWORK = 7, the relative flow angle and relative total 
pressure loss cooffic ; cnt distributions are specified by 
means of an off-design analysis procedure. "Reference", 
"stalling", and "choking" relative inlet angles are 
specified. The minimum loss coefficient varies para- 
bolically with the relative inlet angle so that it is twice 
the minimum value at the "stalling" or "choking" values. 

A maximum value of 0. 5 is imposed. "Reference" 
relative outlet angles and the rate of change of outlet 
angle with inlet angle are specified, and the relative 
outlet angle varies linearly from the reference value 
with the relative inlet angle. NLOSS should be set to zero. 

If NLOSS = 1, the relative total pressure loss coefficient 
distribution is specified. 

If NLOSS = 2, the iscntropic efficiency (for compression) 
distribution is specified. 

If NLOSS - 3, the entropy rise distribution is specified. 
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NL1 


NL2 


NEVAL 


NCURVE 


If NLOSS = 4, the total pressure loss coefficient distribution 
is specified by use of curve-set NCURVE of the NSET2 
families of curves giving the fraction of final (trailing 
edge) loss coefficient. 

The station from which the loss(in whatever form NLOSS 
specifics) is measured, is NL1 stations removed from the 
station being evaluated. NL1 is negative to indicate an 
upstream station. 

When a relative total pressure loss coefficient is used to 
specify losses, the i dative dynamic head is taken NL2 
stations removed from the station being evaluated. NL2 
may be positive, zero, or negative; a positive vai c 
indicates a downstream station, a negative value indicates 
an upstream station. 

If NEVAL = 0, no action is taken. 

If NEVAL > 0, curve-set number NEVAL of the NSET1 
families of curve giving diffusion loss parameter as a 
function of diffusion factor will be used to re-estimate 
the relative tota* pressure loss coefficient. NLOSS must 
be 1, and NL1 and NL2 must specify the leading edge of 
the blade. See also NOEL. 

If NEVAL 0, curve-set number NEVAL is used as 
NAVAL 0, except that the re-estimation is only 
made after the overall computation is completed 
(with the input losses) . The resulting loss 
coefficients are displayed but not incorporated 
into the overall calculation. See also NDEL . 

When NLOSS = 4 , curve-set NCURVE of the NSET2 
families of curves, specifying the fraction of 
trailing-edge ;pss coefficient as a function 
of meridional chord is used. 
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NLITER 


NDEL 


N0UT1 

NOUT2 

NOUT3 


When NEVAL >0, up to NLITER re-estimations 
of the loss coefficient will be made at a 
given station during any one pass through 
the overall iterative procedure. Less than 
NLITER re-estimations will be made if the 
velocity profile is unchanged by re-estimat- 
ing the loss coefficients. (See discussion 
of tolerance scheme in Section VI, Ref .) 
When NEVAL * 0, set NDEL to 0. When 
NEVAL f 0, and NDEL > 0 , a component of the 
re-estimated loss coefficient is a shock 
loss. The relative inlet Mach number is 
expanded (or compressed) through a 
Prandtl-Meyer expansion on the suction 
surface, and NDEL is the number of points 
at which the Prandtl-Meyer angle is given. 

If NDEL = 0, the shock loss is set at zero. 
Must satisfy 05 NDEL121, and also the sum 
of NDEL for all stations < 100. 

Set NOUTl = 0 
Set NOUT2 = 0 

Thir data item controls the generation of 
NASTRAN - compatible temperature and 
pressure difference output for use in 
subsequent blade stress analyses. For 
details of the triangular mesh that is 
used, see the Output Description in 
Section 
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NBLADE 


N0UT3 « XY , where 

If X ■ 1, the station is at a blade 
leading edge. 

If X * 2, the station is at a blade 
trailing edge. 

If Y = 0, then both temperature and 
pressure data will be generated. 

If Y * 1, then only pressure data will 
be generated. 

If Y = 2, then only temperature data 
will be generated. 

If N0UT3 = 0, the station may be between 
blade rows, or within a blade row for 
which output is required, depending upon 
the use of NOUT3 / 0 elsewhere. See 
also description of NBLADE below . 

This item is used in determining the 
pressure difference across the blade. The 
number of blades is I NBLADE I . If NBLADE 
is positive, "three-point averaging" is 
used to determine the pressure difference 
across each blade element. If NBLADE is 
negative, "four point averaging" is used. 

(See the Output Description in Section .) 

If NBLADE is input as zero, a value of +10 
is used. At a leading edge, the value for 
the following station is used: elsewhere the 
value at a station applies to the interval 


- 34 - 


SPEED 


DELC 


DELTA 


V/ BLOCK 


upstream of the station. Thus by varying 
the sign of N BLADE, the averaging method 
used for the pressure forces may be varied 
for different axial segments of a blade 
row. 

This card is omitted if NDATA = 0. The 
speed of rotation of the blade. At a blade 
leading edge, it should be set to zero. 

The product SPDFAC times SPEED has units of 
revolutions/ (T x 60). 

The coordinate at which Prandtl-Meyer 
expansion angles are given. It defines 
the angle as a function of the dimensions 
of the leading edge station, in the manner 
specified by NDIMEN for the current, that 
is trailing edge station. Must increase 
monotonical ly . For dimensional cases, 
units are L. 

The Prandtl-Meyer expansion angles. A positive value 
implies expansion. If blade angles are given at the leading 
edge, the incidence angles are added to the value specified 
by DELTA. Units; A. (Blade angles are measured in 
the cylindrical plane.) 

A blockage factor that is incorporated into the continuity 
equation to account for annulus wall boundary layers. It 
is expressed as the fraction of total area at the computing 
station that is blocked. If NBL = 1, values (except at 
Station l) are revised during computation, involving data 
items VISK and SHAPE. 
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BBLOCK, 

BDIST 


NDIFF 

DIFF 

FDHUB 

FDMID 

FDTIP 

NM 

NRAD 

TFRAD 


DM 


WFRAC 


A blockage factor is incorporated into the continuity 
equation that may be used to account for blade wakes or 
other effects. It varies linearly with distance along the 
computing station. BBLOCK is the value at mid-station 
(expressed as the fraction of the periphery blocked), and 
BDIST is the ratio of the value on the hub to the mid- 
value. 

When NSET1>0, there are NDIFF points defining loss 
diffusion parameter as a function of diffusion factor. 

Must satisfy 1 £ NDIFF s 15. 

The diffusion factor at which loss parameters are specified. 
Must increase monotonicaliy. 

Diffusion loss parameter at 10 per cent of the radial blade 
height. 

Diffusion loss parameter at 50 per cent of the radial blade 
height. 

Diffusion loss parameter at 90 per cent of the radial blade 
height. 

When NSET2> 0, there are NM points defining the fraction 
of trailing edge loss coefficient as a function of meridional 
chord. Must satisfy 1 £ NM £ 11. 

The number of radial locations where NM loss fraction/ 
chord points are given. Must satisfy 1 £ NRAD £ 5. 

The fraction of radial blade height at the 
trailing edge where the following loss fraction/ 
chord curve applies. If NRAD =1, it has no 
significance. 

The location on the meridional chord where 
the loss fraction is given. Expressed as a 
fraction of meridional chord from the leading 
edge. Must increase monotonicaliy. 

Fraction of trailing edge loss coefficient 
that occurs at location DM. 
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The fraction of the total flow that is to 
occur between the hub and each streamline. 
The hub and casing are included, so that 
the first value must be 0.0, and the last 
(NSTRM) value must be 1.0. 

Estimated streamline radius. (These data 
are input from hub to tip for the first 
station, from hub to tip for the second 
station, and so on.) Units? L. 

Estimated axial coordinate at intersection 
of streamline with computing station. 
Units: L. 

Estimated distance along computing station 
from hub to intersection of streamline 
with computing station. Units: L. 

Station and streamline number. These 
are merely read in and printed out to 
give a check on the order of the cards. 



1.15.3.2 AERODYNAMIC OUTPUT DATA 


1. ANALYTIC MEANLINE SECTION 

Printed output may be considered to consist of four sections; a print- 
out of the input data, details of the blade sections on each stream surface , a 
listing of quantities required for aerodynamic analysis, and details of the 
manufacturing sections determined on the constant-z planes. These are 
briefly described below. In the explanation which follows, parenthetical 
statements are understood to refer to the particular case of the_jdouble- 
circular-arc blade (ISECN = 2). 

The input data printout includes all quantities read in, and is self- 
explanatory. 

Details of the strcamsurface blade sections are printed if IPRINT - 
0 or 1. Listed first are the parameters defining the blade section. These 
are interpolated at the stream surface from the tables read in. Then follow 
details of the blade section in "normalized" form. The blade section geometry 
is given for the section specified, except that the meridional projection of 
the chord is unity. For this section of the output, the coordinate origin is 
the blade leading edge. The following quantities are given: blade chord; 

stagger angle; camber angle; section area; location of the centruid of the 
section; second moments of area of the section about the centroid; orienta- 
tion of the principal axes; and the principal second moments of area of the 
section about the centroid. Then are listed the coordinates of the camber 
line, the camber line angle, the section thickness, and the coordinates of 
the blade surfaces. NPOINT values are given. 

A lineprinter plot of the normalized section follows. The scales for 
the plot are arranged so that the section jvist fills the page, so that the 
scales will generally differ from one plot to another. "Dimensional" details 
of the blade section are given next. The normalized data given previously is 
scaled to give a blade section as defined by IFCORD and CORD. For this 
section of the output, the coordinates . re with respect to the blade stacking 
axis. The following quantities are giv< •n: blade chord; radius and location 
of center of leading (and trailing) edge(s); section area, the second moments 
of area of the section about the centroid and the principal second moments 
of area of the section about the centroid. The coordinates of NPOINT points 
on the blade surfaces are then listed, followed by the coordinates of 31 
points distributed at (roughly) six degree intervals around the leading (and 
trailing) edges. Finally, the coordinates of the blade surfaces and points 
around the leading (and trailing) ed^e(s) is (are) shown in Cartesian form. 
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The quantities required for aerodynamic analysis are printed at ail 
computing stations specified by the IFANGS parameter. The radius, blade 
section angle, blade loan angle, blade blockage, and relative angular 
location of the camber line are printed at each stream surface intersection 
with the particular computing station. The blade section angle is measured 
in the cylindrical plane, and the blade lean angle is measured in the constant- 
axial-coordinate plane. 

Details of the manufacturing sections are printed if IPRINT = 0 or 2. 

At each value of z specified by ZINNER, ZOUTER, and NZ, section 
properties and coordinates are given. The origin for the coordinates is the 
blade stacking axis. The following quantities are given: section area; the 

location of the centroid of the section; the second moments of area of the 
section about the centroid; the principal second moments of area of the 
section about the centroid; the orientation of the principal axes; and the 
section torsional constant. Then the coordinates of NPOINT points on the 
blade section surfaces are listed, followed by 3 1 points around the leading 
(and trailing) edge(s). 

If NAERO = 1, the additional input and output required for, and 
generated by, the interface are also printed. (Apart from the input data 
printout, this is the only printed output when IPRINT = 3.) 

If the NASTRAN parameter PGEOM / -1 then cards are punched 
that may be used as input for the NASTRAN stress analysis program. 
For the purpose of stress analysis, the blade is divided into a 
number of triangular elements, each defined by three grid points. 
The intersections between computing stations and streamsur faces 
are used as the grid points and the grid points and element number 
ing scheme adopted is illustrated in Figure 1. 

TheNASTRAN input data format includes cards identified by 
the codes GRID, CTR1A2 and PTRIA2. The data are fully described 
in Reference 7, but briefly, the GRID cards each define a grid 
point number and give the coordinates at the grid point, the 
CTRIA2 cards each define an element in terms of the three 
appropriate grid points (by number, and in a significant order), 
the PTRIA2 cards each give an average blade thickness for an 
element . 
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CASinc, 



Figure i. NASTRAN Grid Point and Element Numbering 
Scheme, 

2. AERODYNAMIC SECTION 

a. Regular Printed Output 

The input data are first printed out in its entirety, and the results 
for each running point follow. The output is generally self-explanatory and 
definitions are given here for some derived quantities. Tabular output is 
generally not started on a page unless it can be completed on the same page, 
according to the maximum number of lines permitted by the input variable 
NPAGE. 


The results of each running point are given under a heading 
giving the running point number. Any diagnostics generated during the 
calculation will appear first under the heading. (Diagnostics are described 
in the following section. ) Then, a station-by- station print out follows for 






each station through to the last station, or to the station where the calcu- 
lation failed, if this occurred. One or more diagnostics will indicate the 
reason for the failure, in this event. Included in the mcshpoint coordinate 
data is the distance along the computing station from the hub to the inter- 
ception of the streamline with the station (L), and the station lean angle 
(GAMA). Where the radius of curvature of a streamline is shown as zero, 
the streamline has no curvature. The whi^l angle is defined by 


*t*an o( 


V* 

V* 


(1) 


For stations within a blade, or at a blade trailing edge, a rel- 
ative total pressure loss coefficient is shown. The loss of relative total 
pressure is computed from the station defined by the input variable NL1. If 
a loss coefficient was used in the input for the station (NLOSS = 1 or 4, or 
NWORK = 7), the input variable NL2 defines the station where the normalizing 
relative dynamic head is taken; otherwise, it is taken at the station defined 
by NL1. If the cascade solidity is given as anything but zero, it is used in 
the determination of diffusion factors. The following definition is used: 

^ Lr Ve, r V d2r 

■ v~ + TTvC 


Inlet conditions (subscript 1) are taken from the station 
defined bv the input variable NL1. 

The last term in Equation 2 is multiplied by - 1 if the blade 
speed is greater than zero, or the blade speed is zero and the preceding 
rotating blade row has negative rotation. This is necessary because 
relative whirl angles are (generally) negative for rotor blades and for 
stator blades that follow a rotor having "negative" wheel speed. Incidence 
and deviation angles are treated in the same way, so that positive and 
negative values have their conventional significance for all blades. 

If annulus wall boundary layer computations were made (NBL 
= 1), details are shown for each station. Then, an overall result is given, 
including a statement of the number of passes that have been performed and 
whether the calculation is converged, unconverged, or failed. When the 
calculation is unconverged, the number of me sh^points where the meridional 
velocity component has not remained constant to within the specified 
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tolerance (TOLNCE) on the last two passes is shown as IVFAIL. 
Similar ily, the number of streamtubes, defined by the huo and 
each streamline in turn, where the fraction of the flow is not 
within the same tolerance of the target value is shown as 
IFFAIL. If these numbers are small, say less than 10% of the 
maximum possible values, the results may generally be used. 
Otherwise, the computation should be rerun, either .for a 
greater number of passes, or with modified relaxation factor 
constants. The default option relaxation constants will 
generally be satisfactory but may need modification for some 
cases. If insufficient damping is specified by the constants, 
the streamlines generated will tend to oscillate and this 
may be detected by observing a relatively small ladius of 
curvature for the mid-passage streamline that also changes 
sign from one station to the next. This may be corrected by 
rerunning the problem (from scratch) with a lower value input 
for RCONST, say, of 4.0 instead of 6.0. When the damping is 
excessive, the velocities will tend to remain constant while 
the streamlines will not adjust rapidly to the correct locations. 
This will be indicated by a small IVFAIL and a relatively large 
IFFAIL. For optimum program performance, RCONST should be 
increased, and the streamline pattern generated thus far could 
be used as a starting point. The second constant XMMAX (the 
maximum value of the square of Mach number used in the relaxation 
factor) is incorporated so that in high subsonic or supersonic 
cases the damping does not decrease unacceptably. The default 
value of 0.6 may be too low for rapid program convergence in 
some such cases. 

If the generation of blade pressure load data for the 
NASTRAN program is specified (by the input variable NOUT3) , a 
self-explanatory printout is also made. The blade element 
numbering scheme is the same as that incorporated into both 
blading sections of the program, and illustrated in Figure 1. 
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If the loss coefficient re -estimation routine has been used for 
any bladerow(s) (NCVAL 4 0), a printout summarizing the computations 
made will follow. A heading indicating whether the re -estimation was 
incorporated into the overall iterative procedure or whether it was merely 
made "after the event" is first printed. Then follows a self-explanatory 
tabulation of various quantities involved in the redetermination of the loss 
coefficient on each streamline. 

b. Diagnostic Printed Output 

The various diagnostic messages that may be produced by the 
aerodynamic section of the program are all shown. Where a computed value 
will occur, "x" is shown here. 

JOB STOPPED - TOO MUCH INPUT DATA 

The above message will occur if the sun. of NSPEC or NDATA 
or NDEL for all stations is above the permitted limit. Execution ceases. 

STATIC ENTHALPY BELOW LIMIT AT xxx. xxxxxExxx 

The output routine (subroutine UD0311) calculates static 
enthalpy at each n.eshpoint when computing the various output parameters and 
this message will occur if a value below the limit (HMIN) occurs. The limiting 
value will be used, and the results printed become correspondingly arbitrary. 
HMIN is set in the Program UD03AR and should be maintained at some 
positive value well below any value that will be validly encountered in 
calculation. 


PASSxxx STATIONxxx STREAMLINE xxx PRANDTL-ME YER 
FUNCTION NOT CONVERGED - USE INLET MACH NO 

The loss coefficient re -estimation procedure involves iteratively 
solving for the Mach number in the Prandtl-Meyer function. If the calculation 
does not converge in 20 attempts, the above message is printed, and as 
indicated, the Mach number following the expansion (or compression) is 
assumed to equal the inlet value. (The routine only prints output following 
the completion of all computations and printing of the station-by -station 
output data. ) 

PASSxxx STATIONxxx ITERATIONxxx STREAMLINExxx 
MERIDIONAL VELOCITY UNCONVERGED VM = xx. xxxxxxExx 
VM(OLD) = xx. xxxxxxExx 
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For "analysis" cases, that is at stations where relative flow 
angle is specified, the calculation of meridional velocity proceeds 
iteratively at each meshpoint from the mid-streamline to the case and then 
to the hub. The variable LPMAX (set to 10 in Subroutines UD0308 and 
UD0326 ) 1 im its the maximum number of iterations that may be made at a 
streamline without the velocity being converged before the calculation 
proceeds to the next streamline. The above message will occur if ail iter- 
ations are used without achieving convergence, and the pass number is 
greater than NFORCE. Convergence is here defined as occurring when the 
velocity repeats to within TOLNCE/5.0, applied nondimensionally . No 
other program action occurs. 

PASSxxx STATIONxxx MOMENTUM AND/OR CONTINUITY 
UNCONVERGED W/W SPEC r xx. xxxxx VM/VM (OLD) HUB - 
xx. xxxxxMID-xx. xxxxx TIP = xx. xxxxx 

If, following completion of all ITMAX iterations permitted 
for the flow rate oi meridional velocity, the simultaneous solution of the 
momentum and continuity equations profile is unconverged, and the pass 
number is greater than NFORCE, the above message occurs. Here con- 
verged means that the flow rate equals the specified value, and the 
meridional velocity repeats, to within TOLNCE/5.0, applied nondimension- 
ally. If loss coefficient re -estimation is specified (NEVAL> 0), an 
additional iteration is involved, and the tolerance is halved. No further 
program action occurs. 

PASSxxx STATIONxxx VM PROFILE NOT CONVERGED WITH 
LOSS RECALC VM NEW/VM PREV HUB = xx. xxxxxx MID = 
xx, xxxxxx CASE = xx. xxxxxx 

When loss re - estimation is specified (NEVAL> 0), up to 
NLITER solutions to the momentum and continuity equations are completed, 
each with a revised loss coefficient variation. If, when the pass number is 
greater than NFORCE, the velocity profile is not converged after the 
NLITER cycles of calculation have been performed, the above message is 
issued. For convergence, the meridional velocities must repeat to within 
TOLNCE/5. 0, applied nondimensionally. No further program action occurs, 

A further check on the convergence of this procedure is to 
compare the loss coefficients used on the final pass of calculation, and thus 
shown in the station-by- station results, with those shown in the output from 
the loss coefficient re -e stimation routine, *vhich are computed from the 
final velocities, etc. 
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PASSxxx STATlONxxx ITERATIONxxx STREAMTUBExxx STATIC 
ENTHALPY BELOW LIMIT IN MOMENTUM EQUATION AT 
xxx. xxxxxExxx 

The static enthalpy is calculated (to find the static temperature) 
during computation of the "design" case momentum equation, that is, when 
whirl velocity is specified. If a value lower than HMIN (see discuss : on of 
second diagnostic message) is produced, the limiting value is inserted. If 
this occurs when IPASS > NFORCE, the above message is printed. If this 
occurs on the final iteration, the calculation is deemed to have failed, 
calculation ceases, and results are printed -u< rough to this station. 

PASSxxx STATlONxxx ITERATIONxxx STREAMTUBExxx LOOPxrtx 
STATIC H IN MOMENTUM EQUN. BELOW LIMIT AT xxx. xxxxxSxxx 

This corresponds to the previous message, but for the 
"analysis" case. For failure, it must occur on the final iteration and loop. 

PASSxxx STATlONxxx ITERATIONxxx STREAMTUBExxx 
MERIDIONAL MACH NUMBER ABOVE LIMIT AT xxx. xxxxxExx 

When Subroutine UD0308 is selected (NEQN = 0 or 1), the 
meridional Mach number is calculated during computation of the design 
momentum equation, and a maximum value of 0. 99 is permitted. If a 
higher value is calculated, the limiting value is inserted. If this occurs 
when IPASS > NFORCE, the above message is printed. If this occurs on 
the final iteration, the calculation is deemed to have failed, calculation 
ceases, and results are printed through to this station. 

PASSxxx STATlONxxx ITERATIONxxx STREAMTUBExxx LOOPxxx 
MERIDIONAL MACH NUMBER ABOVE LIMIT AT xxx. xxxxxExxx 

This corresponds to the previous message, but for the "analysis" 
case. For failure, it must occur at the final iteration and loop. 

PASSxxx STATlONxxx ITERATIONxxx STREAMTUBExxx 
MOMENTUM EQUATION EXPONENT ABOVE LIMIT AT xxx. xxxxxExxx 

An exponentiation is performed during the computation of the 
design case momentum equation, and the maximum value of the exponent is 
limited to 88. 0. If this substitution is required when IPASS > NFORCE, 
the above message is printed. If it occurs on the final iteration, the calcu- 
lation is deemed to have failed, calculation ceases, and results are printed 
through to this station. 
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PASSxxx STATIONxxx ITERATIONSxxx STREAMLINExxx 
(MERIDIONAL VELOCITY) SQUARED DELOW LIMIT AT 
xxx. xxxxxExxx. 

If a meridional velocity, squared, of less than 1. 0 is calcu- 
lated during computation of the design-case momentum equation, this limit 
is imposed. If this occurs when IPASS>NFORCE, the above message is 
printed. If this occurs on the final iteration, the calculation is deemed to 
have failed, calculation ceases, and results are printed out through to this 
station. 

PASSxxx STATIONxxx ITERATlONxxx STREAMLINExxx LOOPxxx 
(MERIDIONAL VELOCITY) SQUARED BELOW LIMIT AT 
xxx. xxxxxExxx. 

This corresponds to the previous message, but for the 
"analysis'* case. For failure, it must occur on the last iteration and loop. 

PASSxxx STATIONxxx ITERATlONxxx STREAMTUBExxx 
STATIC ENTHALPY BELOW LIMIT IN CONTINUITY EQUATION 
AT xxx, xxxxxExxx. 

The static enthalpy is calculated during computation of the 
continuity equation. If a value lower than HMIN (sec discussion of second 
diagnostic message) is produced, the limiting value is Imposed. If this 
occurs when IPASS>N FORCE, the rbove message is printed. If this 
occurs on the final iteration, the calculation is deemed to have failed, 
calculation ceases, and results are printed out through to this station. 

PASSxxx STATIONxx ITERATlONxxx STREAMLINExxx 
MERIDIONAL VELOCITY BELOW LIMIT IN CONTINUITY AT 
xxx. xxxxxExxx. 

If a meridional velocity of less than 1. 0 is calculated when the 
velocity profile is incremented by the amount estimated to be required to 
satisfy continuity, this limit is imposed. If this occurs when IPASS ^ 
NFORCE, the above message is printed. If this occurs on the final iteration, 
the calculation is deemed to have failed, calculation ceases, and results are 
printed through to this station. 

PASSxxx STATIONxxx ITERATlONxxx OTHER CONTINUITY 
EQUATION BRANCH REQUIRED 

If when IPA SS>N FO RCE, a velocity profile is produced that 
corresponds to a subsonic solution to the continuity equation when a super- 
sonic solution is required, or vice versa, the above message is printed. If 
thi s occurs on the final iteration, failure is deemed to have occurred, calcu- 
lation ceases, and results are printed out through to this station. 
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PASSxxx STATIONxxx ITERATlONxxx STREAMLINExxx 
MERIDIONAL VELOCITY GREATER THAN TWICE MID VALUE 


During integration of the "design" momentum equations, 
no meridional velocity is permitted to be greater than twice the value on 
the mid-streamline. If this occurs when IPA3S>NFORCE. the above 
message is printed. If this occurs on the final iteration, the calculation 
is deemed to have failed, calculation ceases, and results arc printed 
through to this station. In the event that this limit interferes with a valid 
velocity profile, the constants that appear on cards $08$. 272, $08$. 279, 

$26$ 229, and $26$. 236 may be modified accordingly. Note that as the 
calculation is at this point working with the square of the meridional 
velocity, the constant for a limit of 2. 0 tin.es the mid-streamline value, 
for instance, appears as 4.0. 

PASSxxx STATIONxxx ITERATlONxxx STREAMLINExxx 
LOOPxxx MERIDIONAL VELOCITY A130VE LIMIT xxxxxExx 
LIMIT = xxxxxExx. 

During integration of the "analysis" momentum equations, no 
meridional velocity is permitted to be greater than three times the value 
on the mid- streamline. If this occurs when IPASS>NFO RCE , the above 
message is printed. If this occurs on the final loop of the final iteration, 
the calculation is deemed to have failed, calculation ceases, and results 
are printed through to this station. In the event that the limit interferes 
with a valid velocity profile, the constants that appear on cards 
$035,398, $08$. 409, $26$. 323, $26$. 334, and $26$. 329 may be modified 
accordingly. In each case except that of the last card noted, the program 
is working with meridional velocity squared, so that a lin.it of, for instance, 
3.0 times the mid-streamline value appears as 9.0. 

PASSxxx STATIONxxx STREAMLINExxx LIMITING MERIDIONAL 
VELOCITY SQUARED = xxxxxExx. 

In the Subroutine UD0308 (NEQN= 0 or 1), a maximum 
permissable meridional velocity (equal to the speed of sound) is established 
for each streamline at the beginning of each pass. The calculation yields 
the square of the velocity, and if a value of less than 1. 0 is obtained, a 
value of 6250000.0 is superimposed (which corresponds to a meridional 
velocity of 2500. 0). If this’ occurs when IPASS >NFORCE, the above message 
is printed, and the calculation is deemed to have failed. Calculation ceases 
after the station computations are made, and results are printed through 
to this station. 


PASSxxx STATIONxxx ITERATIONxxx STREAMLINExxx 
MERIDIONAL VELOCITY ABOVE SOUND SPEED VM = 
xxxx. xx A = xxxx. xx. 

In Subroutine UD0308 (NEON - 0 or 1), no meridional velocity 
is permitted to be larger than the speed of sound. The above message will 
occur if this limit is violated during integration of the "design" momentum 
when IPASS > NFORCE. If the limit is violated at any point when IPASS> 
NPOF.CE and on the last permitted iteration (last permitted loop also in 
the case of the "analysis" momentum equation), the calculation is deemed 
to have failed. Calculation ceases, and the results arc printed through to 
this station. 


MIXING CALCULATION FAILURE NO. n 

The above message occurs when flow mixing calculations are 
specified, and the computation fails. The overall calculation is halted, and 
results are printed through to the station that is the upstream boundary for 
the mixing interval in which the fa 'ure occurred. The integer n takes on 
different values to indicate the specific problems as follows. 

n = 1 In solving for the static pressure distribution at the upstream 

boundary of each mixing step, the average static enthalpy is 
determined in each streamtube (defined by an adjacent pair 
of streamlines). This failure indicates that a value less than 
HMIN was determined. 

n = 2 Calculation of the static pressure distribution at the upstream 

boundary of the mixing step is iterative. This failure indicates 
that the procedure was not converged after 10 iterations. 

n = 3 The static enthalpy on each streamline at the mixing step 

upstream boundary is determined from the static pressure 
and entropy there. This failure indicates that a value less 
than HMIN was determined. 

n = 4 The axial velocity distribution at the mixing step upstream 

boundary is determined from the total enthalpy, static enthalpy, 
and tangential velocity distributions. This failure indicates 
that a value less than VMIN was determined. 

n = 5 In solving for the static pressure distribution at the downstream 

boundary of each mixing step, the average static enthalpy is 
determined in each streamtube (defined by an adjacent pair 
of streamlines). This failure indicates that a value less than 
HMIN was determined. 
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ii = 6 Calculation of the static pressure distribution at the down- 

stream boundary of the mixing step is iterative. This failure 
indicates that the procedure was not converged after 10 
iterations. 

n = 7 The static enthalpy distribution at the mixing step downstream 

boundary is found from the total enthalpy, axial velocity, and 
tangential velocity distributions. This failure indicates that 
a value less than HMIN was determined. 

n - 8 In order to satisfy continuity, the static pressure level at the 

mixing step downstream boundary is iteratively determined. 

This failure indicate? that after 15 attempts, the procedure 
was unconvereed. 

c. Aerodynamic Load and Temperature Output 

Four output options may result in cards being produced 
by the aerodynamic section of the program. Use of the input item 
NOUT3 gives "PLOAD2 and Temperature - Cards" punched in a format 
compatible with the NASTRAN stress program. For the purposes of 
stress analysis, the blade is taken to be composed of a number of 
triangular elements. Two such elements are formed by the quadrilateral 
defined by two adjacent streamlines and two adjacent computing stations. 
The way that each quadrilateral is divided into two triangles, and 
the element numbering scheme that is used, are illustrated in 
Figure 1, The pressure difference for eachelement is given by 
an average of either three or four values at surrounding mesh- 
points. The pressure difference at each meshpoint is computed 
from the equation 

A, ^ TT f p $ D Q T+‘ 1S +. V -2. i( rV “)l m 

il{> = fj ' Jt + r fa J 

and as follows. At the blade leading edge a forward difference is 
used to determine the meridional gradients. At the blade trailing 
edge the pressure difference is taken to be zero. At stations 
with the bladerow (following a leading edge) , mean central 
differences are used to determine the meridional gradients. When 
the input item NBLADE is positive (or zero) for a particular 


blade axial segment , then three-point averaging is used. For 
instance, for element number 1 in Figure 1, pressure differences 
at grid points 1, 6, and 7 would be used. If NBLADE is negative 
four-point averaging is used. For instance for element number 1 
pressure differences at grid points 1, 2 , 6 and 7 would be used. 
The same average would also apply to element number 2. Relative 
total temperatures are output at the grid points on the blade. 

A TEMPD value is also output using the average temperature at 
the blade root for the grid points on the rest of the structure. 

1.15.4 ^jMtsitJlc^Proble^ J 


USER’S MANUAL UPDATES 


1.15.4 Sample Problem 

The Static Aerothermoel astic Design/Analysis procedure 
for the bladed disc of an axial flow compressor rotor is 
illustrated by this sample problem. As explained in Section 
1.15,3 the Design and Analysis steps are carried out only at 
the design operating point of the compressor bladed disc - the 
"as manufactured 1 ' structure being only "analyzed" at off-design 
operating points. The Design ojr Analysis mode of the Displace- 
ment Rigid Format 16 is selected by the PARAMETER SIGN. The 
present example uses the Design mode (SIGN = -1) of the rigid 
forma t . 

The finite element model of a sector of the bladed 
disc is shown in Figure 1. The blade grid is specified in the 
Basic coordinate system located on the axis of rotation as shown 
in the figure. The hub is specified in a cylindrical coordinate 
system with the origin and the z-axis respectively coincident 
with the origin and the x-axis of the Basic system. A schematic 
of the aerodynamic model used is shown in Figure 2 wherein the 
aerodynamic mesh is generated by the intersection of 4 streamlines 
and 5 computing stations, three of which lie on the blade. Two 
additional computing stations have been used for the aerodynamic 
section (see Section 1.15.3.1), one each upstream and downstream 
of the blade to enable flow description in these regions. The 
NASTRAN deck for the use of the rigid format is listed in Figure 3. 
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Axis of Rotation 
(x-axIs) 


Figure 1. Finite Element Model of an Axial Flow Compressor 

Bladed Disc Sector, and the Basic Coordinate System 
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Tbe Executive Control Deck consists of cards from ID 
to CEND. S0L 16 and APP Displacement are used for the Steady 
Ar hermoelastlc Design/Analysis problem. CPU time (In minutes) 
is estimated on the TIME card. DIA6 (optional) Is used to 
request diagnostic output. 

The Case Control Deck Is used to select the boundary 
conditions imposed on, and the loads applied to the structure. 

'j extent and the form of the output desired is also selected 
in this deck. In this problem, SPC set 500 is used to restrain 
the hub-shaft attachment degrees of freedom from moving in the 
axial and tangential direction. MPC sot 600 Is used to define 
the blade-hub attachment and the relative motion of the corres- 
ponding grid points on the two sides of the cyclic sector. 

Two subcases must be defined for this rigid format. Subcase 1 
is for the linear solution based on the elastic stiffness while 
Subcase 2 solution includes the differential stiffness effects. 

The OUTPUT (PLOT) packet requests the plots, and is explained 
in Section 4. of the User's Manual. 

The blade is idealized by 12 CTRIA2 plate elements 
while 4 CHEXA1 solid elements are used to model the hub. The 
aerodynamic data describing the blade geometry (blade angle, 
chords, stagger angles etc.) and the operating conditions (flow 
rate, speed, losses etc.) are specified in the ALGDB data block 
input via the DTI bulk data cards. The geometry, material and 
constraint bulk data are as discussed in previous sections of this 
manual. Parameters APRESS = 1 and ATEMP = 1 enable the inclusion 
of the aerodynamic pressure and thermal loads. FXC00R, F Y C (9(9 R 
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and FZC00R parameters each equal to 0.3 Indicate that. In 
this design example, three tenths of the displacements obtained 
(both linear and non-linear) are used to redefine the blade 
geometry. Parameters IPRTCF = 1 and IPRTCI ■ 1 are used for 
a detailed printout from the ALG module upon final and initial 
entries. IPRTCL * 0 requests a summary from the ALG module 
during the differential stiffness loop (see Section 18 of the 
Theoretical Manual). PGE0M = 3 causes the GRID, CTRIA2, PTRIA2 
and DTI bulk data cards to be punched out during the final pass 
through the ALG module. These cards represent the final blade 
geometry and the operating conditions. Parameter STREAML = -1 
suppresses the output of STREAML1 and STREAML2 bulk data cards, 
while Z0RIGN = 0 only is currently permitted. STREAML1 cards 
identify the grid points defining the blade. 

Results are presented in the Demonstration Problems 

Manual. 
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1.15. 5 


Modal, Flutter 6 Subcritical Roots Analyses 

Cyclic symmetric flow is assumed while analyzing the 
turbomachinery rotor/stator. Due to rotational cyclic symmetry, 
only one-bladed disc sector is modeled. The harmonic number 
dependent cyclic normal modal analysis of such structures is 
described in Section 1.12 of the User's Manual. In the present 
development, the results of the normal modes analysis using 
cyclic symmetry have been appropriately integrated with unsteady 
cascade aerodynamic theories and the existing k-method of modal 
flutter analysis. The Mach number parameter has been conven- 
iently replaced by the interblade phase angle parameter for 
blade flutter problems. The discussion that follows is to 
bring out the features pertinent to bladed disc analysis. 

In a compressor or turbine, an operating point 
implies an equilibrium of flow properties such as density, velocity, 
Mach number, flow angle, etc., that vary across the blade span. 

Blade properties like the blade angles, stagger angle, chord,, etc., 
also, in general, change from the blade root to the tip. The 
resulting spanwise variation in the local reduced frequency and 
the relative Mach number must be accounted for in estimating 
the chordwise generalized aerodynamic forces per unit span at 
each streamline. Integration of these forces over the blade 
span yields the blade generalized aerodynamic force matrix. 

In order to nondimensionalize this matrix, the flow and blade 
properties at a referenced streamline are used. The reference 
streamline number, IREF, is specified cn a PARAM bulk data card. 

Since the relative Mach number varies along the blade 
span, necessitating the use of either the subsonic or supersonic 
cascade theories, parameters MAXMACH and MINMACH are used 
respectively to specify the upper and lower limits below and above 
which the subsonic and supersonic unsteady cascade theories are 
applicable. For streamlines with relative Mach numbers between 
the limits MAXMACH and MINMACH, linear interpolation is used. 

No transonic cascade theories have been incorporated. 


It should be noted that for a given interblade phase 
angle and reference reduced frequency, chordwise generalized 
aerodynamic matrices corresponding to local spacing, stagger 
and Mach number at the selected operating point will be generated 
for each streamline on the blade. This is an expensive operation 
and should be carefully controlled to reduce the computational 
work. The aerodynamic matrices are, therefore, computed at a 
few interblade phase angles and reduced frequencies, 
and interpolated for others. These parameters are selected on 
the MKAER01 and MKAER02 bulk data cards. Matrix interpolation 
is an automatic feature of Rigid Format Aero 9. Additional 
aerodynamic matrices may be generated and appended to the previous 
group on restart with new MKAER01 cards, provided the rest of the 
data used for the matrix calculation remain unaltered. 

To save further computational time, the chordwise 
generalized aerodynamic matrices are rirst computed for 
"aerodynamic modes" (see the Theoretical Manual, Section ). 

The aerodynamic matrices for chordwise structural modes are 
then determined from bilinear transformations along each 
streamline prior to the spanwise integration to obtain the 
complete blade generalized aerodynamic matrix. This permits 
a change in the structural mode shapes of the same or a 
different harmonic number to be included in the flutter analysis 
without having to recompute the modal aerodynamic matrices for 
aerodynamic modes. This can be achieved by appropriate ALTERS 
to the Rigid Format. 

For non-zero harmonic numbers, the normal modes 
analysis using cyclic symmetry results in both "sine" and 
"cosine" mode shapes (Section 1.12). The BCD value of the 
parameter MTYPE on a PARAM bulk data card selects the type of 
mode shapes to be used in flutter calculations. It is immaterial 
which is selected. 

The method of flutter analysis is specified on the 
FLUTTER bulk data card. The FLUTTER card is selected by an 
FMETH0D card. At the present time, only the k-method of flutter 
analysis is available. This allows looping through three sets 


of parameters: density ratio (P/p ref ,*p re ^ is given on AER0 
card); interblade phase angle (r); and reduced frequency, (k.) 
Por example, if the user specifies two values of each, there 
will be eight loops in the following order. 


L00P 

DENS 

<r 

RFREQ 

1 

1 

1 

1 

2 

2 

1 

1 

3 

1 

1 

2 

4 

2 

1 

2 

5 

1 

2 

1 

6 

2 

2 

1 

7 

1 

2 

2 

8 

2 

2 

2 


Values for the parameters are listed on FLFACT bulk data cards. 
Usually, one or two of the parameters will have only a single 
value . 

A parameter VREF may be used to scale the output 
velocity. This can be used to convert from consistent units 
(e.g., in/sec) to any units the user may desire (e.g., mph) , 
determined from V Qut = V /V REF . Another use of this parameter 
is to compute flutter index, by choosing V REF = bu>g T. 

If physical output (grid point deflections or element 
forces, plots, etc.) is desired rather than modal amplitudes, 
this data recovery can be made upon a user selected subset of 
the cases. The selection is based upon the velocity; the method 
is discussed in Section 3.23.3. 

V 1 ^ 1.576 &«dj|>le Problem j 
* > (To be wac4-. tten ) 


USER'S MANUAL UPDATES 


1.15.6 Sample Problem 

The problem of determining the complete, unstalled 
flutter boundaries of a compressor or turbine bladed disc 
Involves each member set of an appropriate whole series of 
harmonic families of modes of the cyclically symmetric bladed 
discs, and effects of Interblade phase angle, over an adequate 
set of operating points (flow rates, speeds, pressure ratios, 

Implied Mach numbers, etc.). This sample problem, therefore. 

Is only to Illustrate the procedure to obtain typical data 
leading to the definition of flutter boundaries. 

The finite element model of the compressor bladed 
disc sector is shown in Figure 1. The aerodynamic model (see | 

Section 1.15.2) with 4 streamlines and 3 computing stations is 
shown in Figure 2. The first four of the zeroth harmonic 
family of natural modes and frequencies are chosen for flutter 
investigation via the PARAMeters LM0DES * 4 and KINDEX = 0. 

Operating point conditions of 73.15 lb m/sec flow rate, 

16043 rpm, and 1.84 total pressure ratio are selected so as to 
demonstrate the use of the total stiffness matrix, for cyclic 
modal analysis, saved from the Static Aerothermoel asti c Analysis 
at this operating point (see Demonstration Manual examples 9-5-1 
and 16-1). For this, the Parameter KGGIN is set equal to 1. 

The k-method of flutter analysis is used which is the only 
method currently permitted. The NASTRAN deck used Is listed In 
Figure 3. 


The Executive Control Deck, cards ID through CEND, 
selects the Cyclic Modal Flutter Analysis Rigid Format via 
the S0L 9 and APP AER0 cards. An estimated CPU TIME of 
20 minutes Is Indicated for this example. The DIAG 14 card 
Is optional and lists the Rigid Format. 

The Case Control Deck Is used to select constraints, 
methods and output. In this problem, SPC set 500 Is used to 
constrain the hub-shaft attachment degrees of freedom to move 
only In the radial direction. MPC set 600 Is used to define 
the blade-hub connection. A METH0D card must select an EIGR 
bulk data card for real eigenvalue analysis. An FMETH0D card 
must be used to select a FLUTTER data card for flutter analysis 
k CMETH0D card must select an EIGC data card for complex 
eigenvalue extraction. For a flutter summary printout, the 
parameter PRINT Is set to YESB. The XYPAPERP10T request shown 
will plot V-g and V-f split frame "plots" on the printer output 
To produce plots, it is necessary to specify a plotter, request 
a plot tape, and specify XYPAPERPL0T VG. The "curves" refer 
to the loops of the flutter analysis, and in this example the 
9 loops have been arranged with 3 loops to each frame. 

The blade and the hub are respectively modeled by 
12 CTRIA2 and 4 CHEXA1 elements. The geometry, material and 
constraint bulk data are as discussed in previous sections of 
this manual, and there are no special rules for aeroelastic 
flutter analysis. CYJ0IN data card specifies the pairs of 
corresponding grid points on the two sides of the cyclic sector 
INV method of real eigenvalue extraction is selected on an EIGR 
card wherein five mode shapes and frequencies are requested. 


- 2 - 


0^ these, the first four (Parameter LM0DES - 4) modes are used 
to form the modal flutter equations. The AER0 bulk data card 
Is used to specify the reference chord and reference density. 

For bladed disc flutter analysis, the other two parameters on 
the AER0 card are of no significance. The MKAER01 data card 
causes the aerodynamic matrices to be computed for three Inter- 
blade phase angle-reduced frequency pairs, l.e. («-■ 180°, 
k = 0.3), (180°, 0.7) and (180°, 1.0). 

The FLUTTER bulk data card selects the presently 
permitted k-method of flutter analysis and refers to the FLFACT 
cards specifying density ratios, interblade phase angles, and 
reduced frequencies. The analysis loops through all combinations 
of densities, interblade phase angles and reduced frequencies, 
with density on the inner loop and Interblade phase angle on 
the outermost loop. In this example, 3 density ratios, 1 Inter- 
blade phase angle and 3 reduced frequencies (on FLFACT cards) 
result in (3 x 1 x 3 =) 9 loops. Both linear and surface splines 
are available for interpolation of aerodynamic matrices to 
intermediate values of interblade phase angle and reduced 
frequency. The E16C card is required and the HESS method is 
used. The number of complex eigenvectors to be extracted must 
be specified, and will usually agree with the number of modes 
saved for output specified on the FLUTTER data card. 

For bladed discs, STREAML1 and STREAML2 data cards are 
required. The grid points on each streamline on the blade are 
identified on the STREAML1 card. The flow and blade geometry 
is specified for each streamline on the STREAML2 cards. It 
should be noted that at least 3 streamlines per blade (including 
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the root and the tip) and 3 grid points per streamline must be 


selected for 
Res 


cyclic modal flutter analysis, 
ults are presented In the Demonstration Problems 
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Figure 3. HASTRAN DECK FOR CYCLIC MODAL FLUTTER ANALYSIS OF BLADED DISCS 
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NASTRAN DATA DECK 


T1HE K Required . 

K -- Maximum allowable execution time 


In minutes. 


S gt Kl [ ,K1 ] or S QL An f ,K1 1 Required when using a rigid format (see Section 3.1 
for available options). 

Kl -- Solution number of Rigid Format (see table below and Section 3.1). 

Kl -- Subset numbers for solution Kl , default value * 0. Multiple subsets 
may be selected by using multiple Integers separated by commas. 

An -- Name uf Rigid Format (see table below) 


Displacemen t Appt oac h Rigid Formats 
Kl An 


1 STATICS 

2 INERTIA RELIEF 

3 MODES or NORMAL MODES or REAL EIGENVALUES 

4 DIFFERENTIAL STIFFNESS 

5 BUCKLING 

6 PIECEWISE LINEAR 

7 DIRECT COMPLEX EIGENVALUES 

8 DIRECT FREQUENCY RESPONSE 

9 DIRECT TRANSIENT RESPONSE 

10 MODAL COMPLEX EIGENVALUES 

11 MODAL FREQUENCY RESPONSE 

12 MODAL TRANSIENT RESPONSE 

13 NORMAL MODES ANALYSIS WITH DIFFERENTIAL STIFFNESS 

14 STATICS CYCLIC SYMMETRY 

15 MODES CYCLIC SYMMETRY 

16 STATIC AEROTHERMOELASTIC ANALYSIS WITH DIFFERENTIAL STIFFNESS 


Heat Transfer App roach _R i d Formats 
Kl An 


1 STATICS 

3 STEADY STATE 

9 TRANSIENT 


Aeroelastic Approach Rigid Format 
Kl An 


9 COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 

10 MODAL FLUTTER ANALYSIS 

Subset Numbers 

1 . Delete loop control . 

2. Delete mode acceleration method of data recovery 
(modal transient and modal frequency response). 

3. Combine subsets 1 and 2. 

4 . ChecL all structural and aerodynamic data without execution of the 
aeroel a Stic probl em . 


2 . 2- 3a (9/30/78) 


I 

L 


L 

t 


k 


5, Check only the aerodynamic data without execution of the aeroelas- 
tlc problem. 

6, Delete checkpoint Instructions. 

7, Delete structure plotting and X-Y plotting. 

8, Delete Grid Point Height Generator. 

9, Delete fully stressed design (static analysis). 


I 


ij 


j 

1 


| 

i 


s 

7 


I ‘ 
? 


J 


2 . 2- 3b (9/30/78) 


CASE CONTROL DECK 


The following cards are o ed to make cufnnt requests for i,hn calculated response of com- 
ponents in the $3LU!IflN set (components in the direct or modal formulation of the general K 
system) for dynamics problems: 

1. SA C CEL THAT 1?\ - requests the acceleration of the Independent components for a selected 
set of points or nodal coordinates, 

2. SDJSPl ACt V r NT - requests the displa* em-nts of the independent components for a selected 
set of points or nodal coord mites or the temperatures of the independent components for 
a selected set of points in heat transfer. 

3. SVELDCIT Y - loqucsts the velocities of the independent components for a selected set of 
points or nodal coordinates or the change in tentperature with respect to time of the 
independent components for a selected set of points in heat transfer. 

4. NlLPAp - requests the nonlinear loads for a selected set of physical points (grid points 
and extra points introduced for dynamic analysis) intransient response oroblems. 

The following cards are used to make output requests for stresses and m^cs* as well as the 
calculated response of degrees of freedom used in the model: 

1. ELE gRCE - requests the forces in a set of structural elements or the temperature gradients 
and fluxes in a c et of structural or heat elements in heat transfer. 

2. STRESS - requests the stresses in a set of structural elements or the velocity components 
in a fluid element in acoustic cavity analysis. 

3. SPCFPRCES - requests the single-point forces of constraint at a set of points or the ther- 
mal power transmitted at a selected set of points in heat transfer. 

4. 0L0AD • selects a set of applied loads for output, 

5,. ACCEL E RAJ Jpj - requests the accelerations for a selected set of PHYSICAL points (grid, 

scalar and fluid points plus extra points introduced for dynamic analysis). 

6. OWLACEi! * requests the displacements for a selected set of PHYSICAL points or thc- 

tempera tures for a selected set of PHYSICAL points in heat transfer or the pressures 
for a selected set of PHYSICAL points in hydroelastic i ty . 

7. VELOCITY - requests the velocities for a selected set of f’HYSICAL points or the change in 
temperatures with respect to time for a selected set of PHYSICAL points in heat transfer. 

8. HARMONICS - controls the number of harmonics that will be output for requests associated 
with the conical shell, axfsyrivctnc solids and hydros last ic problems. 

9. CSE - requests structural element strain energies in Rigid Format l. 

10. G PFft RCF - requests grid point force balance due to element forces, forces of single point 
constraint, and applied loads in Rigid Format 1. 

11. THERMAL - requests temperatures for a set of PHYSICAL points in heat transfer. 

12. PRESSURE - requests pressures for a set of PHYSICAL points in hydroelastic i ty . 

13. CS£ - selects contact surface points to be output 
2.3.3 Subcas e Def i ni t ion 

in general, a separate subcase is defined for each loading condition. In statics problems 
separate subcases are also defined for each set of constraints. In complex eigenvalue analysis 


2 . 3-3 ( 9 / 30 / 78 ) 


CASE CONTROL DECK 


Case Control Data Card CSP - Contact Surface Point Selection 

Description : Selects the Interface contact surface points for a static 

aeroelastlc analysis. 

format and Examples : 

CSP • n 
CSP • 31 


Option : Meaning 

n Set Identification number of a CSP card (Integer > 0). 


Remarks : 

1. The normal displacement difference will be output for the selected 
Interface contact surface points. 

2, This card should select only those points of the Interface contact 
surfaces where "contact" constraint conditions were not invoked. Use 
the GPF0RCE Case Control Card to select points for which "contact" 
constraint conditions were invoked. 


2.3-1 la (9/30/78) 




Input Data Card CSJP 


SULK DATA DECK 


Contact Surface Points 


Description ; Defines Interface contact surface points for use In static 
aer ©elastic problems. 


Form at and Example : 
2 3 


♦ CSP1 


MeXg 
SID 

GA1 , GB 1 



CSP 

S 1 o 

GA1 

GB I 

GA2 

GB2 

GA3 

GB3 

CSP 

13 

5 

9 

10 

12 

13 

23 


♦ABC 

GA4 

GB4 

GA5 

GB5 

-etc- 


zn 



♦ ABC 


♦ CSP1 



Contents 

Identification number of contact surface set (Integer > 0) 

Grid point Identification numbers of node point pairs at 
interface contact locations (integer > 0). 


Remarks : 

1. Contact surface sets must be selected In the Case Control Deck (CSP * SID) 
to be used by NASI RAN 

2 . The normal displacement difference between each GA1 and GB1 pair will be 
output if this SID is selected. 

3. Only those points where “contact 11 constraints were not Invoked should be 
selected here. Contact surface points where “contact” constraints were 
Invoked should be selected by a GPF0RCL data card to output element forces 
at the contact locations. 


Interface contact surfaces represented by node pairs (GAl t GB1), (GA2, GB2), 
(GA3 , GB3 ) and ( GA4 , GB4) 


2 . 4- 66b (9/30'/8) 
























BULK DATA DECK 


Input Data Card FLFACT Aerodynamic Physical Data 

Description : Used to specify densities* Mach numbers or Interblade phase angles, 

and reduced frequencies for flutter analysis. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

e 

9 

10 

FLFACT 

SID 

FI 

F2 

F 3 

F4 

mm 

F 6 

OH 

ABC 

FLFACT 

97 

.3 

. 7 

3.5 







♦ BC ! r 

F 9 

-etc- 







4 










M ejjj contents 

SID Set Identification number (unique Integer > 0). 

FI Aerodynamic factor (real). 

Remarks : 

1. These factors must be selected by a FLUTTER data card to be used by NASTRAN. 

2 . Imbedded blank fields are forbidden. 

3. Parameters must be listed in the order in which they arc to be used within 
the looping of flutter analysis. 


2.4-1 16a (9/30/78) 

















BULK DATA DECK 


Input Oat# Card FLUT T ER Aerodynamic Mutter Oata 
Description : Oeflnes data needed to perform flutter analysis. 
format and Example : 


1 

2 

3 

4 

S 

6 

7 

8 

9 

10 

FLUTTER 

SID 

METHOD 

DENS 

MACH 

RFREQ 

IMETH 

NVALUE 



FLUTTER 

19 

K 

119 

219 

319 

S 

5 

. 




Field 

SID 

METHOD 

DENS 

MACH 

RFREQ 

IMETH 

NVALUE 


Set Identification number (unique Integer > 0). 

Flutter analysis method, “K M for k-method (BCD). 

Identification number of an FLFACT data card specifying 
density ratios to be used In flutter analysis (integer ^ 0) . 

Identification number of an FLFACT data card specifying 
MACH numbers or interblade phase angles (m) to be used in 
flutter analysis (Integer 0). 

Identification number of an FLFACT data card specifying 
reduced frequencies (k) to be used in flutter analysis 
( integer > 0) . 

Choice of interpolation method for matrix Interpolation 
(BCD: L - linear, S ■ surface, default is S). 

Number of eigenvalues for output and plots (integer > 0). 


Remarks : 

1. The FLUTTER data card must be selected in Case Control Deck (FMLTHCfD - SID). 


2 . The density is given by P • RHdREF where P Is ihe density ratio given on the 
FLFACT data card and RH0REF is the reference density given on the AER0 data 
card . 


3. The reduced frequency is given by k * (RE FC • w/2 * V ) , where REFC i, given on 
the AER0 data card, u is the circular frequency and V is the velocity. 


2.4*1 16c (9/30/78) 



BULK DATA OECK 


Input Data Card HKACRfM 


Mach Number - frequency Table 


Descr i p t i on : Provides a table of Macn numbers or interblade phase angles (m) 

and reduced frequencies (k) for aerodynamic matrix calculation. 


Format and Example : 


2 


MKAER01 


MKAER01 



♦ BC 

k i 

♦ BC 

. 3 



Field 

m. 


Contents 

List of Mr.ch numbers or interblade phase angles (Real, 

1 < 1 < 8 ). 

List of reduced irequencies (Real, 1 < j < 8). 


Remarks : 

1. Blank fields end the list, and thus cannot be used for 0.0. 

2. Combinations of (m,k) will be used. 

3. The continuation card is required. 

4. Marh numoers are input for wing flutter and interblade phase angles for 
blade flutter. 
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Input Data Card HKAERfl 2 


BULK DATA DICK 


Hach Number - frequency Table 


Description : Provides a list of Hach numbers or Interbladc phase angles (m) and 

reduced frequencies (k) for aerodynamic matrix calculation. 


Format and Examp le: 


NKAER0? 1-10 .30 .10 .60 .70 



Freld 
m ^ • k ^ 


Contents 

List of pairs of Mach numbers 0 ; Interblade phase angles (Real) 
and reduced frequencies (Real) (imbedded blank pairs are 
ski pped ) . 


Remarks : 

1. This card will cause the aerodynamic matrices to be computed for a set of 
parameter pairs. 

2 . Several MKAER02 cards may be in the deck. 

3. Mach numbers are input for wing flutter and interblade phase angle for blade 
flutter . 
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NASTRAN OATA DECK 


PARAM (Coni.) 


y. KMAX - optional In static analysis with cyclic symmetry (rigid format 
TTJ7 The Integer value of this parameter specifies the maximum value 
of the harmonic index. The default value is ALL which is NSEGS/2 for 
NSEGS even and (NSEGS-U/2 for NSEGS odd. 

z. KINDEX - retired in normal modes with cyclic symmetry (rigid format 
TSy. The integer value of this parameter specifies a single value of 
the harmonic index. 

a«. N0OJE - optional in modal flutter analysis. A positive integer of this 
parameter indicates user supplied downwash matrices due to extra points 
are to be read from tape via the 1NPUTT2 module in the rigid format. 

The defau 1 1 value Is - 1 . 

ab. PI , P2 and P 3 - required in modal flutter analysis when using N0OJC 
parameter . ~Tee Section 5.3.2 for tape operation parameters required by 
INPUT! 2 module. The defaults for PI, P2 and P3 are -1, 11 and TAPE10, 
respectively . 

ac. VKEf - optional in modal flutter analysis. Velocities are divided by 
The real value of this parameter to convert units or to compute flutter 
indices. The default value is 1.0. 

ad. PRINT - optional In modal flutter analysis. The BCD value, N0, of this 
parameter will suppress the automatic printing of the flutter summary 
for the k method. The flutter summary table will be printed if the BCD 
value Is YES for wing flutter, or YESB for blade flutter. The default 
is YES. 

ae. APRES S - optional in static ae ro thermoe ) as 1 1 c analysis. A positive 
TnTeger value will generate aerodynamic pressures. A negative value 
(the default) will suppress the generation of aerodynamic pressure loads. 

af . A TEMP - optional in static aerothermoelastic analysis. A positive 
Thteger value will generate aerooynamic temperature loads. A negative 
value (the default) will suppress the generation of aerodynamic thermal 
1 oads . 

ag. STREAML - optional in static aero thermoe! as t i c analysis. STREAML 8 ! 
causes the punching of STREAML1 bulk data cards. STREAML 8 2 causes the 
punching of STREAML2 bulk data cards. STREAML C 3 causes both STREAML 1 
and STREAML? cards to be punched. The default value, -1, suppresses 
punching of any cards. 

ah. * optional in static aerothermoelastic analysis. PGffSM * 1 
causes the punching of GRID bulk data cards. PGE0M * 2 causes the 
punching of GRID, CTRIA2 and PTRIA2 bulk data cards. PGE0M * 3 causes 
the punching of GRID cards and the modified ALGD8 table on DTI cards. 

The default, -1, suppresses punching of any cards. 

al. I PR T - optional In static aerothermoe 1 as t i c analysis. If I PR T > 0, 

t¥en intermediate print will be generated in the ALG module based on the 
print option in the ALGDB data table. If IPRT « 0 (the default), no 
intermediate print will be generated. 
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NASTRAN DATA OECK 


PARAH (Cont . ) 


ij. SIGN - optional In static aorothermoelastlc analysis. Controls the 
type of analysis being performed. SIGN * 1.0 for a standard analysis 
SIGN ■ -1.0 for a design analysis. The default is 1.0. 

ak. gjRlGN, FXCflgR, FYCggR, FZCjlg R - optional In rtetic *e ro th^rmoe 1 a s 1 1 c 
analys i s . fRcse are modi f feat ion factors. The defaults are 
20RIGN ■ 0.0, FXCUflR • 1.0, FYC00R ■ 1.0, and FZCMR *- 1.0. 

<1. H1NMACH - optional in blade flutter analysis. This is the minimum 
Mach number above which the supersonic unsteady cascade theory Is 
valid. The default is 1.01. 

am. HAXMACM - optional in blade flutter analysis. This is the maximum 
Mach number below which the subsonic unsteady cascade theory is valid. 
The default value is 0.80. 

an. I RE F - optional In blade flutter analysis. This defines the reference 
streamline number. IREF must be equal to a SIN on a STREAML2 bulk 
data card. The default value, -1, represents the streamline at the 
blade tip. If IREF does not correspond to a SLN, then the default 
will be taken. 

ao. HTYPE - optional In cyclic modal blade flutter analysis. This controls 
whTch components of the cyclic modes are to be used in the modal 
formulation. MTYPE • SINE for slnr components and MTYPE ■ COSINE for 
cosine components. The default BCj . lue is COSINE. 

ap. KTffUT - optional in static aero thermoe 1 a s 1 1 c analysis. A positive 

1 n teger of this parameter indicates that the user wants to save the 
total stiffness matrix on tape (GIN0 file INPT) via the 0UTPUT1 module 
In the rigid format. The default is -1. 

aq. KGG I N - optional in compressor blade cyclic mc^al flutter analysis. 
A“'pos i t 1 ve integer of this parameter indicates that the user supplied 
stiffness matrix is to be read from tape (GIN0 file INPT) via the 

I NPUTT 1 module in the rigid format. The default is -1. 
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BULK DATA OECK 


Input Data Card STREAML1 


Blade Streamline Data 


Description : Defines grid points on the blade streamline from blade leading 

edge to blade trailing edge. 


Format and Example: 


STREAML 1 

SLN 

G1 

STREAML 1 

3 

2 



G4 

G5 

8 

10 




68 69 -etc- 


Alternate Form: 


STREAML 1 

SLN 

GID1 

"THRU" 

GID2 

STREAML 1 

5 

6 

THRU 

12 



Field 

SLN 

6 1 , 6101 
Remarks : 


Contents 

Streamline number (integer > 0). 

6r1d point Identification numbers (Integer > 0) 


1. This card Is required for blade steady aeroelastic and blade flutter 
prob 1 ems . 

2. There must be one STREAML1 card for each streamline on the blade. 

For blade flutter problems, there must be an equal number of STRLAML1 
and STREAML2 cards. 


3. The streamline numbers, SLN, nust Increase with increasing radial 
distance of the blade section from the axis of rotation. The 
lowest and the highest SLN, respectively, will be assumed to 
represent the blade sections closest to and farthest from the axis 
of rotation. 

I. All grid points should be unique. 

5. All grid points referenced by GID1 through GID2 must exist. 

6. Each STREAML1 card must have the same number of grid points. The 

nodes must bt input from blade leading edge to the blade 

trailing edge ir the correct positional order. 
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BULK DATA DECK 


Input Data Card 5TREAHL2 Blade Streamline Data 


Description : Define aerodynamic data for a blade streamline. 


Format and Examp le: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

STREAML2 

$,N 

NSTNS 

STAGGER 

CHORD 

RADIUS 

BSPACE 

MACH 

DEN 

♦ abc 

STREAML2 

2 

3 

23.5 

1.85 

6.07 

.886 

.934 

.066 



♦ abc 

VEL 

FLOWA 








♦ ABC 

1014.2 

55.12 

- 








Field Contents 

-SIN Streamline number (Integer >0) 

NSTNS Number of computing stations on the blade streamline. 

(3 < NSTNS < 10, Integer) 

STAGGER Blade stagger angle (- 90 . 0 < s tagger <90.0, degrees) 

CHORD Blade chord (real >0.0) 

RADIUS Radius of streamline (real >0.0) 

BSPACE Blade spacing (real >0.0) 

MACH Relative flow mach number at blade leading edge 

(real >0.0) 

DEN fas density at blade leading edge (real >0.0) 

VEL Relative flow velocity at blade leading edge (real >0.0) 

FIOWA Relative flow angle at blade leading edge 

(-90.0 <FL0WA <90.0, degrees) 
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Remarks : 


1. At least three (3) end no more then ten (10) STREAM12 cards ere 
required for e blede flutter analysis. 

2. The streamline number, SIN, must be the seme e$ Its corresponding 
SLN on e STREAML1 cerd. There must be e STREAML1 cerd for eech 
STREAMLZ cerd. 

3. It is not required thet ell streemlines be used to define the 
eerodyneir c matrices used on blede flutter. 
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RIGID FORMATS 


Tht following rigid formats for structural analysis are currently Included 
In NASTRAN: 

1 . Static Analysis 

2 . Static Analysis with Inertia Relief 

3. Normal Mode Analysis 

4. Static Analysis with Differential Stiffness 

5. Buck 1 Ing Analysis 

6. Piecewise Linear Analysis 

7. Direct Complex Eigenvalue Analysis 

8. Direct Frequency and Random Response 

9. Direct Transient Response 

10. Modal Complex Eigenvalue Analysis 

11. Modal Frequency and Random Response 

12 . Modal Transient Response 

. 13. Normal Modes Analysis with Differential Stiffness 

14. Static Analysis with Cyclic Symmetry 

15. Normal Modes Analysis with Cyclic Symmetry 

16. Static Aerothermoelas t 1c Analysis with Differential Stiffness 

The following rigid formats for heat transfer analysis are deluded In 
NASTRAN: 

1. Linear Static Heat Transfer Analysis 
3. Nonlinear Static Heat Transfer Analysis 
9. ’transient Heat Transfer Analysis 

The following rigid formats for aeroelasttc analysis are Included in 
NASTRAN: 

9. Compressor Blade Cyclic Modal Flutter Analysis (Subsonic and Supersonic) 

10. Modal Flutter Analysis (Subsonic) 

3.1 .1 Input File Processor 

The Input File Processor operates in the Preface prior to the execution of 
the DMAP operations In the rigid format. A complete description of the operations 
In the Preface Is given In the Programmer's Manual. The main interest here Is to 
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Indicate the source of data blocks that are created In the Preface and hence appear 
only as Inputs In the OMAP sequences of the rigid formats. None of the data blocks 
created by the Input File Processor are checkpointed, as they are always regenerated 
on restart. The Input File Processor Is divided Into five parts. The first part 
(1FP1) processes the Case Control Deck, the second part (1FP) processes the 6uH Data 
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COMPRESSOR BLAOE MESH GENERATOR 


3.21 COMPRESSOR BLAOE MESH GENERATOR 
3.21.1 OMAP ScQucnce for Compressor Blade Mesh Generator 
RIGID FORMAT DMAP LISTING 
SERIES 0 

DMAP APPROACH. COMPRESSOR BLAOE MESH GENERATOR 

LEVEL 2.0 NASTRAN OMAP COMPILER - SOURCE LISTING 

OPTIONS IN EFFECT' GO ERR*2 NOLIST NOOECK NOREF NOOSCAR 

1 BEGIN $ 

2 ALG CASECC , , , , ALGDB , , / CASECCA .GE0M3A / C.N.-I / C.N.-I / 

V,Y,STREAMl*l / V,Y,PG£0M»2 / V.Y.IPRTH $ 

3 END $ 
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RIGID FORMATS 




1 

| 3.21.2 Description of DMAP Operations for Compressor Blade Moth Generator 

2. ALG generates GRID, CTRIA2 , PTRIA2 end STRlAMt ’ bulk data cards. These 
! cards are output via the system card punch. The GRID and CTRIA2 cards 

represent a compressor blade mesh. The aerodynamic input data is 
checked by performing an aerodynamic analysis. 
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COMPRESSOR BLADE MESH GENERATOR 


3.21.3 Output for the Compressor Blade Mesh Generato r 

The GRID, CTRIA2, PTRIA2 and STREAML1 bulk date cards are punched. 
Aerodynamic output is printed. 

3.21.4 Case Control Deck, DTI Table and Parameters for the Compressor Blade 
Mesh Generator 

1. Only TITLE, SUBTITLE and LABEL cards are processed, all other case 
control cards are Ignored. 

Z. The only required Input Is the ALGBO data table. This data block Is 
Input via Direct Table Input (DTI) bulk data cards. ALGDB contains 
all the aerodynamic Input necessary for the ALG module. For a 
detailed description of the ALGDB data block input see Section 1.15.3.1 
of the User 1 s Manua 1 . 

The following user parameters are used by the Compressor B1»de Mesh Generator. 

1. STRf AML - Optional - A value of 1 casues the punching of STREAMLl bulk 
data cards. A value of 2 causes the punching of STREAML2 bulk data 
cards. A value of 3 causes the punching of both STREAMLl and STREAML2 
cards. The default value, -1, suppresses the punching of all cards. 

2. PGEOH - Optional - A value of 1 causes the punching of GRID bulk data 

cards. A value of 2 causes the punching of GRID, CTRIA2 and PTRIA2 
bulk data cards. PGE0M » 3 causes the punching of GRID cards and the 
modified ALGDB table on DTI cards. The default value, -1, suppresses 
the punchir.g of all cards. 

3. I PRT - Optional - a non-negative value of this parameter will allow 
intermediate print to be generated by the ALG module based on the 
print option in the ALGDB data table. The default value, 0, 
suppresses all intermediate print. 
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STATIC AEROTHERMOELAST IC ANALYSIS MITH DIFFERENTIAL STIFFNESS 

3.22 Stattc Aerothermoelasttc Analysts with Differential Stiffness 

3.22.1 OHAP Sequence for Static Aerothermo a lasttc Analysis with 

Differential Stiffness. 

RIGID FJRMAT OMAP LISTING 
SERIES J 

OISPIACLMENI APPROACH, RIGID FORMA I 16 

LEVEL 2.0 NAS1RAN DMAP COMPILER - SOURCE LISTING 

DPT IONS IN EFFECT* GO -RR*2 NOLIST NOOECK NOREF NOOSCAR 


1 

BEGIN 

NO. 16 STATIC AERO IHEHMOE L ASTI C WITH DIFFERENTIAL STIFFNESS 

2 


GE31 l.GtUM2./GPL.E JC XI N.GPO T ,CS TM .BGPDT , SI 1/ V , N.LUS ET/ V, 
NOGPOT S 

3 

SAVE 

LUSET,NUGPDT A 

4 

ConcT^ 

erruri.nogpot s 

5 

CHKPNT 

GPL .EOEXIN.GPDT.C SIM.UGPOT.SIL » 

6 


GFOM2.E0EXIN/ECT S 

7 

ChKPNI 

ECT F 

8 

PAR AML 

PCOB//C »N,PRE S/C«N,/C,N,/C ,N,/V ,N,NOPCDB S 

9 

(gAW AMR) 

// C.N. COMPLEX / / V.V.SIGN / C.N.0.0 / V.N.CSIGN * 

10 

purge 

PLTSETX.PLTPAK.GPStT S,EL SE T S/N0PCD8 * 

11 

(£und^> 

P liNOPCOU S 

12 


PCOO.EJCXlNtlCT/PLISE TX ,PL TP AR ,GPSE T S , ELSETS/ V, N. NS IL/ V 

JUMP PLOT «- 1 S 

13 

SAVE 

NSIL i JUMPPLUT S 

14 

(pg hs£) 

PLISEIX// S 

15 

PAHAi 

//C,N,MPV/V,N ,PLIFL5/C,N.I /C,N,l S 

16 

PARA N 

//C,H,MPY/V,N,PF|lE/C,N,0/C,N,0 S 

17 

CtPNo3 

P 1, JUMPPLOT V 

18 

<£oT^ 

PLTPAR.GP SETS, ELSE IS .CASECC,' BGPDT ,E QEX I N ,S I L, , , ,/ P k UTXI/ 
NSIL /V, N.LUS? T/ V, N , JUMPPLOT/ V,N, PLTF LG/ V , N» Pr I L £ S 

19 

SAVE 

JUMPPLOT, PL TF LG ,PFILE A 

20 

(gjTnsj) 

PL0TX1// S 

21 

L ABEL 

P 1 $ 

*2 

CHKPNI 

PL I P AR , GP SE T S , F 1. SF T S $ 
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RIGIO FORMATS 


RIGIO FURKAT OMAR LISTING 
SIR I IS 0 


DISPLACEMENT APPROACH, RIGID FORMAT 16 

LEVEL 2.0 NASTR 4 N OMAP COMPILER * SOURCE LISTING 

t 


23 


26 

SAVE 

25 

PAKA* 

26 

Cmkpnt 

27 

(Sai 3 

28 

SAVE 

29 

<£ ONU^> 

30 

PURGE 

31 

chkpnt 

32 

PAMAH 

33 



36 

SAVE 

35 

CHKPNT 

36 

(CO N O 

37 


38 

CHKPNT 

39 

Lie EL 

60 

(gUNO ^ 

61 


62 

CHKPNT 

63 

LABEL 

66 


65 



GEDM 3t EtiEXINt GE0M2/SL l «CPT 1/ V «N,NOGRAV » 

NUGRAV I 

//C,N,ANO/V,N,NONGG/V,N,NOGRAV/V, V.GROPNT* -1 » 

SL T , GP T T $ 

ECT.EPI.OGPOT fSlL.GP TT.C STM/E ST, GEI .GPECT ,/ V,N, LUSET/ V,N, 
NOSIMP/C ,N , 1/ V,N , NOGFNL/ V, N,GENE L $ 

NOSIMP.NJGENL ,GCNEl $ 

ERROR l.NOSIMP t 

UCPST/GENEL » 

EST.GPEC T.GFI ,OGPST » 

/ A. , N »A00/ V *N tNOKGGX /C ,N,t/C ,N,0 $ 

E ST, C STM, MPT, 01 T.GEUM2 ,/KELM.KOICT ,MELM, MOICT , ,/V ,N,NOKGGX/ V, 
N ,NOMGG/C ,N ,/C ,N,/C,N,/C,T, COUP MASS/C, V , CPB AR/ C ,V ,C PMCO/C, V , 
CP3UAD1/C, »,;P0UAD2/C ,T,CPTR| At /C • V ,CPTK I A2/ C • V ,CPTUBE/C, V • 
CPOOPLT/C,»,:iMRPl l/C«Y«CPTROSC » 

NOK GG X ,NOMGG » 

K ELM , Kl) ICT , M; L M,MDIC T t 

JMPKGG.NOKGGX % 

GPECT »K0 1C! ,K ELM/RGG X,GP SI » 

K Gu X » CP ST $ 

JMPKCG $ 

JMPMGG f WUMGG % 

GPFCTrIOICT.MUM/PCC /C |V .WT HASS- 1.0 $ 

MGG % 

JHPHGG % 

9 

IBL UGROPNT $ 

ER«3R6,N0HGG % 
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STATIC AEROTHERMOELAST IC ANALYSIS WITH DIFFERENTIAL STIFFNESS 


aigio format omar listing 

SERIES 0 

DISPLACEMENT APPROACH, PI* 10 FORMA I lb 



I.CVEL l 

46 

c£pw<r~^ 

4? 


46 

LAO a 

49 

(EouiO 

SO 

CM*. Pi* I 

si 

<£ONU D 

5 1 

«m*T} 

51 

CHKPN 1 

54 

LABEL 

55 

PARAM 

56 


57 

SAVE 

56 


59 

PURGE 

60 

CWCPNT 

61 

C£uno3 

6 2 

JUMP 

63 

LAtiFL 

64 

CgONlT^) 

65 

<&£zE^> 

66 

SAve 

67 

CtpNO 3 

66 



ttGPO I ,C$TM,E9EXIN ,K»G /U* PMG / V t ‘ ,GROPNf/C, V.taTMASS » 

OGPMG. ,,,,// » 

LBl l » 

K GGX , KGG/NOGE NL » 

K GG % 

LBL I l.NOGCNL S 

GEI,KGGX/KGG/ V,N, LOSE !/V,N,NOGE NL/V.N.NOS IMP » 

KGG S 
L BL 1 < 

//L.N.MP V/V,N *NSX I P/C ,N ,0/C , N ,0 % 

CASECC,l»fcOMA,EUE Xl,», SI L ,*POT , BGPOT ,CST N/ RG, VS ,US ET , AS El /V,N, 
LUSE I/V,N,MP:F I /V ,N,MPCF//V,N,SI NGLE/ V ,N , OMIT/ V , N, REACT/ V,N • 

N SK |P/V,N,k6PEAT/V,N , T/V,N,NCL/V.N,NUA/C, V.SUBIO S 

MPCFl.MPCF 2, SINGLE ,OMI I ,REACT,NSKI P,REPE AT ,NOSF.I .NOL.NOA » 

ERROR SiNUL t 

GN/MPCF 1 /GO ,R00 »LOU« PO , UOOV »RUU V/OMI T/ PS ,KFS ,KSS , QG/S INGLE/ 
U0U3V/OMIT/VB S ,PB S , K6F S, <0 SS • KOF S.KOSS/S INGLE A 

GM,R G ,G0 ,KUO, LOO, PO, UOOV • ROOV « VS ,PS «KF S , KSS «USET • AS ET , 0 0 00V, 

VBS.PBS.KBF S, KBSS.KUF S.KOSS.OG S 

L BL 40, HE AC I » 

error; s 

L BL AO » 

LBL A, G-NElt 

GPL ,*PSt ,USCI , SIL/QGPST/V,N,NOGPSI S 
NOGPST % 

LBLA.NUGPSr » 

OGPST, » 
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RIGID FORMATS 


RICIO FOMA? DMAP LISTING 
SERIES 0 

DISPLACEMENT APPROACH, H IG 10 FORMAT 16 

LEVEL i . 0 NASTKAN ONAP COMPILER - SOURCE LISTING 


69 

LABEL 

LBLL A 

?0 

<fouiv3 

KGG,KNN/MPCF| $ 

7 I 

CHKPNT 

KNN > 

n 

CCONQ 3 

LN.2.MPCF2 A 

1 3 

<$ceT^> 

USE?, KG/CM A 

16 

CHKPN T 

CM 1 

15 

<gcTn> 

uslt,gm,kgg», ,/KNN,, , A 

16 

CHKPNT 

KNN A 

7? 

LAOEL 

LBL2 A 

76 

<guu iO 

KNN, KFE /SINGLE A 

79 

CHKPNT 

KFF > 

•0 

(JuniT^ 

LW. 3, SINGLE $ 

6! 


USE T, KNN,,, /KFE ,KFS,KSS,,, A 

62 

CHKPN? 

Kl S,KSS,KEF A 

63 

LAOEL 

LBLJ A 

64 

(gjautQ 

KFF.KAA/OMIT A 

85 

CHKPNT 

KAA I 

66 

<|UNO 3 

LttLS.OHII A 

87 

C|mpP> 

USE T ,KE F , , , /G 0 , KA A ,KQO ,L00 , , , , , A 

68 

CMLPN T 

G0,KAA,KOO,L3O A 

69 

LAOEL 

LttL 5 A 

90 

<g0MGT^» 

KAA/LLL A 

91 

CHKPNT 

LLl A 

92 

dscT^) 

SL T «UCPOT«C STM , SI L,E ST,MPT,GPTT ,EOT ,MGG, CASECC, OIT/ PGNA / V,N, 
L USE T/C ,N, l A 

93 

CHKPN T 

PGNA ; 
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i 


k ICIO FQRH AT OH AP I I STING 
f C$ 0 

OISPLACEHEN7 APPROACH, R IS I U f GRHA f 
LEVEL 2.0 NASTRAN 0H4P COMPILER 


16 

* SOURCE LISTING 


94 PAR AH 

95 c SjniT3 

96 < £T c 3 


91 

Coho 

90 

PAR AH 

99 

CONO 

100 

GP3 

101 

r«KPHT 

102 

qscT^) 


10) 

CHKPNT 

104 

(ACC__^) 

105 

LABEL 

106 

<fbu I 

107 

CHKPNT 

100 

(gouiv^) 

109 

CHKPNT 

110 

C^ONO~^> 

111 

CSscT^) 

112 

CHLPNT 

113 

LABEL 

114 

CfscT^ 

115 

SAVE 

116 

CHKPNT 

117 

<£ono 


//CtN.AHO /V, NrAIOAO /V,V,APRESS /V.V.ATEHP $ 

NOAL , AlOAO $ 

C ASECC • « EQE Xl Nt ,A LGOd » * / CA SECCAl ,GE OK) A1 /S , V , APRtSS/S »V, 
ATEiP/C,N,-l/C ,N,-1/V. V, I P«TCI/$,N,lFAIt I 

F IX I S# IEA IL t 

//C,N,A HO / Vi W #AL GAO /V ? V,APRESS /V, V, Af EHP I 
NOAL t ALOAO S 

GCUM 3 Al,EQEX| N.GE 0 H 2 /SL TAl ,GPT TAl / V,N.NOGRAV $ 

SlTAl.GPTTAl % 

Sl.TAl tOGPOr,: STH, SI t ,E S 1 ,HPV »GPT TAl • EOF » MGG, CAS ECCA1»0IT / 

PGA1 / V »N » LU St T / C »N»1 S 

PGA l % 

PGXA,PGAl / PG $ 

NUAl % 

P GNA • PG/ALOAO % 

P 0 % 

PG.Pl/’IOSfcT I 
PI % 

l 61 I Of NO SE T I 

USE T # G^rYS,KF SiGO, ,PG/»PO,PS,Pt » 

PU»PS»PL $ 

L 61 10 S 

t LI ,KAA tPLfCOOfKOC#PO/UlVfUOO^iRULV,RUOV/V,N,OHU/V,Y f l RES"-!/ 
CtN, l/V«NfEPSI % 

EPSI $ 

Ul V, UCUV.RULV ,KUUV t 
L BL 9 , IKES $ 


i 


RIGID FORMATS 


RIGID FORMAT JMAP LISTING 
SEA I ES 0 

DISPLACEMENT APPROACH, RIGID I CRH A I 16 

LEVEL 2.0 NASTRAN OMAP COMPILER - SOURCE LISTING 


tie 

(Jutgph) 

i 19 

(jjATGPj) 

1 20 

LABEL 

121 


i n 

chkpnt 

123 

3 

24 

PARAH 

25 


i 26 

SAVE 

27 

<£0N0 D 

20 

(gLOT 3 

29 

SAVE 

30 


31 

LABEL 

32 


33 

(gSMG~r> 

3 V 

CHKPNT 

35 

(gOMO 

36 

CTuui Q 

37 

LABEL 

38 

PAR AM 

39 

PAM AM 

40 

PARAMK 


GPL ,USE T «SIL,RULV//C ,N,L A 
GPL , USE T , SIL .RUOV//C ,N ,U A 
L Bt 9 > 

USET, .ULV.UOOV, VS, GO.GM.PS, KF S.KSS./UGV.PGl ,QG/C,N, l/C.N.OSO A 

UGV, UG A 

CASfcCCit STM, MPT, 01 T.EQE'.l N.SI L.GPTT ,EOT, 8GP0I , , QG.UGV ,EST , , PG/ 
OPGl.OuGI.OUJVl.UESl ,OEFl ,POG VI /C ,N,l)SO I 

//C,I,1PV/V,N,CAAONO/C,N,0/C ,N,0 t 

3UGV1,JPG1,U3G1,UEFI ,Otil , // V.N.CARDNO A 

C AH O N L A 

P2.JUMPPL0T A 

PLTPAR.GPSETS.ELSE I S ,CA*ECC .BGPDT ,EQEX I N ,S I L.PUGVl, ,GPECT,0ES1/ 
PLUTX2/V,N»NS1L/V,N,LLSE T/V «N tJUMPPLQT/V *N»PLTF LG/V , N, PF ILE A 

PFILE A 

PL0TX2// A 

P 2 » 

F.CT.EPT.BGPDT ,SIL ,GP T I ,C STM/ XI , X2 ,ECPT .GPCT/ V.N.LUS ET / V,N, 

.NOS IMP/C ,N.O/V.N,NOGENL/V,N,GENEL A 

CASECC, GPTT,SIL,EOT,UGV,CSTM,MPT ,ECPT ,GPCT ,011/ KOGG/ V.N, 
OSCOSETA 

KOGG A 

NUALOtALOAJ A 
P GNA ,PG A 
NUALO A 

//c,n,aoo/v,n,shifi/c.n,-i/c,n,o a 

//C,N,AOO/V,V,COUM/V,N,ALMAVS»1/V,N, NEVER* I A 
//C,A,40l>/V,N,USEPSI /c,n,o. o/c.n.o.o a 


3 . 22-6 ( 9 / 30 / 78 ) 


STATIC AEROTHERMOEIASTIC ANALYSIS KITH DIFFER! IAL STIFFNESS 


I 

I 

I 

I 


RIGIO F3 A A AT OMAP LISTING 
SERIES 0 

DISPLACEMENT APPROACH. RIGIO FORMAT 16 

LEVEL 2.0 NASTAAN OMAP COMPILER - SOURCE LISTING 


I 

I 


m 

142 

U3 

144 

145 

146 

147 
146 

149 

150 

151 

152 

153 

154 

155 

156 

157 
156 
1'9 
160 
161 
162 

163 

164 

165 


PAR AML 

CJUMP^ 



1 S//C.N* NULL/C tNt/CfN'/C »N#/VtN,NQV$ ft 
OUTtP TOP I 
OLITLPTOP * 

POtPGl/NUYS ft 
PGl ft 

/ /C • N • KL OC K « V » N » f 0 I 
K DOG I KONM/MP: F 2 ft 
K0V9 $ 

L 6L 20 1 MPCF 2 1 
USETiGMfKDGGi tt/KONNttt ft 
KONN ft 
L6L2U ft 

KDNN , KDFF /SINGLE ft 
KOFF ft 

L 6L 3D« SINGLE ft 

USE T « KONN iff/ KOFF §KOF $« KD SS » » » ft 
KDFFfKOFSfKObS ft 
L6L30 ft 

K OFF # KDAA /OMI T ft 
KCAA ft 

L BL 50 «0M IT ft 

USET t GO f KOFF/KOAA ft 

KOAA e 

L6L50 ft 

KAA v KOAA / K3 LL / C • N. (1 • 0.0. 0) / V*N.CS1GN ft 


3.27-7 (9/30/78) 




RIGID FORMATS 


RIGID fHA® LISTING 

SERIFS 0 

01 $PL Af f MFN T AP»*01CM f «JG!0 fCRRAT 16 

IEV ®5 ?.C .V.STRAN D*AP C 0 4 PILCR - SO'JPCE LISTING 


166 

167 

169 

169 

173 

in 

17 ? 
173 
7 A 
i 75 

76 

77 
79 
79 

,eo 

8 ? 

83 

e<. 

85 

66 

87 

89 

89 

90 




EQUlV 

COKC 


L AOH 



t ADEL 



SAVE 


CWK* M T 

<^TpTt> 

< ^P TP 

<gr*r^ 

O'e* Q 

PAf A * 

< £p'^ 3 
< 5 TT^ 


KFi.KOfS/ KOFS / C,K,U.3,0.0I / V.N.CSIGN * 

<«S«0<$/ KflSS / C>K»(1.0|0.0) / V.N.CSIC.N » 

PGIK.NDYS i 

KTS’l.TS, /»>SS .N,0/C,N,l/C,N,l/C,N,l I 
KBFA.YS, fOfi/r. ,*|,0/C ,N,l/C,N,l/C.N, 1 » 
lSf:T ,PF S.PSS/Pf./C .N,K/C ,N ,*/C ,K,S » 

PN,P*,X/MPCF2 « 

L9LE f ) ,*<PCP2 * 

US" T ,PN, /l»GY/ r ,N,0/C ,N,N/C ,N f * » 

IBIV^ % 

PG<,PG/P'»G/f. ,\,T -1.0.0. 0> t 
PG~.,PG1/At. WAYS $ 

PGGK » 

1,/PGJ/ l 
UGv / A UG V t 

<9Ll/lPLL/V,N,PnwfP/V t N,0ET $ 

PcT.PCt»"P i 
IBll t 

//G.fJ.C/CiK.OET » 

/ r. ,'4,0/C ,N,PrKET i 
ItMLP’CP \ 

IMPTCP % 

//r,,N ,KLrCK/V,K,T| » 

NOALl, Atcan » 

ca',rc.c.FTT,pQ e sn.atjGv.atGr.j.cSTM.nr.pnT / iaskca.cfomja /SiV, 

APR e SS/S,Y,ATf PP/C.N.-l/C ,N.-1/V,Y, tPPTCl/S.N, I f 1 1 1/ V . V , <H CN/ V, 


3 . 22-8 ( 9 / 30 / 78 ) 


STATIC AEROTHERMOELASTIC ANALYSIS WITH DIFFERENTIAL STIFFNESS 


RIC1C FORMAT OMAP LISTING 
SERIES 3 

01 SPl * r EMEN T APPROACH, RIGID FORMAT 16 

LEVU. 2*0 NAS TRAN UMAP COMPILER - SOURCE LISTING 


191 (£ UNd3 

192 PARAS 

193 PARAM 
19* CONO 

195 GP3 

196 ( fSGi 3 


15 7 

198 

199 

200 
201 
202 



L AbfcL 
EQUI V 
CHKPNT 



20 3 
20 A 

205 

206 
207 
206 

209 

210 
211 
212 
213 


SAVE 

CHKPNT 



chkpnt 



y.ior ign/v,y,fxco: . ,f ycoor/v.y.fzcoor % 

DONE , IF A It S 

//C,N,Mpy / V, V, IPRTC L /C,N,0 I 
/ /C , N ,ANO / V, N , ALOAD /V,Y,APRESS /V,Y,ATEMP $ 

NOALliALOAO $ 

GEOM 3A,£OLAlN,GEOM2/SLTAfGPTT A/ V ,N, NOASL/ V , N, NJGRAV/ V , N, NOATL $ 

SL TA, $ GP OtfCSTMtSILfEST, HP TfGPTTA, EOT, MG G, CAS EC CA,D!T / PGA /V, 

NtLUSEI /C , N , IS 
PG1,PGA / PG2 % 

NJAL l $ 

PGi,PG2 / ALOAO * 

PG2 » 

USLr,GH,YStK3FSfGU,,PG2 / , PdO ,P0S, PBL % 

t 61 L ,K0LLfPBL , , , /Ul V , ,R UttL V • /C , N , -l f V ,Y • I RES/ V » N, NOSK I P/Vt N, 
EPSI » 

EPSI % 

UdLVtRUdLV $ 

L0L9J,|RC$ $ 

GPL » L SE T i 51 Lt RUBL V//C • N, L > 

L 01 93 $ 

USE I, tuHLV, ,V StGO,GM,POS,KflFSfR8SSt/UBGVt,OBG/C|N,l/C,N,OSl S 
UOGV,OBG I 
NUAC 2 , ALOAD % 

UdGViAUGV % 

HUAI 2 * 

U07V, UGV/DUGV/C ,N , I- l . 0, 0. 01 $ 


3.22-9 (9/30/78) 


RIGID FORMATS 


RIGID FORMAT OMAP LISTING 
SIR I E$ 0 

DISPLACEMENT APPROACH. RIGID FORMAT 16 

LEVEL 2.0 NASTRAN OMAP COMPILER * SOURCE LISTING 


21* (gSMGT^) 


2 IS 
216 
217 
2 1H 


CFXPNT 



219 

220 
221 
222 
223 
22 * 

225 

226 
227 
226 

229 

230 

231 

232 

233 
23* 

235 

236 

237 


SAVE 


CMKPNT 



CASECCi GP T T # S I L .6 0 I tOUG V ,C ST M ,MPT .ECPf ,G PCI # 01 T / OKDGG/V , N, 

0 SCO SET 6 

DR DUG % 

0K0GG.U3CVfPG0/PGI l/CiNtO/CiN*! /CiNtl/CfNil 6 
PGll.PGA / PG 1 2 i 

PG2 f PGl2#UaGV //C tY,EPSI0»l.E-5 /V,N#DSEPSI / C.VfNT*IO /V#N* 
TO / V »N r T I /V #N »OONE /V,N, SHIFT /V t N .COUNT/ 1, Y » GET AO* * S 

OSEP S I »30NE tSHlFT «C0UNT $ 

DONE# DONE $ 

SHIFT. SHIFT I 

PGfP G 1/HE VER $ 

PGt UPG1/ALWL VS % 

PC|,PGIi/NE VER $ 

INLP TOP, 1000 $ 

PGI 1 # PG l f PG • • // $ 

SHI F T % 

QKOGGfKOGG/lOGGl/C ,N , < - 1 . 0 , 0 . Q) $ 

KOGGl $ 

UBGVt UGV/ALW* V S /KDGG 1 » KOGG /ALWAYS $ 

KUGG % 

K OGGi K3GG l/NE VCR / LG V , UdG V/NE VER % 

OUTL P TJP t i 000 * 

KDGG l *<OGG • IG V t » / / % 

OONt $ 

//CtNtNOP / V # Y |K TOUT % 

JMPK rOUT»KTOUT % 


3 . 22-10 ( 9 / 30 / 78 ) 


I 


STATIC AEROTHERMOELASTIC ANALYSIS WITH DIFFERENTIAL STIFFNESS 


RIGID FORMAT DM/P LIS ING 
SERIES 0 


I DISPLACEMENT APPROACH, H ID ID FORMAT 16 

' LEVEL 2.0 NASTR4N UMAP COMPILER - SOURCE LISTING 

I 


230 ( ACC 

239 CuTPuT y) 

240 ( Output )- 

241 LABEL 

242 CHKPNT 
24 3 (SuT 


244 ( ?UR2 3 

245 ( QfP 3 

246 SAVE 

247 C jDk l 3 



249 
2 50 
251 



252 SAVE 


253 ( PRTiS$ .) 

254 LABEL 
265 ( JUMP 
25b LABEL 
257 ( gkTPAjO fr 
25B LABEL 
259 (g RTPAR ^) 


KGO.KOGG / KTOTAL / C.N, 11.0,0.01 / V.N.CSIGN » 

KTDTAL.i.. // C .Y.LOCATION— 1 / C , V , l NPTUNI T»0 » 

.... // C.N.-3 / C.N.O » 

JMPK TOUT % 

C STM » 

CASECC.EDT ,E9E XIN.UttGV.ALGDB . CST R.BGPDT / C ASECCB.GEUNJB /C.N. 
-l/C.N.-l/V.V.STKEAML/V.V.PGECM/V.V.IPRTCF/S.N, IFAIL/V.V.SIGN/ 
V./.20HIGN/V. Y.FXCUUH/V, V.F VCOOR/V.Y.F2COOR > 

CASECC.C STM.MPT.31 T.E0EK1N.SI L.GPT T .EOT, BGPOT , .OBG.UOGV , EST . . /, 
OUBGl.OUBGVl.OE SB 1 .OEFBt .PUBGV1 /C.N.0S1 % 

OUUGV 1 ,0(JBG1, OEFB 1 ,0E SOI , , // V.N.CARDNO » 

CAR DNU » 

USET.PG2.UBLV..YS.G0 .GM.PB S.KBF St KBSS , / AUDGV , APGG , A (JOG /C.N. 

1 /L.N.DSt » 

CASECC . AUOGV, KE LM ,K0 1C T ,EC T ,EUE XI N.GPECT , APGG. AQBG /ONRGVl, 
0GPFU1 .* C.N. STATICS » 

ONRGVl.OGPFBl .... // % 

P3. JUMPPLOT t 

PLTPAK.GP SETS,ELSETS,CASECC,BGP0T,E0EXIN,SIL,PUDGV1,.GPECT, 

OESB 1/PL0TX3/ V.N, MSI L /V, N, LOSE T/ V .N . JUMP PLOT/ V, N. PCTFLG/V, N, 
PULE t 

PFILr t 

PL0TX3// $ 

P 3 t 

FHIS t 

ERROR 1 % 

//C.N,- 1/C.N.DFrF STIF » 

ERROR 2 $ 

//C.N,- 2/C, N.DIFF STIF 6 


l 


3.22-11 (9/30/78) 


RIGID FORMATS 


RIGID FORMAT OMAR LISTING 
SERIES 3 

DISPLACEMENT APPROACH. R IS 10 FORMAT 16 

LEVEL 2.0 NASTRAN OMAP COMPILER - SOURCE LISTING 


260 

261 

262 

263 

266 

265 


LABEL ERROR* $ 

<^RTPAK$) //f .N 6/C .N.O IFF S II F $ 

LABEL ERROR 5 I 

^RTPAR?J> //C.N t~ 5/C .N.O IFF S 71 F $ 

LABEL FINIS % 

END > 

• •NO ERRORS FOUNO - EXECUTE NASTRAN PROGRAMS 


3 . 22-12 ( 9 / 30 / 78 ) 


6P1 generates coordinate system transformation matrices, tables of grid 
point locations, and tables for relating internal and eiternal grid point 
numbers . 

Go to DMAP No. 256 If no grid point definition table. 

GP2 generates Element Connection Table with Internal indices. 

PARAMR sets CS!GN*(S!GN, 0.0), where SIGN is ♦l.O or -1.0 for analysis or 
design type run. 

Go to DMAP No. 21 if no plot package is present. 

PtTSET transforms user input into a form used to drive structure plotter. 

PRTMSG prints error messages associated with structure plotter. 

Go to OMAP No. 21 if no undeformed structure plot request. 

PL0T generates all requested undeformed structure plots. 

PRTMSG prints plotter data and engineering data for each undeformed plot 
generated. 

GP3 generates Static Loads Table and Grid Point Temperature Table. 

TA1 generates element tables for use in matrix assembly and stress 
recovery. 

Go to DMAP No. 256 and print error message if no structural elements. 

EMG generates structural element matrix tables and dictionaries for later 
assembly. 

Go to DMAP No. 39 if no stiffness matrix is to be assembled. 

EMA assembles stiffness matrix [K* ] and Grid Point Singularity Table. 

y y 

Go to OMAP No. 43 If no mass matrix Is to be assembled. 

EMA assembles mass matrix [M^]. 

Go to OMAP No. 48 If no weight and balance request. 

5. Go to OMAP No. 260 and print error message If no mass matrix exists. 

6. GPWG generates weight and balance Information. 

0FP formats weight and balance information and places It on the system 
output file for printing. 

Equivalence [K*g] to [K^] if no general elements. 

Go to OMAP No. 54 If no general elements. 

SMA3 adds general elements to [K*g] to obtain stiffness matrix [K^]. 

GP4 generates flags defining members of various displacement sets (USET)» 
forms multipoint constraint equations [R Q ]{u } ■ 0 and forms enforced 
displacement vector { Y } . 9 9 


3.22-13 (9/30/78) 


kIGIO formats 


58. 

61. 

6'4. 

65. 

67. 

68 . 

70. 

72. 

73. 
75. 


78. 


80. 


81. 


83. 

86. 

87. 


Go to DHAP No. 262 and print orror message If no Independent dogroct of 
freedom ore dtflnod. 

Go to DMAP No. 63 If no free-body supports supplied, otherwise go to OMAP 
No. 258. 

Go to DHAP No. 67 If general elements present. 

GPSP determines If possible grid point singularities remain. 

Go to DHAP No. 69 If no Grid Point Singularity Table. 

6FP formats table of possible grid point singularities and places It on 
the system output file for printing. 

Equivalence [K^] to [K„ n ] If no multipoint constraints. 

Go to DHAP No. 77 If HCE1 and HCE2 have already been executed for current 
set of multipoint constraints. 

MCE1 partitions multipoint constraint equations [R g ] • [R m ;R n ] and solves 
for multipoint constraint transformation matrix [G„] * -[R_]"'[R_]. 

MCE 2 partitions stiffness matrix 


[K 


99 


] « 


1 



and performs matrix reduction 

■ t«„„> * “lit"..) * 

Equivalence [K nf) ] to [K^] if no single-point constraints. 
Go to OMAP No. 83 If no single-point constraints. 

SCE1 partitions out single-point constraints. 


["n.I 


K ff * K fs 

6<T‘» 


Equivalence [ ] to [K a# ] if no omitted coordinates. 
Go to DMAP No. 89 if no omitted coordinates. 

SMP 1 partitions constrained stiffness matrix 


[Kff] 

solves for transformation matrix [G Q ] 
and performs matrix reduction [K ] 

Q Q 


1^0. I'.oJ 

C«..I * 


3.22-14 (9/30/78) 




STATIC AEROTHERMOELASTIC ANALYSIS WITH DIFFERENTIAL STIFFNESS 


90. RMBG2 decomposes constrained stiffness matrix [K |4 J ■ [L t( ][U tt ]. 

Ui 

92. SS61 generates non-aerodynamlc static load vectors {P£ ). 

95. Go to DHAP No. 105 If no aerodynamic loads. 

96. ALG generates aerodynamic load data. 

102. S5G1 generates aerodynamic load vector (P*). 

Hk k 

104. Add {P } and {P*} to form total load vector {P }. 

9 NA 9 

106. Equivalence (P^) to (P£ ) If no aerodynamic loads. 

108. Equivalence (P # ) to (P ( ) If no constraints applied . 

110. Go to DMAP No. 113 If no constraints applied. 

111. SSG2 applies constraints to static load vectors 


. (P f ) - {P f }-[K fs ]{Y s ) . 

and {P t > - (P 4 )*[gJ]{P 0 ) . 

114. SSG3 solves for displacements of Independent coordinates 

{u t> • C'a.J'^V • 

solves for displacements of omitted coordinates 

<»;> ■ tK 00 r . 

calculates residual vector ( RUL V ) and residual vector error ratio for 
Independent coordinates 

(«P t ) • (P t ) - [K aa ]{u t > 

(u!h«P,) 
c l . 

(pj)<u t > 

and calculates residual vector ( RU0V ) and residual vector error ratio 
for omitted coordinates 

' <P 0 > * tK 0 C ]{uJ) , 

{ Uq } ( 6 P 0 ) 

c ° ' < tu 2> 

0 0 


<v 


I 1 

(pj * {- 


(p f ) 


3.22-15 (9/30/7 ) 


RIGID FORMATS 


117. Go to OMAP No. 120 If residual vector* ere not to be printed. 

118. Print residual vector for Independent coordinates (RUIV). 

119. Print residual vector for omitted coordinates (RU#V). 

121. SDR1 recovers depend! at displacements 

<“ 0 > • t G 0 i {u t> ♦ <«;> . 



and recovers single-point forces of constraint 

<9 $ > « -1P S ) ♦ CKj s ]{u f ) ♦ [K ss ]{¥ s }. 

122. S0R2 calculates element forces and stresses (0£F1, 0ES1) and prepares load 
vectors, displacement vectors and single-point forces of constraint for 
output (RPG1 , 0UGV1 , PUG Vi , 0OG1 ) . 

125. PFP formats tables prepared by SDR2 and places them on the system output 
file for printing. 

127. Go to DMAP No. 131 If no static deformed structure plots are reguested. 

128. P10T generates all requested static deformed structure plots. 

130. PRTMSG prints plotter data and engineering data for each deformed plot 
generated. 

132. TA1 generates element tables for use In differential stiffness matrix 
assembly. 

133. 0SMG1 generates differential stiffness matrix [K^„]. 

99 

135. Go to DHAP No. 137 If no aerodynamic loads. 

136. Equivalence (P^) to {P^} to remove aerodynamic loads from total load 

vector before entering (inferential stiffness loop. New aerodynamic loads 
will be generated In iou*. 

142. Go to next DHAP Instruction If cold start or modified r ‘art. RUTIPTRP 
will be altered by tne Executive System to the proper 1 Ion Inside the 
loop for unmodified restarts within the loop. 

143. Beginning of outer loop for differential stiffness Iteration. 

144. Equivalence (P^) to {P^} If no enforced displacements. 

147. Equivalence t K ggl t0 tK„ n 3 if no multipoint constraints. 


3.22-16 (9/30/78) 


STATIC AEROTHERMOELAST 1C ANALYSIS WITH DIFFERENTIAL STIFFNESS 


149. 

150. 


1 S3. 

155. 

1 56. 


1 59. 
161. 
162. 


165. 

166. 

167. 

168. 

169. 

170. 

171. 
174. 

176. 

177. 


6o to OHAP No. 1S2 If no Multipoint constraints. 
NCE2 partitions differential stiffness Matrix 


[K 






and perforMS Matrix reduction uj n ) ■ [kJ„3 ♦ 

Equivalence [K d n ] to [K d ^3 IT no single-point constraints. 
Go to DNAP No. 158 If no single-point constraints. 

SCE1 partitions out single-point constraints 


[K d 1 

1 nn J 


k ff 

,d 


I K d . 


fs 

d 


llsf | K ssJ 


Equivalence [K d ^] to [K 44 3 If no omitted coordinates. 

Go to DHAP No. 164 If no omitted coordinates. 

SMP2 partitions constrained differential stiffness matrix 


t'Jf) 


I IC** 

*aa *ao 

~ + ~ 


oa 


oo 


and performs matrix reduction [K?_] ■ [K^ 1 ♦ [K*]*[G A ] 

do A A 0 A O 

• |C 0 ) T [< 0 )t6 0 J. 

AD° [K 4# ] and [K aa ].CSIGN to form [K^J. 

,r '0 [K fs ] and [K d s ].CSIGN to form [K^ $ ] . 

ADD [K $s ] and [K d s ].CSIGN to form C s 3 • 

Go to ONAP No. 178 If no enforced displacements. 

MPYAD multiply [Kj $ ] and { V 8 ) to form {P $$ ). 

NPYAD multiply [Kj $ 3 and (Yj) to form { P } . 

UHERGE expand { P f s > and { P $ s } to form { P n } . 

UNERGE expand ( P n } to form {P*}. 

ADD - { P* > and (P.) to form ( P } . 

9 9 99 

Equivalence (P^) t0 ^ P gl^‘ 


♦ 


Ca 0 ] T Cx 



3.22-17 (9/30/78) 



RIGID FORMATS 


ADD (Pg)) and nothing to create (Pg # ). 

Copy (u^) to (Ug) to Initialize aerodynamic displacements . 

RBHG2 decomposes the combined dlfferentlel *t1ffne»* metrlx end elestlc 
stiffness metrlx. 

PRTPARM prints the seeled value of the determlnent of the combined dlfferen- 
tlel stiffness metrlx end elastic stiffness metrlx. 

PRTPARM prints the scele factor (power of ten) of the determinant of the 
combined differentia) stiffness metrlx end the elastic stiffness matrix. 

Go to next DMAP instruction if cold start or modified restart. 1NLPTRP will 
be altered by the executive system to the proper location inside the loop 
for unmodified restarts within the Hop. 

Beginning of inner loop for differentia) stiffness iteration. 

Go to DMAP No. 194 if no aerodynamic loads. 

ALG generates aerodynamic load data. 

Go to DMAP No. 235 if ALG fails to converge while generating aerodynamic 
load data. 

SSG1 generates aerodynamic load vector (P*). 

ADD (P .} and {P*} to form total load vector (P -). 
gl g gZ 

SSG2 applies constraints to static load vectors 


(P 


92 ' 





■ <p$> • t*; s ]{v s }, 

■ <pJ> ♦ [gJhpJi . 


SSG3 solves for displacements of Independent coordir tes for current differ- 
ential stiffness load vector 


( uf) • [K^r’fpji 


i 

1 

i 

5 

s 

* 

l 

3 

1 

i 

1 

\ 

i 


i 


\ 


3.22-18 (9/30/78) 


STATIC AEROTMEAMOELASTIC ANALYSIS WITH DIFFERENTIAL STIFFNESS 


205. 

206. 
208. 


210 . 
211 . 

213. 

214. 
216. 


217. 

218. 
220 . 
221 . 

222. 


•Ad calculates residual vacte. (RBULV ) *nd residual vector orror ratio for 
current differential stiffness load vector 

(6PJ) - (Pj) - t »},](«{) . 
b . <^) t up; > 

4 <"p}) T {u$) 

Go to DMAP No. 207 If rasldual vector for current d 1 f f e»*en 1 1 « 1 stiffness 
solution Is not to be printed. 

Print residual vector for current differential stiffness solution. 

S0R1 recovers dependent displacements for current differential stiffness 
solution 



[GoKuj) ♦ ( U ; b ) 

. < u i> ■ 



• (U b > . 



and recovers single-point forces of constraint for current differential 
stiffness solution 


iqj) • *{p£MK$ f 3{u5WK b f ]{Y b ) . 
Go to DMAP No. 212 If no aerodynamic loads. 

Equivalence { u ® } to lu*) . 

9 9 

ADD -(U b ) and (U ) to form { > . 

9 9 9 

DSMG1 generates differential stiffness matrix [4Kgg] . 
MPYAD form load vector for Inner loop Iteration. 


(P 0 )- kJjiuJmp 1 

9 j , 99 9 90 


ADD (P„ ) and {P*} to form (P } . 

9,1 9 9 I2 

DSCHK performs differential stiffness convergence checks. 

Go to DMAP No. 235 If differential stiffness Iteration Is complete. 

Go to DMAP No. i27 If additional differential stiffness matrix changes are 
necessary for further Iteration. 

Equivalence breaks previous equivalence of (Pg) to {Pg,} . 
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223 . 

< 4 , 

tn. 


equivalence (P g ^) t0 ^gt^’ 

Equlvtlcnce breaks previous equivalence of {F.«) to (P„ J. 

•' 9 n 

Go to OMAP Mo >87 for additional Inner loop differential stiffness 
1 teratlcn. 


226. TABPT tabic prints vectors (P n ), { P _ . } » and (PI. 

9j) 9' 9 

228. A00 -[«J g ] and [kj g ] to form [nj gl 2. 

230. Equivalence (U g ) to (U g ) and [K^ g j] to [K gg ]. 

23* Equivalence breaks previous «* iu 1 va lence of [K*J 1 to [K^„] and {U„ } to (U b ). 

99 99 9 9 

233. Go to DMAP No. 143 for additional outer loop differential stiffness Iteration. 

234. TABPT table prints [K< gl] , tK J gJ ^ (|| ^ 

237. Go to DMAP No. 241 If the total stiffness matrix Is not to be saved on tape. 

238. ADD [K gg 3 and [Kjj ] to form [KTOTAl] . 

239. 0UTPUT1 outputs [KTOTAL] to tape. 

240. OUTPUT! prints the names of the data blocks on the output tape. 

243. AlG generates final aerodynamic results and generates GRID and STRCAML2 bulk 
data cards on the system punch, If requested. 

244. SDR2 calculates element forces and stresses (0EFB1, 0ESB1) and prepares 
displacement vectors and single-point forces of constraint for output 
(0UBGV1, PUBGV1 , 0OBG1) for all differential stiffness solutions. 

245. 0FP formats tables prepared by S0R2 and places them on the system output 

file for printing. 

247. SDR1 recovers dependent displacements after differential stiffness loop for 
grid point force balance. 

248. GPFDR calculates for requested sets the tirld point force balance and element 
strain energy for output. 

249. OFP formats the tables prepared by GPFDR and places them on the system output 
file for printing. 

250. Go to OMAP No. 254 If no deformed differential stiffness structure plots are 
requested. 

251. PL0T generates all requested deformed differential stiffness structure plots. 

253. PRTMSG prints plotter data and engineering data for each deformed plot 
generated. 

255. Go to DMAP No. 264 and make normal exit. 

257. STATIC ANALYSIS WITH DIFFERENTIAL STIFFNESS ERR0R MESSAGE N0. 1 - N0 
STRUCTURAL ELEMENTS HAVE BEEN DEFINED. 

259. STATIC ANALYSIS WITH DIFFERENTIAL STIFFNESS ERR0R MESSAGE N0. 2 - FREE B0DY- 
SUPP0RTS N0T ALLOWED. 
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261 . STATIC ANALYSIS WITH DIFFERENTIAL STIFFNESS ERR0R MESSAGE N0. 4 - NASS 
MATRIX REQUIRED FOR WEIGHT AND BALANCE CALCULATIONS. 

263 STATIC ANALYSIS WITH DIFFERENTIAL STIFFNESS ERROR MESSAGE NO. 5 - NO 
INDEPENDENT DEGREES OF FREEDOM HAVE BEEN DEFINED. 
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3.22.3 Automatic Output for Static Aerothermoelastlc Analysis with 
Differentia 1 ! stiffness 

The value of the determinant of the sum of the elastic stiffness and the 
differential stiffness Is automatically printed for each differential stiffness 
loading condition. 

Iterative differential stiffness computations are terminated for one of five 
reasons. Iteration termination reasons are automatically printed In an information 
message. These reasons have the following meanings: 

1. REASON 0 means the Iteration procedure was Incomplete at the time of exit. 
This Is caused by an unexpected interruption of the Iteration procedure prior to 
the time the subroutine has had a chance to perform necessary checks and tests. 

Not much more has happened other than to Initialize the exit mode to REASON 0. 

2. REASON 1 means the iteration procedure converged to the EPS 1 0 value 
supplied by the user on a PARAM bulk data card. (The default value of EPS 10 
Is 1.0E-5.) 

3. REASON 2 means iteration procedure is diverging from the EPS 1 0 value 
supplied by the user on a PARAM bulk data card. (The default value of E P S 1 0 
Is 1.0E-5.) 

4. reason 3 means Insufficient time remaining to achieve convergence to the 
CPS 1 0 value supplied by the user on a PARAM bulk data card. (The default value 
of EPS10 is 1.0E-5.) 

5. REASON 4 mtans the number of Iterations supplied by the user on a PARAM 
bulk data card has been met. (The default number of Iterations Is 10.) 

Parameter values at the time of exit are automatically output as follows: 

1. Parameter DUNE: -1 Is normal; ♦ N Is the estimate of the number 

of Iterations tequired to achieve convergence. 
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2. Parameter SHIFT: +1 Indicates a return to the top of the Inner loop 

was scheduled; *1 Indicates a return to top of the outer loop was scheduled 
following the current Iteration. 

3. Parameter DSEPSI: the value of the ratio of energy error to total 

energy at the time of exit. 
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3.22.4 Cast Control Deck DTI Table and Parameters for Static Aerothermoelastic 
Analyst! with Differential Stiffness 

The following Items relate to subcase definition and data selection for 
Static Aerothermoelastic Analysis with Differential Stiffness: 

1. The Case Control Deck must contain two subcases. 

2. A static loading condition must be defined above the subcase level with 
a LOAD, TEMPERATURE(LOAO), or DEFORM selection, unless all loading is 
specified by grid point displacements on SRC cards. 

3. An SPC set must be selected above the subcase level unless all constraints 
are specified on GRID cards. 

4. Output requests that apply only to the linear solution must appear in the 
first subcase. 

5. Output requests that apply only to the solution with differential stiff- 
ness must be placed in the second subcase. 

6. Output requests that apply to both solutions, with and without differen- 
tial stiffness may be placed above the subcase leve. 

7. Aerodynamic input for the Aerodynamic Load Generator (ALG) module is 
input via data block ALG08. This data block must be input using Direct 
Table Input (DTI) bulk data cards. For a detailed description of the 
ALGD6 data block Input see Section 1.15.3.1 of the User's Manual. 
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The following output may bo requested for Static Aerothermoelastlc Analysis with 

Differential Stiffness: 

1. Nonzero Components of the applied static load for the linear solution at 
selected grid points. 

2. Displacement and nonzero components of the single-point forces of con- 
straint, with and without differential stiffness, at selected grid 
points . 

3. Forces and stresses In selected elements, with and without differential 
st 1 f f ness . 

4 . Undeformed and deformed plots of the structural model. 

I. 6RDPNT - optional - a positive Integer value of this parameter will cause the 
Grid Point Weight Generator to be executed and the resulting weight and 
balance Information to be printed. 

•• WTMASS - optional - the terms of the mass matrix are multiplied by the real 
value of this parameter when they are generated In ENG. 

I. IRES - optional - a positive Integer value of this parameter will cause the 
printing of the residual vectors following the execution of SSG3. 

k CflUPMASS - CPBAR. CPR0D, CP0UA01 . CPQUAD2 . CPTRIA1 . CPTR1A2. CPTUBE. CPQDPLT. 


CPTRPLT, CPTRBSC - optional - these parameters will cause the generation of 
coupled mass matrices rather than lumped mass matrices for all bar elements, 
rod elements, and plate elements that Include bending stiffness. 

BETAD - optional- the Integer value of this parameter Is the assumed number of 
Iterations for the Inner loop In shift decisions for Iterated differential 
stiffness. The default value Is 4 Iterations. 

NJ[ - optional - the Integer value of this parameter limits the maximum number 
of Iterations. The default value Is 10 Iterations. 
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7. IFSII - optional • the real valut of this paras* tar la used to tost th* con* 
vergenc* of Iterated differential stiffness. The default value It 10"*. 

S. APRESS - optional In static aerothermoelastlc analysis. A positive Integer 
value will generate aerodynamic pressures. A negative value (the default) 
will suppress the generation of aerodynamic pressure loads. 

9. ATEMP - optional In static aerothersoelastlc analysis. A positive Integer 
value will generate aerodynamic temperature loads. A negative value (the 
default) will suppress the generation of aerodynamic thermal loads. 

10. STREAMl - optional In static aerothermoelastlc analysis. STREAML*1 causes 

the punching of STREAML1 bulk data cards. STREAML ■ 2 causes the punching of 
STREAML2 bulk data cards. STREAML-3 causes both STREAML 1 and STREAML2 cards 
to be punched. The default value, -1, suppresses punching of any cards. 

11. PGE0H - optional In static aerothermoelastlc analysis. PGE0M«1 causes the 
punching of GR10 bulk data cards. PGE0M-2 causes the punching of GRID, 

CTRIA2 and PTR1A2 bulk data cards. PGE0M a 3 causes the punching of GR10 cards 
and the modified ALGDB table on DTI cards. The default, *1, suppresses punch* 
1 tg of any cards . 

12. 1 1 > RT - optional In static aerothermoelastlc analysis. If I PRT > 0, then 
Intermediate print will be generated In the ALG module based on the print 
option In the ALGOB data table. If IPRT ■ 0 (the default), no Intermediate 
print will be generated. (IPRTCI, IPRTCL, IPRTCF) 

13. SIGN - optlona' In static aerothermoelastlc analysis. Controls the type of 
analysis being performed. SIGN ■ 1.0 for a standard analysis. SIGN * *1.0 
for a design analysis. The default 1$ 1.0. 

14. Z0RIGN, FXC0PR, FYCQQR, FZCOQR - optional In static aerothermoelastlc 
analysis. These are modification factors. The defauls are Z0RIGN >0.0, 
FXCOOR • 1.0, FVCOOR - 1.0, and FZCOOR • 1.0. 
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KTOUT - optional In static eerotheraoelastlc analyses. A positive Integer 
of this parameter Indicates that the user wants to save the total stiffness 
Matrix on tape (6IN0 file INPT) via the OUTPUT! Module In the rigid format. 
The default 1$ *1. 


3 . 22- 27 ( 9 / 30 / 78 ) 


COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 

3.23 COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 

3.23.1 DMAP Sequence For Compressor Blade Cyclic Mod<1 Flutter Analysts 

KICII) FORMAT U.NAP LISTING 
SERIES 0 

AERO APPROACH, AIGIU FORMA I 9 

LEVEL 2.0 NASIRAN DMAP COMPILER - SOURCE LISTING 


OPTIONS IN EFFECT* GO ERR-2 NOLI ST NOOECK NOREF NOUS CAR 

1 BEGIN AERO NO. 9 COMPRESSOR BLAOE CYCLIC MOOAL FLUTTER ANALYSIS A 

2 FILE PHIML-APPENO/AJJL'APPENJ/F SAVE* APPENO/CASEVV* APPEND/ CLAMAL" 

A»'P£NO/UVG*APPtNJ/gHHL«APPENO A 

3 GE39l,GtOM2,/GPL,E0t XI N,GPO T ,C S 1 H ,DGPOT , SI L/ V • N,LUS ET/ V,N, 
NOGPiJT A 

A SAVE LUSET.NOGPOT A 

5 C gONU ^ ERROR I.NOGPOT S 

6 CHLPN T GPL , LOL X IN ,GPOT ,C STM ,BGPO T ,SI L A 

7 PURGE OIJE, 02 JE /NO) JE * 

8 GLCM2,LUEXIN/ECT A 

9 CHRPNT CCT A 

0 < feP3 3 jEUMi,LJEXlN»GE0M7/#GPT!/V«N ,NOGR AV A 

11 ChRPNT GPTT $ 

12 ECTtEPI , UGPOT «SIL,v»PTT,CSTM/EST , GE I ,GPEC T ,/ V , N, LUSET/ V,N, 

NO SI HP /C ,N,i/ V ,N, KOGE NL/V, A, GENE L A 

NOGENL.NOSIMP ,GENEL A 

ERRJK I iNOSIMP A 

UGPSI/GENtL A 

EST.GPLC T.GEI .OGPST A 

//C»fl»ADD/V»'i tNUKGo X/C ,N,I /C , N, 0 A 

//C,N « AOO/VtN , NOMf.G/C ,N,l/C ,N,0 A 

// C ,N , NOP / V, V, KGG I N»-l A 

JMPKGG IN , RGGI N A 

/ /C , N • AOl) / V, N »NORGG X /C,N,-1 /C.N.O A 
/ KGGIN, ,,,/: ,V,LOCATION»“l /C . V , I NPTUNI T*0 A 



3.23-1 (9/30/78) 


RIGID FORMATS 


RIGID FORMAT OMAR LISTING 
SERIES 0 

ACRU APPROACH, RIGID FORMA T 9 

LEVEL 2.0 NASIKAN OMAP COMPILER - SOURCE LISTING 


23 

2A 

25 

26 


CgouTo 


CMKP*0 

LABEL 



2? 

2R 

29 

30 

31 

32 

33 
3A 

35 

36 

37 
3d 
39 
AO 
A 1 
A 2 
A3 
AA 
AS 
A6 


SAVE 

CFUPNT 



CHKPNT 

LABEL 



LAbfcl 

(^uTT) 

CWPUI 


<£ED> 

Chkpni 
L ABEL 

cgZ^> 


KGGIN.KGGX A 
KGGX A 
JHPKGCtN A 

E Sf ,CSTM,*PT,r.lT,GEOM2,/KtlM,fcl)ICt ,MELM,MOICT , , / V ,N, NOKGGX/ V 
N,NU.1GG/C,N,/C ,N, /C,N,/C.V,CUUPMA$S/C, Y,CPUAH/C*Y»CPKUIV l. V 
CPJU4Jl/C,Y,:PJUftO<VC,V,CPJRI At /C,V,CPTA1A2/C,V tCPtU DC/ C,V. 
CPODPC t/C, r,:PIRPl l/C, V.CPTAtfSC a 

NOKGGX.NUMGG A 

K ElM *KJ |C T ,ME LM,HU1C 1 A 

JMPKGGX i N0KG3 X $ 

f 

GPEC l *KJIC I iXELN/KGG XtGP ST A 
KGGX.GPST A 
JMPKGGX A 
EKHUfU.NUMGG A 

GPEC I,MU1C1,MEIM/MGG ,/C ,N,-I /C, V ,WT MASS* 1.0 I 
MGG A 

LGPMG.GKOPNT t 

BGPUl ,C SIM.EJE KIN ,MGG/UGPWG/V ,Y , G«OPNt >«-l/C , V , W! MAS S t 
OGPrfG,, ,,,// A 
LGPMG A 

k’cgx , KGO/NOGE NL I 
KGG A 

t at 1 1 , NOGENL $ 

GEI.KGCX/KGG/ V,N,I USE 1 / V , N , NOGF NL/ V «N ,NOS IMP $ 

KGG A 
l Bl 1 1 A 

C ASECC • GFUM*, EQE XI M, SU.GPOT ,BGPl)I,CSf P/ KG, , USE T, ASET/ V,N, 
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RIGID FJRM4? D9RP L T S 1 1 No 
SERIES 0 

AERO APPROACH? RIGID FORMS I 9 

LEVEL 2.0 NAS ERA N OMAP COMPILER * SOURCE LISTING 


Hi 

40 

49 

50 

51 

5 2 

53 

54 

55 

56 

57 
50 

59 

60 
61 

6 2 
63 
69 
fcb 
66 

67 

68 

69 

70 


LOSE t/V,N,*p;f l/V,N,MPLT//V,N,SI NGLE/V,N,OMIT/V,N,AEACT/C,N,0/ 
V , 9 * R LP E A f / V, N , NO St I / V.N ,NOL/ V, N ,NOA/t , T ,SUd(0 I 

N'C F I , S INGLE , OKI T ,ft£ AC T, NOSE t ,H'CF 2 .AC '£ AT , ACL » NOA A 

//C«N,NJ1/V,N,fttACi>ATA /V,N, REACT A 

EHAORS.HLACuA TA A 

GH,GNO/NPCT I/GO, GOO/ONI T/KFS«QPC /SINGLE A 
GEOM<,«EOEXtN, USE! /CKO/ V.V.CIV'E / V,N,NOG0 A 
NUGO A 
CVCO A 

LHk JR 6,'40G0 A 
L»L*,uE4EL A 

GPL, OP ST, USE I ,SIL/03PST/V,N,N0GPSI A 

NOGPSJ A 

L Ul 6 ,NOGP ST A 

OGPSI,,, ,,// A 

L8L9 A 

RGG.KNN/HPCEl /NGG ,HNN/MPCFl A 

KNN, INN A 

L8L/.MPCH A 

USE 1 ,KG/GN A 

ON A 

USE I , GM , KGG ,MGG , , /KNN ,NNN, , A 
KN.N,NNN A 
L8L / A 

KNN.KFF/SINGLE/MNN.MFF/SINGLE a 
K FT * M F F A 
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rigid Fork ai diap listing 

SERIES 3 

AEMU APPROACH, RIGID FDR HR I 9 

LEVEL 2. 0 NASIRAN DHAP COMPILER • SODRCE LISTING 



BA 

65 

66 
87 



88 

B9 

90 

91 


SAVE 

ChKPNT 



92 SAVt 

93 CHIP 9 T 


L DL i , SINGLE » 

USET,RNN,MNN, ,/RFF ,RF S,,HFF,, » 

RFF,RFS,HfF A 
L DL 3 6 

RIF kAA/OMIT/ HFF.HAA/GHIT I 
K AA , rt AA A 
l BL 5,l)H| T A 

USET ,RFF,,,/GU«RAA ,K(0 ,LOU, , , • • $ 

GO , K A A A 

0 3ET,GU,NFF/HA4 a 
9 AA A 
L BL 5 A 

0Y9AH ICS, GPL, SI l , USE l/GPLO.SI LO , USE 1 0 , TF POOL, , , , , « ECO, ECOVN/V • 
N,lUStt/V,N,LUi>t TD/V.N.NUTI l / V, N , NUOL I / V »N» NUPS OL/V,N,NOFRL/V# 
N.M'JNLF l/V,N,NUTRL/V«N*NUEtO/C,N«/V , N, NODE A 

LUSE I J.NUUfc ,NUUO A 

ERR JR 2,N0tEU A 

GO , GJD/.90UE /J H ,GND /NCUE A 

C VCO.KAA ,HAA, , , /RKR..9KR,,, / C.N.FORE / V,V,NSEGS-~I /V,V, 

R INOE X«- 1 / V,Y,C YCSEQ*-l / C,N,l / w.N.NOGO A 

NUGD A 

RRK.MRR A 
ERRJK6.N0G0 A 

KRR ,HRR, , ,EEJ « ,CA SECC / LAHR,PMIR, .OEIGS / C,N, MOOES /V,N, 
NbIGV A 

NEIGV A 

L AMR ,PH|R, HI, OEIGS A 
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COMPRESSOR BLADE CYCLIC NODAL FLUTTER THESIS 

AICIO FORMAT ONAP LISTING 
SERIES 0 

AERO APPROACH. RIGID FORMAT 9 

LEVtL 2.0 NASTRAN OMAP COMPILIK - SOURCE LISTING 


44 

PAR 44 

*5 


94 

SAVt 

47 


40 

CCvuD 

49 

SAVt 

100 

CHK^Nf 

lot 

(gjHU 3 

ua 

<sso 

103 

(4^0 

104 


104 

SAVfc 

106 

<£pcu^) 

107 

SAVE 

106 

CnKPN 1 

104 

CPZmtT) 

110 

<^mp7ao} 

111 

(JjTkTTn) 

ill 

SAVE 

tli 

PURGE 

114 


116 

CMCPNr 

1 16 

(g^ 


/'C.N.MPV / V.N.CAMONO / C.N.O / C.N.O A 
J E I G S .L AMR . • . . // V.N.CARONO A 
CARONO > 
gRAU* A.NCIGV A 

CVCO. ...PMIK.LAMK /» » ,P1l A .LAMA / C.N.BACK / V.V.NSEGS /V.V . 

R IN JL X / V.V.CVCSEU /C.N.l / V.N.NUGO A 

S0G3 A 

lama.phia » 

ERROR t , N03C> A 

USE I. .PillA. «, Cu ,GH, , XF i, . / PHlG.t / C.N.l / C.N.RCIG A 

CASECC,CSTl«1PT,0IT,E0EXlN,SIL,,«BGPDT .LAMA. . PHIG.EST , , / .. 
OPHIG... / C.N.REIG A 

OPHIG, ,,,«// V.N.CARONU A 

CARONO A 

EOT . USE T ,BGPJ T.CSIP.EOEXIN.GM.GO / AERO, ACPI , Ft 1ST . GT KA. PVECT / 
V.N,NK/V,N,NJ/V,Y„1I NMALH/V, Y.MAXMACH/V. V, IREF/V, V.MTVPE/V, M* 
NEIoV/V.V.KlNOEX«-l A 

NK.N J A 

AERO.ACPT.Fll ST.GTKA, PVECT A 
PHIA, PVECT, / PHI A X, , , / C.N.l A 

PHIAX.1AA.PHI AX,, , / Ml / C , N ,3 /C . N , I / C , N, 1 / C. N ,0/C , N, 1 A 

C ASfcCC , MA TI’UJL .EQ0VN,,TFPUCL/K2PP,M2PP,82PP/V,N,LUSETD/V,N, 
N0K2PP/V,N,N0M2PP/V,N,NUtf2PP A 

NOX2PP,NOM2PP .NU82PP A 

X 2D0/N<K2PP/M2OD/N0M2PP/B2OO/NOB2PP A 

N2PP ,.120D/*JSE T/B2PP.H2JO/NOSET/K2PP.K2DO/NOSET A 

K2PP.M2PP.U2PP ,K2DU,N200,B200 A 

USE TO, CM .CO ,, , « ,K2PP .M2PP.B2 PP/ , , ,GMJ,G00,K200, M2 00, 8200/C, N, 
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HIGIO FORMAT 0*AP LISTING 
SCRIES U 

AERO APPROACH, RICH) FORMA T 9 

LEVEL 2.0 NasTRAN OMAP COMPILER - SOURCE LISTING 


CMPLEV/CiN.OI SP/C ,N,MO0AL/C,N,0.0/C,N,0.0/C,N,0.0/V,N,NOK2PP/V, 
N«N0M2PP /V«N« NU62PP/ V,N, HPCF 1/V»N*SINGLE/V, N,ON IT /V » Nt NOOE/ C, 
N»-1/C,N*-1/C , N« - 1/C *N»-l t 


I IT CKKPIII * 200, M2DD, 82)0, GOU.UMJ » 

1 18 USE T J«PM|A A,M I iLAMKfUl T , M200 ,S2D0,K200,CASECC / MMM • 8MM • RMM . 

" PtllOH / V, N , NUUC /C , V , l MODE S*VM99‘)V/C,V,LFXt©*0.0/C«V« MFKEO* 0.0/ 
V.N,NUM2PP/V,N,NOB2PP/V.N.NOR2PP/V,N,NONCUP/V#H.FHOOt/C.V, 

KDAMP — l % 

1 19 SAVE NUNC OP « f MJQE » 

120 CHKPN1 MMM,<1HM,K»M,PH|DH $ 

121 PAKANL PCOB//C,N,PRt S/C , N , /C ,N, /C ,N , /V , N , NUPCOO » 

122 PURGE PL1SE I* , PL TPAK ,GP SE I S.LLSE IS / NOPCOt) » 

123 ( fUN J P2. 4UPCUII % 

12* ( PL 'iSEp PCU8, CGOYN.E: T / PLISE U.PLIPAR, UPSETS, ELSETS / V.N.NSILl /V,N, 
JUMP PLOT »- 1 % 

125 SAVE N SIL 1, JUMPPLOT % 

126 ( gaiHS^ PETSEI* // * 

127 PARAM //C,N,MPV/V,N .PLTFLG/C ,N,1 /C,N,1 » 

128 PAR AN //C « N ,MP Y/ V«N «PF 1 L£ /C «N,0/C ,N,0 I 

129 ( CUNO^ ) P2, JUMPPLOT t 

l 30 ( £lu7~^ ) PL TP AR , uP St r S .ELSE TS <0 A St C C • BGPO f , EQ1 YN , , , , ,/ PLOT X 1/ V,N« NS I L 1/ 
V.V.LUSE T/V.N ,JUMPPLUI/V,N,PLlFLG/V,N,PF ILE A 

131 SAVE JUMPPLOT, PL TFLG,PFILt % 

132 (£•<][*>£) rLUTRI // % 

133 L ABEL P2 % 

135 < guNO~^ > ERROR 2,NUEE0 $ 

135 PARAM //L,J,A00/V,N,0ESTRY/C,N,0/C,N,1 » 

1 36 Ca ro aeru.acp t/ajjl ,skj.oijk,o2 jk/v,n,nk/v,n, nj/v,n, oesi rv » 

137 SAVE OESTRY 1 
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RIGID EOANAT DM AP LISTING 
SERIES U 

AERO APPROACH, RIGID FORMA T 9 

LEVEL 2.0 NASIKAN OMAP COMPILER - SOURCE LISTING 


18 CMLPNT 
39 ( gONO 3 
AO ( Input? /) 

A I LABEL 
A2 P ARAM 

A3 ( Imp 




59 SAVE 

60 CCONU 


AJJL ,SKJ,01J* ,02 JK A 
NOJJ E «NO0 JC $ 

/Ul J C , 02JE , • • /C ,V,POSI Tl ON— l/C ,V,UNS TNUM-1 1/ C»V *OSRl ABEL* 

T AP t I 0 A 

NOUJC A 

//C , N , ADO/ V ,N , XQHH L/C ,N,l /C.N,0 A 

AJJL , SAJ ,01 JK ,0 2JK,G TKA , PH I OH ,01 JE , 02 JC , US ET 0, AERO/ GHHL, QJHL/V , 
N.NOUt/V.N.XJHHL A 

XQilHL A 

OHHL ,QJML A 

//t,N, IP V/V,N ,NOP/C ,N,-l /C , N , 1 A 
//C, M,MP Y/V.N ,NOP/C,N,l/C ,N,I A 
/ /C , N ,MP V/ V ,N , NOH /C , N ,0/C , N , 1 A 
//C,.\I,MPY/V*N,F LOOP/ V, V, NOD JE — l /C, N,0 A 
LOUPTUP A 
LUUPTUP A 

K MH, JHM, P- 'H , JHML , C A SFCC ,FLI ST/E SAVE ,KXHH ,BXHH, MS — V V ,N, FLOOP/V , 
N.IStARI A 

FLCIUP , TSTAKT A 

KXHH, BXHH,MXHH ,FEO ,C A SEC C /PHI H , C LAMA, OCE IGS/ V ,N ,E IOVS A 
tlGVS A 

l UL Z AP , E IGVS A 
LUL 16.NUH A 

CASECC , EQOYN, USE TO ,PH | H ,C LAMA , , /CPH l H , /C ,N, CE IGEN/ C, N, MOOAL /C, 
N, 1 2 3/V, N, NUM /V»N,NOP/V,N,F MODE A 

NOH, NOP A 

LBL 16 «NOH A 
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RtCIO FURMAT OMAP 1 1 STING 
SERIES U 

AERO APPROACH, RIGID FUR HA T 9 

LEVEL 2,0 NASTRAN OMAP COMPILER - SOURCE LISTING 


161 

16 2 SAVE 
U3 LAUEL 
164 ( fA2 


165 

166 
167 
166 
16V 

170 

171 

172 

173 

174 

175 

176 



LABEL 



177 
176 
179 
16 J 
iQl 
182 



183 CTcbT 

184 CHKt'NT 


□ PHI H ( , i , , / /V iN «C AHDNO $ 

CAR UNO I 
L Ul 1 6 % 

PMIH,CLAMA,ESAvE/PH|ML,CLA HAL ,C A St Y V ,OVG/ V , N,T S T ART /C, Y ,VREF« 
1. 0/C, V, PH I NT- YE SB * 

I STAR! » 

PMIHL,CLANAL,CASE YY.OVO » 

CONTINUE, I STAR 1 i 
L BE ZAP i 

CONTINUE ,F LOOP % 

LOOP TUP ,100 $ 

ERROR 3 % 

CONTINUE A 
OVG * 

XVCOU//C.N, PUS/C tN,/C ,N,/C,N*/V,N»NOXYC0B * 

NOXYUU ( , NO XYCOrt i 

X YCJli ,OVU, , , , /XYPLTCE/C.N, VG /C , N ,PSE T/ V, N, PF ILE/V , N , C ARONO > 
PI IL E ,C ARUNO » 

XVPL ICE// t 
NOXYOLI t 

//C,.\.,ANO/V,N,PJUMP/V,N.NOP--l/‘,N, JJMPPLOT S 
F IM I S , P J UMP » 

CASE Y Y,CLAMAL ,PHI HL ,C ASECC , , /CL AMALl , CPM I HI , CAS EZZ, ,/C, N, 

C E I ON $ 

CPH I H 1 , P H I OH/ C PHI 0 » 

CPMIU V 
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RIOIO FORMAT UNAP LISTING 
SCRIES 0 

ACRU APPROACH. A IG 10 FORMAT 9 

LEVEL 2.0 NASTRAN OMAP COMPILER - SOURCE LISTING 


165 

(fuuiQ 

CPM|U,CPM1P/N0A $ 

166 

CfcUNU 2) 

LUL 14 t .M0A $ 

187 

esc Hp) 

OSETOt tCPHIl), , 'GOOfGMDt St t /C PHI P , tQPC/C • Nt 1/ C, N, DYNAMICS % 

IBB 

l ADEL 

l Bl 1 4 % 

16V 

CHKPNT 

CPMIP.JPC $ 

190 

<mID 

CP Ml l) t CP Hi A /i^ UUC $ 

191 


LHLNOfc tNCUf I 

19 2 

<$Tc 

USETl)/BP/CtNtO/CtNtA/C tN.E 1 

193 

<£AklQ 

CPMlD.tHP/CPHlA,, f /C «Mtl/C»N« 3 * 

m 

LABEL 

LbtNUL* $ 

195 


CASEZ2 V CSTM«MPT»3I TtEgOVNtSllOt » tBGPDT tCLAHAll « QPC » CPHl P» EST# » / 
tUOPCl.UCPMlP ,UESC1 tUE/Cl t PC PHI P/C»NtCFI GN * 

196 

CHKPNT 

PCPH1P * 

19/ 


UCPHIPtUOPClt UfcSC ItUEFCl # . // V tN #CARDNO S 

196 


P 3i J UMPPLOT A 

199 


PLTPAR tGP! C TS «E L SE 1 S »C A Sk Z Z t BGPO T tEQDYNt S 1 L0» » PCPHI P« »/PlO! X 3/ 
VtNtNSILl/ x *IUSE T/V»N. JUMP PLOT / V tM# PL TF LG/V# N* PULE $ 

200 


PLUTXJ// $ 

201 

L ABLE 

PJ I 

202 

(JUMP 

MN1S A 

203 

LABEL 

ERROR 1 A 

204 

(tjRTPAK^) 

//Ci.'4»- 1/CtNtP SUB SON $ 

205 

LABEL 

ERROR 2 A 

206 

igjUP AKJ$ 

//Ctl.*2/CtN»FSUBSCN $ 

20/ 

LABEL 

ERROR 3 A 

206 

<0RTPAR?$ 

/ /C # M »- 3 /C » N# P SUB SCN $ 
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RIGID FORMAT UMAR LISTING 
SERIES 0 

AERO APPROACH, RIGID FORMAT 9 

LEVEL 2.0 NASIRAN OMAP COMPILER - SOURCE LISTING 


209 

L AHtL 

ERROR A % 


210 


//C,.4,” A/C 

t N« F SUB SliN % 

211 

LABEL 

ERRORS % 


212 

(^rT7ak$> 

// C »N A 

/ : » n tC vcMoofc s % 

2 1 i 

LABEL 

ERROR 6 % 


214 


U C.N.-S 

f Z # NfC YCMOOES $ 

21* 

LAB1L 

F IN 1 S i 


216 

ENC 

% 



• *^U ERRORS FOUND - EXECUTE NASTRAN PROGRAM** 
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3,23 

3 . 

S. 

8 . 

10 . 

12 . 

14. 

20 . 
21 . 

22 . 

23. 

26. 

29. 

30. 

33. 

34. 

36. 

37. 

38. 

40. 

42. 

43. 

46. 

49. 

51. 

54. 


Description of DMAP Operations for Compressor Blade Cyclic Modal 
Hut ter Analysis" 

GP1 generates coordinate system t ransf orma t ion matrices* tables of grid 
point locations, and tables for relating internal and external grid 
point numbers. 

Go to DMAP No. 203 and print error message if no grid points arc present. 

GP2 generates Element Connection Table with internal indices. 

GP3 generates Static Loads Table and Grid Point Temperature Table. 

TA1 generates element tables for use in matrix assembly and stress 
recovery . 

Go to DMAP No. 203 and print error message if no elements have been 
defined. 

Go to DMAP No. 2b if stiffness matrix is not user input. 

Set parameter NOKGGX * -l so that the stiffness matrix will not 
be generated in DMAP No. 26. 

INPUTT1 reads the user supplied stiffness matrix from tape (G1N0 
file I N P T ) . 

Equivalence [K* g ] to [K gg ]. 

EMG generates structural element matrix tables and dictionaries for 
later assembly. 

Go to DMAP No. 32 if no stiffness matrix is to be assembled, 

EMA assembles stiffness matrix [K*^] and Grid Point Singularity Table. 

Go io DMAP No. 203 and print error message if no mass matrix exists. 

EMA assembles mass matrix [M^]. 

Go to DMAP No. 39 if no weight and balance request. 

GPWG generates weight and balance information. 

flFP formats weight and balance information and places it on the 
system output file for printing. 

Equivalence [K* ] t0 C g 3 IF no general elements. 

Go to DMAP No. 45 if no general elements. 

SMA3 adds general elements to t g ^ t0 obta * n stiffness matrix [K^]. 

GP4 generates flags defining members of various displacement sets 
(USET), forms multipoint constraint equations [R^Hu^} « 0. 

Go to DMAP No. 211 and print error message if free-body supports 
are present. 

GPCYC prepares segment boundary table. 

Go to DMAP No. 213 and print error message If CYJOIN data is Inconsistent. 
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RIGID FORMATS 


Go to DMAP No. 60 If general elements present. 

GPSP determines If possible grid point singularities remain. 

Go to DMAP No. 60 if no grid point singularities remain. 

RFP formats the table of possible grid point s ingular i ties and places 
it on the system output file for printing. 

Equivalence [K 3 to [K ] and [Ml to [Ml if no multipoint 
constraints. 99 99 n 

Go to DMAP No. 68 if MCE1 and MCE? have already been executed for 
current set of multipoint constraints. 

MCE 1 partitions multipoint constraint equations [R g ] « ! R n J and 

solves for multipoint constraint transformation matrix [G^] * 

MCE? partitions stiffness and mass matrices 


1 V * 

K ! K 

nn i nm 

1 

1 

and 

r V - 

*n„ 

^n»" 


*ran * Snm 



v 

M mm 

j 


and performs matrix reductions 


IS,,,! • <"<1 

t"n„J ' ‘«n„l * 

Equivalence I^ nn 3 to [K^3 and [M nn ] to [M^] if no single-point 
const ra 1 n ts . 

Go to DMAP No. 74 If no single-point constraints. 

SCE1 partitions out single-point constraints. 


' K ff i 
* 

K fs' 

and 


"(f ! 

1 

M fs‘ 

— - —f- - 

i 


f 

1 


V ! 

V 



> ! 

M ss_ 


Equivalence [K^j to [K aa ] and to [M aa J if no omitted degrees 

of freedom. 

Go to DMAP No. 8? if no omitted coordinates. 
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f 

! 


78 . 


80 . 


83 . 


8S . 


86 . 


87 . 
90 . 
91 . 


95 . 


97 . 


98 . 


SMP 1 pirtlttons constrained stiffness matrix 


tKff] 


H 


ao 




L oa * 


OOj 


and solves for transformat ion matrix [G 0 ] * -T K oo ]'* [ K q ^] 

and performs matrix reduction [K ] ■ [K ] ♦ f 3f G 1 . 

oo 04 0 

SMP? partitions constrained mass matrix 


ft 


tMffJ 


aa i 




ao 


oo 


and performs matrix reduction 


(»..] ■ t"..J 4 I "!,”'. 3 4 4 


DPD generates flags defining members of various displacement sets used 
In dynamic analysis (USETD), tables relating internal and external 
grid point numbers, Including extra points introduced for dynamic 
analysis, and prepares Transfer Function Pool and Eigenvalue Extraction 
Data . 


Go to DMAP No, 205 and print error message If no Eigenvalue Extraction 
Data . 

Equivalence [G o J to [G^] and [G m ] to [G^j] if no extra points Introduced 
for dynamic analysis. 

CYCT? transforms matrices from symmetric components to solution set. 

Go to DMAP No. 213 and print error message if CYCT2 error was found. 
READ extracts real eigenvalues from the equation 

t‘kk - • 0 • 

and normalizes eigenvectors according to one of the following user 
requests : 

1) Unit value of selected coordinate 
2} Unit value of largest components 
3) Unit value of generalized mass. 

(JFP formats eigenvalues and summary of eigenvalue extraction infor- 
mation and places them on the system output file for printing. 

Go to DMAP No. 209 and exit if no eigenvalues found. 

CYCT2 finds symmetric components of eigenvectors from solution set 
eigenvec tors . 
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101 . 
10 2 . 


103. 

104. 

105. 
109. 

no. 

Ill . 
114. 

116. 


Go to DMAP No. 213 and print error message if CYCT2 error was found 
SDR1 recovers dependent components of the eigenvectors 



SDR? prepares eigenvectors for output (RPHIG). 

9FP formats tables prepared by SDR? and places them on the system 
output file for printing. 

APOD processes the aerodynamic data cards from EDT. AERO and ACPT 
reflect the aerodynamic parameters. PVECT is a partitioning vector 
and GTKA is a transformation matrix between aerodynamic (K) a.:d 
structural (a) degrees of freedom. 

PAR7N partitions the eigenvector into all sine or all cosine 
components . 

SMPYAD calculates modal mass matrix 

r«j • r.;j' [.„] iO 

MTRXiri selects the direct input matrices [K pp ] t ^ M pp3 • t c pp ]- 

Equivalence [H* p ] to [C pp ] to [B^] and [K* p ] to [K d(j ] if no 

no cons t ra 1 n t s applied. 

2 2 

GKAD applies constraints to direct input matrices C K pp 3 • [ M pp J* an< * 
2 2 

and [ Bj ^ 3 (see Section 9.3.3 of the Theoretical Manual) and 
forms [G md ] and [G Qd ]. 


3.23-14 ( f '/30/78) 


COMPRESSOR BLADE CYCLIC MODAL FLUTTER ANALYSIS 


1 16 . 


123. 

124. 

126. 

129 . 

130. 
132. 

134 . 

136. 

139. 

140. 


GKAM selects eigenvectors to form and assembles stiffness, matrices 


and damping matrices 

in modal coordinates: 




• n-f« • 




■ N-:S] • 

^dh^ M dd^*dh^ ’ 



• fr-iS] • 

£*dh^ B dclJt4 dh ] 

where 




KDAMP - 1 


KDAMP - 

-1 (defaul t) 


* modal masses m^ « modal masses 

bj ■ m < 2* f |9<f | ) b < • 0 

k i • m < 4«* f. k 4 * { 1 ♦ 1 g ( f 1 ) ) 4«*fjm < 

Go to OMAP No. 133 If no plot package Is present. 

PLTSET transforms user input into a form used to drive structure 
plotter. 

PRTMSG prints error messages associated with structure plotter. 

GO to OMAP No. 133 if no undeformed aerodynamic structure plot request. 

PL0T gonerates all requested urdeformed structure plots. 

PRTMSG prints plotter data and engineering data for each undeformed 
aerodynamic plot generated. 

Go to OMAP No. 205 and print error message if no Eigenvalue 
Extraction Data. 

AMG forms the aerodynamic materix list £ A j j ] , the area matrix [S^], 

and the downwash coefficients [Dj^] and [0*^]. 

Go to OMAP No. 141 if no user- supp 1 i ed downwash coefficients. 

1NPUTT2 provides the user-supplied downwash factors due to extra 
points ([Dj e ] . 
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143. 


U9. 

150. 

151. 

152 . 


AMP computes the aerodynaMlc matrix list related to the model 
coordinates as follows: 


t*dhJ • 

♦al [ '.e 

♦el 1 ♦ee 

U J 



£G k iJ 

• t G Ie 

f°5h> - 

t»jl l °jiJ 


10)0 

* f°jk 

t°jhJ - 

i°j. i V 1 


10)0 

* lc U 

For each 

(m,k) pair: 






* 

t°jP * <«°jy 



for each 

group: 







■ 


group 



’ki 


fw ■ I s t jHQj h J 
[Q^h^ ■ C**n^tQ^h3 


ki • 


•kh J 




*ih 


PARAM initializes the flutter loop couter (F199P) to zero. 

Go to next DMAP instruction If cold start or modified restart. 
190PT9P will be altered by the Executive System to the proper 
location inside the loop for unmodified restarts within the loop. 

Beginning of loop for flutter. 

FA1 computes the total aerodynamic mass matrix the total 

aerodynamic stiffness matrix and the total aerodynamic 

damping matrix as well as a looping table FSAVE. For 

the K-inethod 

'* « U a /b’)M hh ♦ (p/25 Q> 


M 


hh 


<hh 


'hh 


'hh 


'hh 
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154. CEAO extracts complex eigenvalues fro* the equation 

Kh?' ♦ C" ‘ 0<V • o 

and normal 1 its eigenvectors to unit magnitude of largest component. 

1S6. Go to OMAP No. 166 If no complex eigenvalues found. 

167. Go to DMAP No. 163 If no output request for the extra points Intro* 
duced for dynamic analysis or modal coordinates. 

168. VDR prepares eigenvectors for output, using only the extra points 
Introduced for dynamic analysis and modal coordinates. 

160. Go to DMAP No. 163 If no output request for the extra points 
Introduced for dynamic analysis or modal coordinates. 

161 0FP formats eigenvectors for extra points Introduced for dynamic 

analysis and modal coordinates and places them on the system output 
file for printing. 

164. FA2 appends eigenvectors to PHI HI , eigenvalues to CLAMAL, Case Control 
to CASE Y Y f and V-g plot data to 0VG . 

167. Go to DMAP No. 172 If there Is Insufficient time for another 
flutter loop. 

169. Go to DMAP No. 172 If flutter loop complete. 

171. Go to DMAP No. 207 for additional aerodynamic configuration triplet 

values . 

175. Go to DMAP No. 179 if no X-Y plot package Is present. 

176. XYTRAN prepares the input for requested X-Y plots. 

178. XYPL0T prepares requested X-Y plots of displacements, velocities, 

accelerati ons , forces , stresses , loads or single-point forces 
of constraint vs. time. 

181. Go to DMAP No. 215 if no output requests Involve dependent degrees 
of freedom or forces and stresses. 

182. M0DACC selects a list of eigenvalues and vectors whose Imaginary parts 
(velocity in Input units) are close to a user input list. 

183. DDR1 transforms the complex eigenvectors from modal to physical 
coordinates 

t.$] • ■ 

185. Equivalence [#£] to [$£] If no constraints applied. 

186 Go to OMAP No. 188 if no constraints applied. 
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187 . 


190 . 

191 . 

192 . 

193 . 


195 . 

197 . 

198 . 

199 . 

200 . 

202 . 

204 . 

206 . 


SORI recovers dependent coaponents of eigenvectors 


‘♦d> 









<♦» 


♦ ♦') 
n T e' 


<*i> • ♦ ♦;> • 



snd recovers single-point forces of constraint (q $ J • 

j’q’ | ■ dp • 

Equivalence [♦£] to If no extra points Introduced for dynamic 

analysis . 


Go to OMAP No. 194 If no extra points present. 

VfcC generates a d-slze partitioning vector (RP) for the a and e sets. 
PARTN performs partition of [♦£] using RP, 

<♦£> ■ 

SDR? calculates element forces and stresses (0EFC1, 0ESC1) and 
prepares eigenvectors and single-point forces of constraint for 
output (0CPMIP, 0QPC1). It also prepares PC PH I P for deformed plotting. 



OFP formats tables prepared by SDR2 and places them on the system 
output file for printing. 


Go to DMAP No. 194 If no deformed structure plots are requested. 
PLOT prepares all deformed structure plots. 

PRTMSG prints plotter data and engineering data for each deformed 
plot generated. 

Go to DMAP No. ?15 and make normal exit. 


M0OAL COMPLEX EIGENVALUE ANALYSIS ERROR MESSAGE NO. 1 - HASS HATRIX 
REQUIRED F0R MODAL FORMULATION. 

MODAL COMPLEX EIGENVALUE ANALYSIS ERROR MESSAGE NO. ? - EIGENVALUE 
EXTRACTION DATA PEQUIRED F0R REAL EIGENVALUE ANALYSIS. 
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208. MODAL COMPLEX EIGENVALUE ANALYSIS ERROR MESSAGE NO. 3 
EXECUTE MORE THAN 100 LOOPS. 

210. MODAL COMPLEX EIGENVALUE ANALYSIS ERROR MESSAGE N0. 4 
VALUES REQUIRED FOR MODAL FORMULATION. 

212. NORMAL MODES MITH CVCLIC SYMMETRY ERROR MESSAGE NO. 4 
SUPPORTS NOT ALLOMED. 

214. NORMAL M0OES MITH CVCLIC SYMMETRY ERROR MESSAGE NO. S 
SYMMETRY OATA ERROR. 


- ATTEMPT TO 
• REAL EIGEN- 

- FREE B0OV 

- CYCLIC 
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3.23.3 Output for Compressor Blade Hodaj Flutter Analysts 

Thp Real Eigen value Summary T«P1c and the Real Eigenvalue Analysis 
tuaaary. as described under Normal Mode Analysts, ara automatically prmtaa. 

Alt raal algenvaluas are Included even though all may not ba used in the 
Modal formulation. 

The grid point singularities from the structural model ara also output. 

A flutter summary for each value of the configuration parameters is 

printed out if PRINT»YESB. This shows p, k, l/k,a, ff *v V, g and f 

sound* * 

for each complex eigenvalue. 

V-g and V-f plots may be requested by the XVRUT control cards bv 
specifying the c«»*. ve type as VG. The "points” are loop numbers and the 
•components* are G or F. 

Printed output of the follow! tg types, sorted by complex eigenvalue 
root number (SORT I ) and (m, k, p) may bt requested for all complex eigenvalues 
kept, as either real and Imaginary parts or magnitude and phase angle 
(0* - 360" lead): 

1. The eigenvector for a list of PHYSICAL points (grid points, 
extra points) or SOLUTION points (modal coordinates and extra points). 

2. Nonzero components of the single-point forces of constraint for 
a list of PHYSICAL points. 

3. Complex stresses and forces In selected elements. 

The RFREOUENCY case control cord can select a subset of the complex eigenvectors 
for data recovery. In addition, undeformed and deformed shapes mav be reouested. 
Undeformed shapes may include only structural elements. 
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Case Control Deck and Parameters for Compressor Blade Cyclic 


Modal Flutter Analysis 

1. Only one subcase Is allowed 

2. Desired dvrect Input matrices for stiffness [Y.^ppJ, mass lM*pp], and 
damping [B ] must be selected via the keywords K2PP, M2PP, or 
B2PP. 

3. CMETH0D must be used to select an EIGC card from the bulk Data Oeck. 

4. FMETH0D must be used to select a FLUTTER card from the Bulk Oata Oeck 

5. NETH0D must be used to select an EIGR card that exists in the Bulk 
Oata Deck. 

6. SOAMPING must be used to select a TABDMP1 table if structural damoina 
is desired. 

7. An SPC set must be selected unless the model is a free body or all 
constraints are specified on GRID cards. Scalar Connection Cards or 
with General Elements. 

B. Each NASTRAN run calculates modes for only one symmetry index, K* 

The following user parameters are used in Compressor Blade Cyclic Hodal 

Flutter Analysis. 

1, GRDPNT - optional - A positive integer value of this parameter will 
cause the Grid Point Weight Generator to be executed and the result- 
ing weight and balance information to be printed. All fluid related 
masses are ignored. 

?. WTHASS - optional - The terms of the structural mass matrix are 
multiplied by *he real value of this parameter when they are 
generated in SMA2. Hot recommended for use In hvdroelastic 
problems . 

3. C0UPMASS - CPBAR . CPR0D, CP0UAD1 , CP0UAD2 , CPTRIA1. CPTRIA2, 


CpTUBE . C PQDPLT » CPTRPLT. CPTJIBSL - optional - These parameters 
will cause the generation of coupled mass matrices rather than 
lumped mass matrices for all bar elements, rod elements, and plate 
elements that include bending stiffness. 
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LFRCQ and HFREQ - required unless L HIDES Is used. The reel values 
of these parameters give the frequency range (LFREQ Is lower 
limit and HFREQ Is upper limit) of the modes to be used In the 
modal formulation. To use this option, LMflOES must be set to 0. 

L M|DF $ - used unless set to 0. The Integer value of this parameter 
Is <:he number of lowest modes to be used In the modal formulation. 
The defalult value will request all modes to be used. 

N0OJE - optional In modal flutter analysis. A positive Integer 
of this parameter Indicates that user supplied downwash matrices 
due to extra points are to be read from tape via the 1NPUTT2 
module In the rigid format. The default value Is *1. 

Pi , P2 and P3 - required In modal flutter analysis when using 
N0DJE parameter. See Section 5.3.2 for tape operation parameters 
required by INPUTT2 module. The defaults for PI, P2, and P3 are 
-1, 11 and TAPE1D, respectively. 

VREF - optional In modal flutter analysis. Velocities are 
divided by the real value of this parameter to convert units 
or to compute flutter Indices. The default value Is 1.0. 

PRINT - optional In modal flutter analysis. The BCD value N0, 
of this parameter will suppress the automatic printing of the 
flutter summary for the k method. The flutter summary table 
will be printed If the BCD value Is YES for wing flutter, or 
YESB for blade flutter. The default Is YES. 

CT YPE - required - the BCD value of this parameter defines 
the type of cyclic symmetry as follows: 

(1) R|T - rotational symmetry 

(2) DHL - dihedral symmetry, using right and left halves 

(3) OSA - dihedral symmetry, using symmetric and anti- 
symmetric components 

NSCGS - required - the Integer value of this parameter Is the 
number of Identical segments In the structural model. 
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12. CYCSEQ - optional - tht Integer value of this parameter specifies 
the procedure for sequencing the equations In the solution set. 

A value of +1 specifies that all cosine terms should be sequenced 
before all sine terms, and a value of -1 for alternating the 
cosine and sine terms. The default value Is -1. 

13. KINDEX - required In compressor blade cyclic modal flutter analysis. 
The Integer value of this parameter specifies a single value of 

the harmonic Index. 

14. Ml NMACH - optional In blade flutter analysis. This Is the minimum 
Ma;h number above which the supersonic unsteady cascade theory Is 
val Id. The default Is 1.01. 

15. MAXMACH - optional In blade flutter analysis. This Is the maximum 
Mach number below which the subsonic unsteady cascade theory Is 
valid. The default value Is 0.80. 

16. IREF - optional In blade flutter analysis. This defines the 
reference streamline number. IREF must be equal to a SLN on a 
STREAML2 bulk data card. The default value, -1, represents the 
streamsurface at the blade tip. If IREF does not correspond to 
a SLN, then the default will be taken. 

17. MTYPE - optional In cyclic modal blade flutter analysis. Tills 
controls which components of the cyclic modes are to be used In 
the modal formulation. MTYPE * SINE for sine components and 
MTYPE * COSINE for cosine components. The default BCD value Is 
COSINE. 

18. KGGIN - optional In blade flutter analysis. A positive Integer 
of this parameter Indicates that the user supplied stiffness 
matrix Is to be read from tape (GINO file INPT) via the 
INPUTT1 module In the rigid format. The default Is -1. 
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4.149 FUNCTIONAL MOOULE ALG (AERODYNAMIC LOAD GENERATOR) 

4.149.1 Entry Point : UD0300 

4.149.2 Purpose 


The principal function of ALG Is to generate an aerodynamic pressure and/or 
temperature distribution for compressor blades. The ALG module may also be used 
as a compressor blade mesh generator to punch GRID. CTRIA2 and PTR1A2 bulk data 
cards. Bulk data cards STREAML1 and STREAML2 can also be generated by ALG by 
user request. 

4.149.3 DMAP Calling Sequence 

ALG CASECC, EDT, EQEXIN, {Jjj®}), ALGDB , CJTM, BGPDT / CAJECCA, GE0M3A/ 

S, Y, APREJJ/ %, Y, ATEMP/ V, Y, STREAML/ V. Y, PGE0M/ V, Y, IPRT/ 

S. N. 1FAIL/ V, Y, SIGN/ V, Y. Z0R1GN/ V, Y. FXC00R/ V, Y, FYC00R/ 

V. Y, F2C00R $ 

4.149.4 Input Data Blocks 


CASECC 

EDT 

EQEXIN 

AUGV . 
UBGV ’ 

ALGDB 

CSTM 

BGPDT 


• Case control data table 

• Aerodynamic bulk data cards 

• Equivalence between external grid or scalar numbers and Internal 
numbers 

- Displacement vector matrix giving displacements In the g-set 

- Compressor blade data table 

- Coordinate system transformation matrices 

- Basic grid point definition table 


Notes: 


1. CASECC and ALGDB cannot be purged. 

2. AUGV or UBGV can be purged. 
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FUNCTIONAL MODULE ALG (AERODYNAMIC LOAD GENERATOR) 

3. EQEXIN, CSTM and BGPDT can be purged If AUGV Is purged. 

4. EOT can be purged If AUGV Is purged and parameter STREAML • *1. 

5. ALGDB may be Input via DTI bulk data cards. 

4.149.5 Output Data Blocks 

CASECCA - Revised case control data table 

GE0M3A - Static load and temperature table 

Note : 

1. CASECCA and GE0M3A may not be purged. 

4.149.6 Parameters 

APRESS - Input - Integer - default * -1. If APRESS > 0. then aerodynamic 
pressures will be generated. 

ATEMP • Input - Integer - default ■ -1. If ATEMP > 0, then aerodynamic 

temperatures will be generated. 

STREAML - Input - Integer - default * -1. Controls the punching of STREAML 1 

and STREAML2 cards. STREAML • 1, punch STREAML 1 cards. STREAML ■ 2, 
punch STREAML 2 cards. STREAML • 3. punch both STREAML 1 and STREAML2 
cards. 

PGE0M - Input - Integer - default * -1. Controls the punching of blade 

geometry bulk data cards. PGE0M ■ 1, punch GRID cards. PGE0M » 2, 
punch GRID, CTRIA2 and PTRIA2 cards. PGEOM • 3, punch GRID cards 
and the modified ALGDB table on DTI cards. 

IPRT - Input - Integer - default * 0. If IPRT > 0, then Intermediate 
print will be generated based on the print option in ALGDB data 
table. 
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I FA I L - Output - Integer - default * 0. Set to -1 If there Is a convergence 

fat lure . 

SIGN - Input - real - default * 1.0. Controls the type of analysis 

being performed. SIGN • 1.0 for standard blade analysis. 

SIGN ■ -1.0 for design analysis. 

I0RIGN - Input - real - default ■ 0.0. Modification factor. 

FXC00R - Input - real - default ■ 1.0. Modification factor. 

FYC00R - Input - real - default • 1.0. Modification factor. 

FZC00R - Input - real - default • 1.0. Modification factor. 

4.149.7 Method 

(a) Data block AIGDB contains all the Input needed to generate the aerodynamic 
pressures and tei f eratures on the compressor blade. However, the aero- 
dynamic loads are a function of the blade shape and the data defined In 
ALGDB must first be modified to account for any change in the blade shape 
or Input via the displacement vector matrix AUGV . If AUGV is purged, 
then ALGDD is not modified. The AIGDB data block is read and the aero- 
dynamic loads are calculated for the compressor blade being analyzed. 

(b) The CASECC data block is read and a copy of it is output to CASECCA with 
changes to data items 4 and 7 for all subcases. In CASECCA, word 4 is 
set to 60 If aerodynamic pressure loads were generated, and word 7 is 
set to 70 if aerodynamic thermal loads were generated. 

(c) The GE0M3A data block contains aerodynamic load and temperature data. 

If parameter APRESS > 0, then PL0AD2 cards with set Identification number 
60 are stored on GE0M3A. If parameter ATEMP > 0, then TEMP and TEMPO cards 
with set identification number 70 arc stored on GF0M3A 
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(d) Parameters STREAMl and PC COM control the punching of buH dr', a cards 

STREAML1 , STREAML?, GRID, CTRIA2, PTRIA2 and DTI. The ALG module may be 

used In a one module DMAP program as a compressor blade mesh and 
geometry generator as follows: 

BEGIN % 

ALG CASECC... .ALGDB, ./CASE CCA, GE0M3A/C . N , - 1 /C .N.-l/C.N ,3/C ,N ,2/C ,N,1$ 

END ) 

4.149.6 Subroutl nes Called 

4.149.8.1 Utility subroutines GMMATS, PRETRS and TRANSS are called. 

4.149.8.2 Subroutine Name: UD03PR 

1. Entry Point: UD03PR 

2. Purpose: Modify ALGDB data block. 

3. Calling Sequence: CALL UD03PR (1ERR) 

4.149.8.3 Subroutine Name: UD03PB 

1. Entry Point: UD03PB 

2. Purpose: Identify data fields as being either BCD alpha, real or integer. 

3. Calling Sequence: CALL UD03PB ( 1 DAT , NTYPE) 

4.149.8.4 Subroutine Name: UDO3P0 

1. Entry Pcint: UDO3P0 

2. Purpose: Generate data blocks CASECCA and GE0M3A. 

3. Calling Sequence: CALL UDO3P0 (SCR1) 
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4. 149. 8. 5 Subroutine Name: 


U003AP 


1 . Entry Point: U003AP 

2 . Purpose: Punch the modified AlGDB table data block on 

OTI Bulk Oata cards if parameter PGtflM - 3. 

3. Calling Sequence: CALL UD03AP (IFNAME, IFNfl) 
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UD03AN 

UDQ313 


U003AR 

UDQ3Q 1 -UDQ3Q7 


UD0314 
UD0315 
UD03U 
U0031 7 
UD03U 
UD0319 


UD0309-UD0312 
UD0325 
U0D329 
UD0330 
UDG1-U06 9 

Z = 1 


UDQ308 


UD0326 
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may occur: 3001, 3002, 3003 and 3008. 
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4. ISO FUNCTIONAL MODULE APDB (AERODYNAMIC POOL DISTRIBUTOR FOR BLADES) 

4.150.1 Entry Point : APDB 

4.150.2 Purpose 

Bulk data cards which control the solution of aerodynamic problems are 
processed and assembled Into various blocks for convenience and efficiency in 
the solution of the aerodynamic problem. APDB also generates the transformation 
natrlx [Gj ta ] T (GTKA) and the partitioning vector PVECT. 

4.150.3 DMAP Calling Sequence 

APDB EOT, USET , BGPDT, CSTM, EQEXIN, GM, G0/ AER0, ACPI, (“LIST, GTKA, 

PVECT/ V. N, N K / V* U , NJ/ V, Y, MINMACH/ V, Y , MAXMACM/ V, Y , IRE 
V, Y, MT Y PE/ V, N, NEIGV/ V, Y, K1NDEX ■ -1 $ 

4.150.4 Input Data Blocks 

EDT - Aerodynamic bulk data cards 

USET - Displacement set definition table 

BGPDT • Basic grid point definition table 

CSTM - Coordinate system transformation matrices 

EQEXIN - Equivalence between external points and scalar Index values 
GM - Multipoint constraint transf ormati on matrix 

G0 - Structural matrix partitioning trans f o rma tl on matrix 

Notes : 

1. EDT. USET. BGPDT and EQEXIN cannot be purged. 

2. CSTM may be purged If all points are In the basic system. 
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3. GM and G0 may be purged If there ere no multipoint or no omitted 
points . 

4.150.5 Output Data Blocks 

AER0 - Control Information for control of aerodynamic matrix generation 
and flutter analysis 

ACPT - Information pertaining to each Independent group of aerodynamic 
elements 

FLIST - Contains AER0, FLFACT and FLUTTER cards copied from EDT 

GTKA - Aerodynamic transf ormat Ion matrix - K set to a set 

PVECT - Cyclic modes partitioning vector for matrix PHIA from module CYCT2 

Notes : 


1 . 

AER0 , ACPT, 

FLIST and GTKA cannot be 

purged . 


2 . 

PVECT may be 

purgeJ If there are no 

cyclic modes 

to be partitioned 

4.150.6 

Paraneters 




NK 

- Output 

- Integer - no def aul t . 

Degrees of 

freedom In the NK 


displacement set . 

NJ - Output - Integer - no default. Oegrees of freedom In the NJ 

displacement set. 

MINMACH - Input - real - default * 0.8. This Is the maximum Mach number 
below which the subsonic unsteady cascade theory Is valid. 

HAXMACH - Input - real - default ■ 1.01. This is the minimum Mach number 
above which the supersonic unsteady cascade theory Is valid. 
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IRCF • Input - Integer • default • This defines the streamline number 

of the reference stream surface. IRCf must equal an SLN on 
a STREAML2 card. The default value, *1, represents the stream 
surface at the blade tip. If IRCF does not correspond to an SLN, 
then the default will be taken. 

MYTPE - Input - BCD - default • COSINE . This controls which components of 
the cyclic modes are to be used In the modal formulation. MTYPE * 
SINE for sine components and MTYPE • COSINE for cosine components. 

NEIGV * Input - BCD * no default. The number of eigenvalues found. 

Usually output by the READ module. 

K1NDEX - Input • BCD - default • -1. Harmonic Index number used in cyclic 
analyses . 

4. ISO. 7 Method 


Subroutine APDB Is the main control program for this module. It allocates 
buffers, reads input files, and Initializes output files. APDB creates the AERO, 
ACPT and FLlST tables and generates the PVECT partitioning vector. Subroutine 
APDB 1 generates the GTKA transformation matrix. APDB1 reduces [G^] t0 


mu 


ch like module SSG2, using the following matrix operations: 


- 


KN 


G T 

KM 


K» J • < c m’’ Km’ * "In’ 




*Kf 


*KS 
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ic’,) . 

“I.) • < 6 / K.J • 

At each step where a matrix multiply Is Indicated, the multiply It skipped If the 
result Is known to be aero (l.e., U fl or are null). 

4. 150. 8 Subroutines Called 

Utility routines BISLBC, CALCV, SSG2B, TRANSS and GMMATS all called. 

4. ISO. 8. I Subroutine Name: APDB1 

1. Entry Point: APOBl 

2. Purpose: To generate transformation matrix [gJ ]. 

3. Calling Sequence: CALL APDB 1 (JBUM. IBUF2, NEXT, LEFT, NJTNI , NL1NEI, 

LCJTlS, ACITM, N0DEX, N0DE1, IJILC, XY2B). 

4. 150. 9 Design Requirements 

Open core Is located at /APDB2Z/. APDB uses five scratch files. 

4.150.10 Diagnostic Messages 

System fatal messages 3001, 3002, 3003, 3008 and 3037 may occur. The APDB 
module also generates Its own messages that are not numbered. These messages are 
sel f-e*plana tory . 
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6,1,1 % 16 Algid Format Error Messages for Stw ! c Atrothermoelast ic Analysis with 
Differential Stiffness 

Hp * 1 - flp STRUCTURAL ELEMENTS HAVE BEEN DEFINED. 

The differential stiffness matrix is null because no structural elements 
have been defined with Connection cards. 

NP. 2 - FREE BpDY SUPPORTS N0T ALLOWED . 

Free bodies are not allowed in Static Analysis with Differential 
Stiffness. The SUPPORT cards must be removed from the Bulk Data Deck 
and other constraints applied if required for stability. 

NO l • MASS MATRIX REQUIRED FOR WEIGHT AND BALANCE INSULATIONS. 

The mass matrix is null because either no elements were defined with 
Connection cards, nonstruc tural mass was not defined on a Property card, 
or the density was not defined on a Material card. 

MP. 5 - N0 INDEPENDENT DEGREES OF FREEDOM HAVE DEEN DEFINED. 

Either no degrees of freedom have been defined on GRID, SPOILT or Scalar 
Connection cards, or all defined degrees of freedom have been constrained 
by SPC, UPC, OMIT, or GRDSET cards, or grounded on Scalar Connection cards. 
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6.1.3 Aero Approach Rigid Formats 

The texts of the rigid format error messages arc given In the following 
lection for the aero approach rigid formats. The text for each message Is given 
In capital letters and Is followed by additional explanatory material, including 
suggestions for remedial action. 

6. 1.3.1 Rigid Format Error Messages for Modal Flutter Analysis 

N0 , 1 • MASS MATRIX REQUIRED F0R M0DAL F0RMULAT I0N . 

The mass matrix Is null because either no structural elements were 
defined with Connection cards, nonstructural mass was not defined 
on a Property card or the density was not defined on a Material card, 

NO. 2 - E IGl.. VALUE EXTRACTION DATA REQUIRED F0R REAL EIGENVALUE ANALYSIS 

Eigenvalue extraction data must be supplied on an EIGR card and METHOD 
must select an EIGR set In the Case Control Deck. 

NO. 3 - ATTEMPT TO EXECUTE MORE THAN 100 LOOPS. 

An attempt has been made to use more than 100 different sets of direct 
input matrices. This number can be increased by altering the REPT 
Instruction following FA2. 

NO. 4 - REAL EIGENVALUES REQUIRED FOR MODAL FORMULATION. 

No real eigenvalues were found In the frequency range specified by the 
user . 

6. 1.3. 2 Rigid Format Error Messages for Compressor Blade Cyclic Modal Flutter 
Analysis. 

NO. 1 - MASS MATRIX REQUIRED FOR MODAL FORMULATION 

The mass nutrix is null because either no structural elements were defined 
with Connection cards , nons true tura 1 mass was not defined on a Property card 
or the density was not defined on a Material card. 

NO. 2 - EIGENVALUE EXTRACTION DATA REQUIRED FOR REAL EIGENVALUE ANALYSIS 

Eigenvalue extraction data must be supplied on an EIGR card and METHOD 
must, select an EIGR set in the Case Control Deck. 

NO. 3 - ATTEMPT TO EXECUTE MORE THAN 100 L00PS. 

An attempt has been made to use more than 100 different sets of direct 
Input matrices. This number can be Increased by altering the RfPT instruc- 
tion following F A2 . 

NO. 4 - REAL EIGENVALUES REQUIRED F0R MODAL FORMULATION. 

No real eigenvalues were found in the frequency range specified by the 
user. 
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i - FREE B0DY SUPP0RTS H0T ALL0WEO. 

Free bodies are i'iot allowed In Statics with Cyclic Symmetry* The SUP0RT 
cards must be removed from the Bulk Data Deck and other constraints 
applied If required for stability* 

► « CYCLIC SYMMETRY OATA ERR0R. 

See Section 1.12 for proper modeling techniques and corresponding PARAM 
card requl rements . 
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A 


P 

Parameter value used to control utility module 
HATGPR print of A-set matrices. 

ABFL 


DBM 

rA b * Hydroelastic boundary area factor matrix. 

ABFLT 


DBM 

Transpose of fA b ^1 

ACCE 


IC 

Abbreviated form of ACit; RATION. 

ACCELERATION 

IC 

Output request for acceleration vector. (UM-2.3, 4.2) 

ACPI 


DBT 

Aerodynamic Connection and Property Data. 

Active 

Col umn 

PH 

Column containing at least one nonzero term 
outside the band. 

ADD 


FMrt 

Functional module to add two matrices together. 

ADD 


M 

Parameter constant s d in utility module PARAM. 

ADD5 


FMM 

Functional Module to add up to five matrices together. 

ADUM1 


IB 

Defines attributes of dummy elements 1 through 9. 

AEFACT 


IB 

Used to input lists of real numbers for 
aero elastic analysis. 

AER0 


DBT 

Aerodynamic Matrix Generation Data. 

AER0 


IB 

Gives basic aerodynamic parameters. 

AJJL 


DBML 

Aerodynamic Influence Matrix List. 

ALG 


FMS 

Aerodynamic load generator. 

ALGDB 


DBT 

Aerodynamic Load input for ALG (D-16). 

ALL 


IC 

Output request for all of a specified type of output. 

ALLFDGL 

TICS 

IC 

Request tic marks on all edges of X-Y plot. 

ALOAD 


P 

Set negative if no aerodynamic loads (0*16) . 

ALTER 


IA 

Alter statement for DMAP or rigid format. 

ALWAYS 


P 

Parameter set to -1 by a PARAM statement in the 
Piecewise Linear Analysis Rigid Format ( D * 6 ) . 

AMG 


FMS 

Aerodynamic Matrix Generator 

AMP 


FMS 

Aerodynamic Matrix Processor, 

AND 


M 

Parameter constant used in executive module PARAM. 

ACHJTS 


M 

Indicates restart with solution set output request. 

APD 


FMS 

Aerodynamic pool distributor and element generator. 
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APDB 

APP 

APPEND 

APRESS 

ASET 

ASET1 

ATEMP 

AUT0 

AUT0 


NASTRAN DICTIONARY 

FMS Aerodynamic pool distributor for blades. 

1A Control card which specifies approach (DISP or DMAP). 

H Fiel may be extended (see FILE). 

PU Positive Value generates aerodynamic pressures. 

IB Analysis set coordinate definition card. 

IB Analysis set coordinate definition card 

PU Positive value generates aerodynamic tempera tures , 

IC Requests X-Y plot of autocorrelation function. 

DRT Autocorrelation function table. 
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CQ0MEM2 

IB 

Quadrilateral membrane element connection definition card. 

CQDPLT 

IB 

Quadrilateral bending element connection definition card. 

CQUAD1 

IB 

General Quadrilateral element connection definition card. 

CQUAD2 

IB 

Homogeneous quadrilateral element connection definition card. 

CRfO 

16 

Rod element connection definition card. 

CSNEAR 

IB 

Shear panel element connection definition card. 

CSL0T3 

18 

Triangular slot element connection definition card for acoustic 
analysis. 

CSL0T4 

IB 

Quadrilateral slot element connection definition card for 
acoustic analysis. 

CSP 

tc 

Selects a set of contact surface points. 

CSP 

IB 

Contact surface point set definition. 

CSTM 

OBT 

Coordinate System Transformation Matrices. 

CSTMA 

DBT 

Coordinate System Transformation Matrices - Aerodynamics. 

CTETRA 

IB 

Tetrahedron clement connection definition card. 

CTfRORG 

IB 

Toroidal ring element connection card. 

CTRAPRG 

IB 

Trapezoidal ring element connection card. 

CTRBSC 

IB 

Basic bending triangular element connection definition card. 

CTRIA1 

IB 

General triangular element connection definition card. 

CTRIA2 

IB 

Homogeneous triangular element connection definition card. 

CTRIARG 

IB 

Triangular ring element connection card. 

CTRMEM 

IB 

Triangular membrane element connection definition card. 

CTRPLT 

IB 

Triangular bending element connection definition card. 

CTUBE 

IB 

Tube element connection definition card. 

CTWIST 

IB 

Twist panel element connection definition card. 

CURVIINESVMBDI 

IC 

Request to connect points with lines and/or to use symbols 
for X-Y plots. 

CVISC 

IB 

Viscous damper element connection definition card. 

CMEDGE 

13 

Wedge element connection definition card. 
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F#RCE2 

IB 

Static load definition (magnitude and four grid points). 

FfRCEAX 

IB 

Static load definition for conical shell problem. 

FREEPT 

IB 

Defines point on a free surface of a fluid for output purposes. 

FREQ 

IB 

Frequency list definition. 

FREQS 

N 

Indicates restart with change In frequencies to be solved. 

FREQ1 

IB 

Frequency list definition (linear Increments). 

FREQ? 

IB 

Frequency list definition (logarithmic Increments). 

FREQRESP 

P 

Parameter used In SDR2 to Indicate a frequency response problem. 

FREQUENCY 

IC 

Selects the set of frequencies to be solved In frequency 
response problems. 

FRL 

DBT 

Frequency Response List. 

FRQSET 

P 

Used in FRRO to Indicate user selected frequency set. 

FRRD 

FMS 

Frequency and Random Response - Displacement approach. 

FSAVE 

DBT 

Flutter Storage Save Table. 

FSIIST 

IB 

Defines a free surface of a fluid In a hydroelastic problem. 

Functional Module 

PH 

An Independent group of subroutines that perform a structural 
analysis function. 

FXCOOR 

PU 

Aerodynamic modification factor (D-16). 

FYCOOR 

PU 

Aerodynamic modification factor (D-16). 

FZCOOR 

PU 

Aerodynamic modification factor (D-16). 
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1C 


10 

I PAIL 
IFP 

I FP1 

IFP3 

IFP4 

IKAG 

lhr- L 
I NCI UDC 
INERTIA 

INERTIA RrUET 
1NPT 

INPUT 

Input Data Block 

Input Data Cards 

INPUTT) 

INPUT! 2 
i NPUTT 3 
INPUTT4 
Internal Sort 

INV 

I PRT 
IREF 
IRES 


IC Transient analysis initial condition set selection. 

IA The first card of any data deck is the identiftcaiu n (10) 

card. Tne two data ite.us on this card are BCD values. 

P Set negative by ALG if convergence fails (0-16). 

EH Input File Processor. The preface nodule which processes 

the sorted Bulk Data Dock and outputs various Jata Mocks 
depending on the card types present in the Bulk Data Deck. 

EM Input File Processor 1. The preface nodule which processes 

the Case Control Deck and writes the CASECC, PCOB and IYCDB 
data blocks. 

EM Input File Processor 3. The preface module which processes 

bulk data cards for a conical shell problem. 

EM Input File Processor 4. The preface module which processes 

bulk data cards for a hydroeldstic problem. 

IC Output request for real and imaginary parts of some quantity 

such as displacement , load, single point force of constraint 
element force, or stress, 

P Parameter constant used in executive module PARAM. 

IC Used in set definition for structure plots. 

P Used in printing rioid format error messages for Static 

Analysis with Inertia Relief (0-2). 

IA Selects rigid format for static analysis with inertia relief. 

M A reserved NASTRAN physical unit (Tape) which must be set 

up Dy the user wnen used. 

FMU Generates most of bulk data for selected academic problems. 

PH A d 3 ta block input to a module. An input data block must 

have been previously output from some module and may not be 
written on. 

PH The card input data to the NASTRAN system are in 3 sets, the 

Executive Control Deck, the Case Control Oeck, and the Bulk 
Data Deck. 

FMU Reads data blocks from GIN0-written user tapes. 

FMU Reads data blocks from F(JRTRAN-wri tten user tapes. 

FMX Durmjy user input module. 

FMX Duw\y user input module. 

PH Same order as external sort except when SEQGP or SEQEP bulk 

data cards are used to change the sequence. 

IB Inverse power eioenvaluc analysis option - specified on 

LIGR, EIGG or EIGC cards. 

PU Controls printing of aerodynamic results, 

PU Defines reference streamline for blade flutter. 

P Causes printout of residual vectors in statics rigid formats 

when set nornegative via a PARAM bulk data card. (0-1, 0-2, 
0-4. D-S, 0-6). 
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KOSS 

OBH 

[Kj s ] - Partition of differential stiffness Matrix. 

KFF 

DBM 

[Kff] * Partition of stiffness matrix. 

ICFS 

DBM 

[F fs 3 - Partition of stiffness matrix. 

KGG 

OBM 

[K 1 - Stiffness matrix generated by Structural Matrix 
« Assembler. 

KGGIN 

PU 

Positive value selects KGGX from INP'JTTl. 

KGG IN 

DBM 

Sum of elastic and differential stiffness 
matrices (D-16, A-9). 

KGGL 

DBM 

(K 4 ] - Stiffness matrix for linear elements. Used only in 
99 the Piecewise Linear Analysis Rigid Format (0-6). 

KGGLPG 

P 

Purge flag for KGGL matrix. If set to -1, It implies that 
there are no linear elements in the structural model. (D-6). 

KGGNl 

DBM 

[k" 4 ] - Stiffness matrix for the nonlinear elements. Used 
*9 i.i the Piecewise Linear Analysis Rigid Format only. 

{0-6). 

KGGSIH 

DBM 

Stan of KGGNL and KGGL. Used in the Piecewise Linear Analysis 
Rigid Format only. (D-6). 

KGGX 

DBM 

[K*g] - Stiffness matrix excluding general elements. 

KGGXL 

DBM 

[K ] - Stiffness matrix for linear elements (excluding 
99 general elements). Used In the Piecewise Linear 

Rigid Format only. (0-6). 

KHH 

DBM 

C*C- h ] - Stiffness matrix used in modal formulation of 
n dynamics problems (D-10 thru D-12). 

KLL 

DBM 

[K f# ] - Stiffness matrix used in solution of problems in 
“ static analysis (0-1, 0-2, D-4, D-5, D-G). 

KLR 

DBM 

[K 3 - Partition of stiffness matrix. 

(CNN 

DBM 

LK ) - Partition of stiffness matrix. 

m 

OBM 

[K 3 * Partition of stiffness matrix. 

KRR 

DBM 

[K pr 3 - Partition of stiffness matrix. 

ICSS 

DBM 

[K s$ 3 - Partition of stiffness matrix. 

KTOTAL 

DBM 

Sum of elastic and differential stiffness 
matrices (D-16, A-9) . 

KTOUT 

PU 

Postlve value outputs KTOTAL to 0UTPUT1. 

KXHH 

DBM 

Total modal stiffness matrix - h-set. 
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HHTT5 

IB 

Specifies table references for temperature-dependent, 
anisotropic, thermal material properties. 

MAX 

IB 

Eigenvector normalization option • used on EIGR, EIGB and 
EI6C cards. 

MAXIMUM OE FORMATION 

IC 

Indicates scale for deformed structure plots. 

NAXIT 

P 

Limits maximum nuofcer of Iterations In nonlinear heat transfer 
analysis. 

MAXLINES 

IC 

Maximum printer output line count - default value Is 20000. 

NAXMACH 

PU 

Controls subsonic unsteady cascade calculations. 

MCE 1 

FMS 

Multipoint Constraint Eliminator • part 1. 

MCE? 

FMS 

Multipoint Constraint Eliminator - part 2. 

MOO 

DBM 

[M,J - Mass matrix used In direct formulation of dynamics 
0fl problems (D-7 thru 0-9) . 

MOEMA 

P 

Parameter Indicating equivalence of MOO and MAA. 

MOLCEAD 

P 

Used In printing rigid format error tuessages for modal 
complex eigenvalue analysis (D-10). 

MDLFRRD 

P 

Used in printing riqld format error messages for modal 
frequency response (0-11). 

MOLTRD 

P 

Used in printing rigid format error messages for modal 
transient response (D-12). 

MERGE 

FMM 

Matrix merge functional module. 

METHOD 

IC 

Selects method for real eigenvalue analysis. 

METHODS 

M 

Indicates restart with change In eigenvalue extraction 
procedures . 

MFF 

DBM 

[M^] - Partition of mass matrix. 

MGG 

DBM 

[Mgg] - Mass matrix generated by Structural Matrix Assembler. 

MHH 

DBM 

[M hh ] - Mass matrix used in modal formulation of dynamics 
problems (D-10 thru D-12). 

MI 

DBM 

[m] - Modal mass matrix. 

MINMACH 

PU 

Controls supersonic unsteady cascade calculations. 

MKAER01 

IB 

Provides table of *ich numbers and reduced frequencies (k). 

MKAER02 

IB 

Provides list of Mach numbers (m) and reduced frequencies (k). 

HL 

OBM 

[M lt ] - Partition of mass matrix. 

MLR 

DBM 

[M lr ] - Partition of mass matrix. 

MNN 

DBM 

[Ml - Partition of mass matrix, 
nn 

MOA 

DBM 

[Ml - Partition of mass matrix. 

Od 

MOOA 

FMX 

This module is reserved for user Implementation. 

MOOACC 

FMS 

Mode Acceleration Output Reduction Module. 

HTYPE 

PU 

Controls cyclic mode component selection. 
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P 

Packed Format 

PAPER SIZE 
PARAM 
PARAM 
Parana tar 

PARAM. 

PARAMR 

PARTN 

P8AR 

P8L 

PBS 

PCD3 

PC0NEAX 

POAMP 

PDF 

POT 

PDUM1 

PELAS 

PEN 

PENSIZE 

PERSPECTIVE 

PFILE 

P6 

PC 

PG1 

PGEOM 

PGG 


IC 

FMU 

IB 

PH 


FMU 

FMU 

FMM 

IB 

OBM 

DBM 

OBT 

IB 

IB 

OBM 

DBM 

18 

IB 

IC 

IC 

IC 

P 

DBM 

DBM 

DOM 

PU 

DBM 


Parameter valua used In MAT6PR to print P-set matrices. 

A matrix Is said to be In packed format If only tlia nontaro 
elements of the matrix are written. 

Selects paper size for structure plots using table plotters. 
Performs specified operations on ONAP parameters. 

Parameter definition card. 

A FfRTRAN variable comttunl cated to a OMAP module by the NASTRAN 
Executive System through blank coomon. A parameter's position 
In the OMAP calling sequence to a module corresponds to the 
position of the parameter In blank common at module execution 
time. 

Selects parameters from a user input matrix or table. 

Performs specified operations on real or complex parameters. 
Matrix partitioning functional module. 

Bar property definition card. 

A scalar multiple of the PL load vector. Used only In the 
Differential Stiffness Rigid Format (D-4). 

A scalar multiple of the PL load vector. Used only In the 
Differential Stiffness Rigid Format (0*4). 

Plot control data block (table for use with structure plotter 
functional module PLTSET). 

Conical shell element property definition card. 

Scalar damper property definition card. 

Dynamic load matrix for frequency analysis. 

Linear dynamic load matrix for transient analysis. 

Property definition card for duoiny elements 1 through 9. 

Scalar elastic property definition card. 

Selects pen size for structure plots using table plotters. 

Selects pen size for X-Y plots using table plotters. 

Specifies perspective projection for structure plots. 

Parameter used by PL (IT module. 

Incremental load vector used In Piecewise Linear Analysis 

( 0 - 6 ). 

Statics load vector generated by SSG1. 

Static load vector for Piecewise Linear Analysis (D-6). 

Controls punching of GRID, CTRIA2, PTRIA2 
and DTI cards from ALG. 

Appended static load vector (D-1, D-2). 
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P06V1 

06T 

t 

Displacement components used to plot deformed shape (o-6). 

PUNCH 

IC 

Output iredla request (PRINT or PUNCH). 

PURGE 

EH 

OMAP statement which causes conditional purging of data 
blocks. 

Purge 

PH 

A data block 1$ said to be purged when It Is flagged In the 
FIAT so that It will not be allocated to a physical file 
and so that modules attempting to access it will be signaled. 

PVECT 

DBH 

Partitioning vector for cyclic modes (A-9). 

PVISC 

IB 

Viscous element property definition card. 

PVT 

PH 

Parameter value table. The PVT contains BCD names and 
values of all parameters Input by means of PARAM bulk data 
cards. It Is generated by the preface module IFP and Is 
written on the P 'blem Tape. 

PI 

PU 

INPUTT2 rewind option. 

P2 

PU 

INPUTT2 unit number. 

P3 

PU 

INPUTT 2 tape ID. 

QB0Y1 

IB 

Defines uniform heat flux Into HBOY elements. 

QB0Y2 

IB 

Defines grid oolnt heat flux into HBDY elements. 

QBG 

DBM 

Single point forces of constraint in the Differential 
Stiffness Rigid Format (D-4). 

QDMEM 

IC 

Requests structure plot for all QDMEM elements. 

QDMEM1 

IC 

Requests structure plot for all QDMEM1 elements. 

Q0HEM2 

IC 

Requests structure plot for all QDMEM2 elements. 

QDPLT 

IC 

Requests structure plot for all QDPLT elements. 

QG 

DBM 

Constraint forces for all grid points. 

QHBOY 

IB 

Defines thermal load for steady-state heat conduction. 

QHHL 

DBML 

Aerodynamic matrix list * h-set. 

QJHt 

DBML 

Aerodynamic transformation matrix between h and j sets. 

QP 

DBM 

Constraint forces for all physical points. 

QPC 

DBM 

Complex single point forces of constraint for all physical 
points. 

OR 

DBM 

(q r ) - Determinant support forces. 

QS 

08M 

(q s ) * Single-point constraint forces. 

QUA01 

IC 

Requests structure plot for all QUAD1 elements. 

QUAD2 

IC 

Requests structure plot for all QUA02 elements. 

QVkCT 

IB 

Defines thermal vector flux from distant source. 

QVOL 

IB 

Defines volume heat generation. 
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SIGMA 

P 

Oqflnts Stefan-Boltzmann constant In boat transfer analysis. 

SIGN 

PU 

Controls the typa of static aerothermo- 
elastic analysis performed. 

SIL 

DOT 

Scalar Index List for all grid points. 

SI LA 

OOT 

Scalar Index List - Aerodynamics. 

SILO 

DOT 

Scalar Index List for all grid points and extra scalar 
points introduced for dynamic analysis. 

SILGA 

OBT 

Scalar Index List - Aerodynamic boxes only. 

SINE 

IC 

Conical shell request for sine set boundary conditions. 

SINGLE 

P 

No single-point constraints. 

SKIP BETWEEN FRAMES 

IC 

Request to insert blank frames on SC 4020 plotter for X-V 
plots. 

SKJ 

DBM 

Integration matrix. 

SKPMGG 

P 

Parameter used in statics to control execution of functional 
module SMA2. 

SL8DY 

IB 

Defines list of points on interface between axisymroetric fluid 
and radial slots. 

SL0AD 

IB 

Scalar point load definition. 

SLT 

OBT 

Static Loads Table. 

SMA1 

FMS 

Structural Matrix Assembler - phase 1 - generates stiffness 
matrix [K ] and structural damping matrix [K^]. 

SMA2 

FMS 

Structural Matrix Assembler - phase 2 - generates mass 
matrix [M l and viscous damping matrix [B ]. 

SMA3 

FMS 

Structural Matrix Assembler - phase 3 • add general element 

contributions to the stiffness matrix [«„„]. 

99 

SMP1 

FMS 

Structural Matrix Parti ti oner - part 1. 

SMP2 

FMS 

Structural Matrix Parti t loner - part 2. 

SMPYAD 

FMM 

Performs multiply-add matrix operation for up to five 
multiplications and one addition. 

S0L 

IA 

Specifies which rigid format solution is to be used when APP 
Is DISPLACEMENT. 

Solution Points 

PH 

Points used in the formulation of the general K system. 

SALVE 

FMM 

Solves a set of linear algebraic equations. 

SfRTI 

IC 

Output Is sorted by frequency or time and then by external 10. 

SART2 

IC 

Output is sorted by external ID and then by frequency or time. 

SORT 3 

M 

Output is sorted by individual item or component and then by 
frequency or time. 

SPC 

IB 

Single-point constraint and enforced deformation definition. 
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SRC 

IC 

Single-point constraint set selection. 

SPC$ 

N 

Indicates restart with change In sing1e-po<nt constraint set 
selection. 

SPCt 

IB 

Single-point constraint definition. 

SPCAOO 

IB 

Single-point constraint set combination definition. 

SPCAX 

IB 

Conical shell single-point constraint definition. 

SPCF 

IC 

Abbreviated form of SPCF0RCE. 

SPCFRRCE 

IC 

Single-point constraint force output request. (UN-2. 3, 4. 2) 

Spill 

PH 

Secondary storage devices are used because there Is Insufficient 
main storage to perform a matrix calculation or a data pro- 
cessing operation. 

SPLINE 

D8T 

Spllnlng Data Table. 

SPLINE1 

IB 

Oef nes surface spline. 

SPLINE2 

IB 

Defines beam spline. 

SPRINT 

IB 

Scalar point definition card. 

SSG1 

FMS 

Static Solution Generator • part 1. 

SSG2 

FMS 

Static Solution Generator • part 2. 

SSG3 

FMS 

Static Solution G*neratur • part 3. 

SSG4 

FMS 

Static Solution Generator • part 4. 

SSGHT 

FM4 

Solution generator for nonlinear heat transfer analysis. 

STATIC 

IC 

Requests deformed structure plot for problem In Static Analysis. 

STATICS 

IA 

Selects statics rigid format for heat transfer or structural 
analysis. 

STATICS 

P 

Parameter used In SDR2 to Indicate Static Analysis. 

STEADY STATE 

IA 

Selects rigid format for nonlinear static heat transfer 
analysis. 

STERE0SC0PIC 

1C 

Requests stereoscopic projections for structure plot. 

STREAML 

PU 

Controls the punching of STREAML1 ard 
STREAML2 cards from ALG. 

STREAML 1 

IB 

Gives blade streamline data. 

STREAML 2 

IB 

Gives blade streamline data. 

STRESS 

IC 

Element stress output request. (UM-2.3, 4.2) 

Structural Element 

PH 

One of the finite elements used to represent a part of a 
structure. 

SUBCASE 

IC 

Subcase definition. 

SUBC0M 

IC 

This subcase Is a linear conblnatlon of previous subcases. 

SU8SEQ 

IC 

Specifies coefficients for SUBC0M subcases. 

SUBTITLE 

IC 

Output labeling data for printer output. 

SUP AX 

IB 

Ficticious support for conical shell problem. 
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Z0RI6N 


P(J Aerodynamic Modification factor (C 16). 
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